Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,730,144 A→T 6.6% A170A (GCA→GCT)  sodB → superoxide dismutase, Fe

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,730,1440AT 6.6% 658.3 / 28.7 213A170A (GCA→GCT) sodBsuperoxide dismutase, Fe
Reads supporting (aligned to +/- strand):  ref base A (115/84);  new base T (11/3);  total (126/87)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 213 total aligned reads displayed.

GCTGGACCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGG  >  CP009273/1729999‑1730290
                                                                                                                                                 |                                                                                                                                                  
gctggACCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACg                                                                                                                                                  <  1:1255945/148‑1 (MQ=255)
  tggACCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCTCGTc                                                                                                                                                <  1:1134636/148‑1 (MQ=255)
  tggACCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTc                                                                                                                                                <  1:65471/148‑1 (MQ=255)
  tggACCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTc                                                                                                                                                <  1:2374259/148‑1 (MQ=255)
   ggACCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTcc                                                                                                                                               <  1:2497861/148‑1 (MQ=255)
    gACCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCTCGTCCt                                                                                                                                              >  1:1158555/1‑148 (MQ=255)
        tggctggTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCt                                                                                                                                          >  1:1039819/1‑148 (MQ=255)
        tggctggTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCt                                                                                                                                          <  2:3380388/148‑1 (MQ=255)
         ggctggTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTa                                                                                                                                         >  2:114525/1‑148 (MQ=255)
          gctggTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCTCGTCCTGGCTAt                                                                                                                                        >  1:2612902/1‑148 (MQ=255)
          gctggTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATAGACTATCGCAATGCACGTCCTGGCTAt                                                                                                                                        >  2:2061951/1‑148 (MQ=255)
               tGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGa                                                                                                                                   >  1:3499544/1‑148 (MQ=255)
                  aaaaCAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCa                                                                                                                                <  1:2452211/148‑1 (MQ=255)
                   aaaCAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCAc                                                                                                                               >  1:1365758/1‑148 (MQ=255)
                     aCAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACtt                                                                                                                             <  1:2897897/148‑1 (MQ=255)
                       aGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCTCGTCCTGGCTATCTGGAGCACTTCt                                                                                                                           <  2:185382/148‑1 (MQ=255)
                         cGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTgg                                                                                                                         >  2:2950270/1‑148 (MQ=255)
                                 aaCTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGt                                                                                                                 >  1:996592/1‑148 (MQ=255)
                                 aaCTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGt                                                                                                                 >  2:2402606/1‑148 (MQ=255)
                                  aCTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTg                                                                                                                >  2:1820156/1‑148 (MQ=255)
                                  aCTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTg                                                                                                                >  1:2259391/1‑148 (MQ=255)
                                   cTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGCTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGa                                                                                                               >  1:3245673/1‑148 (MQ=255)
                                    tGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGaa                                                                                                              >  1:1108493/1‑148 (MQ=255)
                                      gCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACt                                                                                                            >  2:2649048/1‑148 (MQ=255)
                                         aTCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTggg                                                                                                         >  1:1009959/1‑148 (MQ=255)
                                           cGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCTCGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGaa                                                                                                       >  1:3034286/1‑148 (MQ=255)
                                           cGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGaa                                                                                                       <  2:2484828/148‑1 (MQ=255)
                                           cGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATAGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGaa                                                                                                       <  2:1249491/148‑1 (MQ=255)
                                            gTTTCAACCTCTAACGCGGCTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAAt                                                                                                      >  1:2328135/1‑148 (MQ=255)
                                                cAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGt                                                                                                  <  1:553561/148‑1 (MQ=255)
                                                  aCCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAg                                                                                                >  1:3259291/1‑148 (MQ=255)
                                                   ccTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGc                                                                                               >  1:3537609/1‑148 (MQ=255)
                                                   ccTCTAACGCGGGTACTCCGATGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTGCTGGGCGCTGGTGAACTGGGAATTCGTAGc                                                                                               >  1:196515/1‑148 (MQ=255)
                                                      ctAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGaa                                                                                            >  2:2687255/1‑148 (MQ=255)
                                                        aaCGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGaaaa                                                                                          <  1:1939092/148‑1 (MQ=255)
                                                        aaCGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGaaaa                                                                                          <  1:801150/148‑1 (MQ=255)
                                                          cgcgGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAAt                                                                                        <  2:589370/148‑1 (MQ=255)
                                                           gcgGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAAtc                                                                                       <  1:414393/148‑1 (MQ=255)
                                                            cgGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAAtct                                                                                      <  2:1773936/148‑1 (MQ=255)
                                                            cgGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAAtct                                                                                      <  2:1475588/148‑1 (MQ=255)
                                                            cgGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAAtct                                                                                      <  2:2774132/148‑1 (MQ=255)
                                                               gTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCTCGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGc                                                                                   >  2:539572/1‑148 (MQ=255)
                                                                tACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCt                                                                                  >  2:405146/1‑148 (MQ=255)
                                                                 aCTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTg                                                                                 <  2:784014/148‑1 (MQ=255)
                                                                  cTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGc                                                                                >  2:1929154/1‑148 (MQ=255)
                                                                  cTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGc                                                                                <  2:814670/148‑1 (MQ=255)
                                                                      gCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCataa                                                                            >  2:3489710/1‑148 (MQ=255)
                                                                        tGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCataata                                                                          <  1:1655071/148‑1 (MQ=255)
                                                                         gACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCAtaataa                                                                         >  2:1438936/1‑148 (MQ=255)
                                                                           ccaccGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACt                                                                       <  2:2455684/148‑1 (MQ=255)
                                                                           ccaccGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACt                                                                       <  1:1323157/148‑1 (MQ=255)
                                                                             accGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGa                                                                     >  1:811508/1‑148 (MQ=255)
                                                                                gATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATgg                                                                  <  2:1108605/148‑1 (MQ=255)
                                                                                 aTGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGc                                                                 >  1:2592228/1‑148 (MQ=255)
                                                                                   gCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCaa                                                               <  2:3314643/148‑1 (MQ=255)
                                                                                   gCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCaa                                                               <  2:327437/148‑1 (MQ=255)
                                                                                   gCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCaa                                                               <  2:1331774/148‑1 (MQ=255)
                                                                                   gCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCaa                                                               >  1:3397892/1‑148 (MQ=255)
                                                                                      aCTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATg                                                            >  2:2185634/1‑148 (MQ=255)
                                                                                            ctgctgACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCAt                                                      <  1:2646369/148‑1 (MQ=255)
                                                                                                tgACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGcc                                                  <  1:252662/148‑1 (MQ=255)
                                                                                                 gACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCt                                                 >  1:2784464/1‑148 (MQ=255)
                                                                                                    cGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGaa                                              <  2:2540861/148‑1 (MQ=255)
                                                                                                      ttGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAgc                                            >  1:3273286/1‑148 (MQ=255)
                                                                                                      ttGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAgc                                            >  2:1680742/1‑148 (MQ=255)
                                                                                                      ttGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAgc                                            <  1:3202590/148‑1 (MQ=255)
                                                                                                       tGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAgcg                                           <  2:3432513/148‑1 (MQ=255)
                                                                                                        gATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAgcgc                                          <  2:1230284/148‑1 (MQ=255)
                                                                                                         aTGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCt                                         >  2:2792654/1‑148 (MQ=255)
                                                                                                         aTGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCt                                         >  1:3364265/1‑148 (MQ=255)
                                                                                                          tGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTa                                        <  2:3160274/148‑1 (MQ=255)
                                                                                                           gTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTAc                                       >  1:775740/1‑148 (MQ=255)
                                                                                                              tGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCt                                    <  1:139993/148‑1 (MQ=255)
                                                                                                               gggAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCtt                                   <  1:3073391/148‑1 (MQ=255)
                                                                                                               gggAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCtt                                   <  1:3489710/148‑1 (MQ=255)
                                                                                                               gggAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAACATTGCCTGAAGCGCTACGCtt                                   >  2:2803300/1‑148 (MQ=255)
                                                                                                                ggAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTa                                  <  2:3004375/148‑1 (MQ=255)
                                                                                                                 gAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTAt                                 >  1:3173228/1‑148 (MQ=255)
                                                                                                                      cGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGc                            >  2:291549/1‑148 (MQ=255)
                                                                                                                        cTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCt                          <  1:2224964/148‑1 (MQ=255)
                                                                                                                        cTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCt                          >  1:3463574/1‑148 (MQ=255)
                                                                                                                         ttattaCATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTa                         >  2:1039259/1‑148 (MQ=255)
                                                                                                                          tattaCATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTAc                        >  1:3098034/1‑148 (MQ=255)
                                                                                                                            ttaCATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCACATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTAcgc                      >  2:2722627/1‑148 (MQ=255)
                                                                                                                             taCATCGACTATCGCAATGCTCGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTAcgcg                     >  2:1163787/1‑148 (MQ=255)
                                                                                                                               cATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGa                   <  2:1935072/148‑1 (MQ=255)
                                                                                                                               cATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGa                   <  2:1306226/148‑1 (MQ=255)
                                                                                                                               cATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGa                   >  1:2308885/1‑148 (MQ=255)
                                                                                                                                aTCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGat                  >  2:1564792/1‑148 (MQ=255)
                                                                                                                                  cGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGTGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGatca                <  2:3187425/148‑1 (MQ=255)
                                                                                                                                    aCTATCGCAATGCTCGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGatcatc              >  1:464317/1‑148 (MQ=255)
                                                                                                                                       aTCGCAATGCTCGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGAt           <  2:2574597/148‑1 (MQ=255)
                                                                                                                                       aTCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGAt           <  2:3343922/148‑1 (MQ=255)
                                                                                                                                       aTCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGAt           >  1:1869679/1‑148 (MQ=255)
                                                                                                                                        tCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATg          <  1:2950270/148‑1 (MQ=255)
                                                                                                                                        tCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATg          >  2:662937/1‑148 (MQ=255)
                                                                                                                                           cAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAg       >  1:2410994/1‑148 (MQ=255)
                                                                                                                                             aTGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCCTATCAGGCCTACGCGGATCATCGATGTAggg     >  1:2028978/1‑147 (MQ=255)
                                                                                                                                               gCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCg   >  1:638137/1‑148 (MQ=255)
                                                                                                                                                cACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCgg  <  2:2328135/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCTGGACCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGG  >  CP009273/1729999‑1730290

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