Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,761,620 T→G 5.6% K641Q (AAG→CAG)  yfjK ← radiation resistance protein; DEAD/H helicase‑like protein; CP4‑57 putative defective prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,761,6200TG 5.6% 2192.3 / 86.5 610K641Q (AAG→CAG) yfjKradiation resistance protein; DEAD/H helicase‑like protein; CP4‑57 putative defective prophage
Reads supporting (aligned to +/- strand):  ref base T (315/260);  new base G (26/9);  total (341/269)
Fisher's exact test for biased strand distribution p-value = 3.40e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

Only 100 of 622 total aligned reads displayed.

GTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAAT  >  NC_000913/2761471‑2761769
                                                                                                                                                      |                                                                                                                                                     
gtcgGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCtt                                                                                                                                                       <  1:1330428/150‑1 (MQ=255)
    gAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACa                                                                                                                                                   >  2:3904001/1‑150 (MQ=255)
     aaaaTTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAg                                                                                                                                                  >  1:8183146/1‑150 (MQ=255)
        aTTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCgt                                                                                                                                               >  1:519857/1‑150 (MQ=255)
        aTTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCgt                                                                                                                                               >  2:3828445/1‑150 (MQ=255)
         ttCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCgtg                                                                                                                                              >  2:2938427/1‑150 (MQ=255)
           ccAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCgtggt                                                                                                                                            >  2:2068607/1‑150 (MQ=255)
           ccAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTGTAAAGCgtggt                                                                                                                                            <  2:5983409/150‑1 (MQ=255)
            cAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCgtggtg                                                                                                                                           >  1:2530446/1‑150 (MQ=255)
              gCAAGTCTACCAGTCGGTGTAACGTCTGATTGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTGTAAAGCgtggtata                                                                                                                                         <  1:2456319/150‑4 (MQ=255)
               cAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCGCGCTGTAAAGCgtggtatac                                                                                                                                        <  2:6904274/150‑5 (MQ=255)
               cAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGttt                                                                                                                                        <  1:7436111/150‑1 (MQ=255)
                 aGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTat                                                                                                                                      <  2:4518616/150‑2 (MQ=255)
                   tCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACt                                                                                                                                    >  1:5479193/1‑150 (MQ=255)
                      aCCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACa                                                                                                                                 >  1:1159254/1‑150 (MQ=255)
                      aCCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACa                                                                                                                                 >  1:5034119/1‑150 (MQ=255)
                       ccAGTCGGTGTAACGTCTGACTGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACAt                                                                                                                                <  2:647097/150‑1 (MQ=255)
                        cAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACAtc                                                                                                                               <  2:3617387/150‑1 (MQ=255)
                        cAGTCGGGGTAACGTCGGAATGGGAATGCCAATCTCGTCTAAGGCATTGACCCCTGGGGGTAGATGGAAACTTTCGACCGCCAAGGTGACATGGGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACAtc                                                                                                                               <  2:3369486/150‑1 (MQ=255)
                         aGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACAtcc                                                                                                                              >  2:240236/1‑150 (MQ=255)
                                  aaCGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAg                                                                                                                     <  2:1225407/150‑1 (MQ=255)
                                  aaCGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAg                                                                                                                     <  1:4888480/150‑1 (MQ=255)
                                      tCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGa                                                                                                                 <  2:7253476/150‑1 (MQ=255)
                                       cTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAg                                                                                                                <  1:3676652/150‑1 (MQ=255)
                                          aaTGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACACCCTCCATCAAGCTGAGCAc                                                                                                             <  1:7211035/150‑1 (MQ=255)
                                            tGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCt                                                                                                           >  1:3786599/1‑150 (MQ=255)
                                            tGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TATGCTGAGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCt                                                                                                           >  2:6169024/1‑150 (MQ=255)
                                                   gCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtag                                                                                                    <  1:4870584/150‑1 (MQ=255)
                                                    ccAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtagg                                                                                                   >  1:1016396/1‑150 (MQ=255)
                                                    ccAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtagg                                                                                                   <  1:3930142/150‑1 (MQ=255)
                                                     cAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAggg                                                                                                  >  1:8421495/1‑150 (MQ=255)
                                                     cAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAggg                                                                                                  <  2:302995/150‑1 (MQ=255)
                                                      aaTCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGt                                                                                                 >  1:6505056/1‑150 (MQ=255)
                                                       aTCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGta                                                                                                >  2:5591145/1‑150 (MQ=255)
                                                        tctcTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGtat                                                                                               >  1:2770957/1‑150 (MQ=255)
                                                              ctAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGa                                                                                         >  2:5515788/1‑150 (MQ=255)
                                                              ctAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGa                                                                                         <  1:6083254/150‑1 (MQ=255)
                                                               tAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGaa                                                                                        >  1:4097894/1‑150 (MQ=255)
                                                                 aGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCAACTTAGGTAGGGTATAGCCGAAg                                                                                       <  2:4997602/150‑1 (MQ=255)
                                                                     aTTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTAtt                                                                                  <  1:5914105/150‑1 (MQ=255)
                                                                     aTTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTAtt                                                                                  <  1:6915228/150‑1 (MQ=255)
                                                                      ttGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTg                                                                                 >  2:3279940/1‑150 (MQ=255)
                                                                      ttGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTg                                                                                 >  2:7360916/1‑150 (MQ=255)
                                                                      ttGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTg                                                                                 >  1:7646094/1‑150 (MQ=255)
                                                                      ttGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTg                                                                                 >  1:2616460/1‑150 (MQ=255)
                                                                       tGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTgg                                                                                >  1:5690428/1‑150 (MQ=255)
                                                                           acCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTGAGGTAGGGTATAGCCGAAGGTATTGGTTAt                                                                            >  2:1363386/1‑150 (MQ=255)
                                                                             cTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGa                                                                          <  2:4910244/150‑1 (MQ=255)
                                                                             cTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGa                                                                          <  1:8324957/150‑1 (MQ=255)
                                                                             cTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGa                                                                          <  1:1908186/150‑1 (MQ=255)
                                                                              tggtggTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAg                                                                         <  2:1771210/150‑1 (MQ=255)
                                                                                gtggtAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCt                                                                       >  2:6631475/1‑150 (MQ=255)
                                                                                   gtAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCa                                                                    >  1:463165/1‑150 (MQ=255)
                                                                                   gtAGATGGAAACTTTCGAACGCCAATTTGACATGTGGGTAATCGAC‑TGTGCTGCGTATACCACGCTGTACAGCGGGGTGTTTAACTACAGCCTCCATCAAGCTGAGCACCTGCGGGCGGGTCTAGCCGCAGGGCTTGGGTCTGAGCtgcc                                                                    >  2:2600948/1‑147 (MQ=255)
                                                                                       aTGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATc                                                                >  1:3052000/1‑150 (MQ=255)
                                                                                       aTGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATc                                                                >  1:4339788/1‑150 (MQ=255)
                                                                                        tGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCa                                                               >  1:7281175/1‑150 (MQ=255)
                                                                                        tGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTAAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCa                                                               >  1:3503861/1‑150 (MQ=255)
                                                                                         ggAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAt                                                              <  2:806186/150‑1 (MQ=255)
                                                                                          gAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAtt                                                             >  2:7561455/1‑150 (MQ=255)
                                                                                             aCTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAttgtt                                                          <  1:2204042/150‑1 (MQ=255)
                                                                                               tttCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGa                                                        <  2:5303829/150‑1 (MQ=255)
                                                                                               tttCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGa                                                        >  1:4781076/1‑150 (MQ=255)
                                                                                                ttCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAg                                                       <  2:1925267/150‑1 (MQ=255)
                                                                                                 tCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGAATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGATGAGCACCTTAGGTACGGTATAGCCGAAGGTATTGCTTATGAGCTTCAAATCATTGTTGAgg                                                      >  2:5165207/1‑150 (MQ=255)
                                                                                                               gACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGt                                        >  2:110966/1‑150 (MQ=255)
                                                                                                                aCATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTa                                       <  1:7915609/150‑1 (MQ=255)
                                                                                                                 cATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTAc                                      <  1:6268499/150‑1 (MQ=255)
                                                                                                                   tgtgtgTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGt                                    <  1:7016181/150‑1 (MQ=255)
                                                                                                                     tgtgtAATCGTC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGTTTATGATCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGttt                                  <  1:5822252/150‑1 (MQ=255)
                                                                                                                     tgtgtAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGttt                                  <  1:7972245/150‑1 (MQ=255)
                                                                                                                       tgtAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTtct                                <  2:2516947/150‑1 (MQ=255)
                                                                                                                         tAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTCTAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTcc                              >  2:8410717/1‑150 (MQ=255)
                                                                                                                          aaTCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTccc                             <  1:1584110/150‑1 (MQ=255)
                                                                                                                           aTCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTcccc                            >  1:7072917/1‑150 (MQ=255)
                                                                                                                            tCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCt                           <  2:718385/150‑1 (MQ=255)
                                                                                                                            tCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCt                           >  1:3389106/1‑150 (MQ=255)
                                                                                                                             cGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCtt                          <  1:7796835/150‑1 (MQ=255)
                                                                                                                              gAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTc                         >  1:4235818/1‑150 (MQ=255)
                                                                                                                               aC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCa                        <  2:1036478/150‑1 (MQ=255)
                                                                                                                                     gctgcGTATACCACGCTTTACAGCGTGGTGTTGAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGCAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGaact                   >  2:8286887/1‑148 (MQ=255)
                                                                                                                                       tgcGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATcc                 >  2:4087468/1‑150 (MQ=255)
                                                                                                                                        gcGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCa                >  2:6047107/1‑150 (MQ=255)
                                                                                                                                         taaacaCCACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCAc               >  2:4657555/6‑150 (MQ=255)
                                                                                                                                         cGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCAc               <  1:6389527/150‑1 (MQ=255)
                                                                                                                                         cGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCAc               <  1:6769997/150‑1 (MQ=255)
                                                                                                                                          gTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTCTTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACt              >  2:891867/1‑150 (MQ=255)
                                                                                                                                           tataCCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTg             >  2:6724141/1‑150 (MQ=255)
                                                                                                                                           aacaCCACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTg             >  1:4502284/4‑150 (MQ=255)
                                                                                                                                               ccACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTAGGGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGt         >  2:2785053/1‑150 (MQ=255)
                                                                                                                                                cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTg        >  2:5721696/1‑150 (MQ=255)
                                                                                                                                                cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTg        >  1:29007/1‑150 (MQ=255)
                                                                                                                                                cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTg        >  2:6834675/1‑150 (MQ=255)
                                                                                                                                                cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTg        >  1:1773280/1‑150 (MQ=255)
                                                                                                                                                cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTg        >  1:6684227/1‑150 (MQ=255)
                                                                                                                                                cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTg        >  2:1179044/1‑150 (MQ=255)
                                                                                                                                                cACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTg        >  2:4841663/1‑150 (MQ=255)
                                                                                                                                                 aCGCTTTACAGCGTGGTGTTTAACTACATCCTCCCTCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGc       <  2:1567330/150‑1 (MQ=255)
                                                                                                                                                      ttACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAAt  <  1:3169432/150‑1 (MQ=255)
                                                                                                                                                      ttACAGCGTGGTGTTTAAATACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTAAGTTTCTCACCTTCAGAAATCCACTGTCGTGCATAAt  >  2:7814712/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAC‑TTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAAT  >  NC_000913/2761471‑2761769

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑