Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,725,887 A→C 5.3% L773L (CTT→CTG plsB ← glycerol‑3‑phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,725,8870AC 5.3% 334.7 / 7.2 112L773L (CTT→CTGplsBglycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (52/54);  new base C (3/3);  total (55/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 116 total aligned reads displayed.

AGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTGCCCGCTGTTGATCA  >  NC_002947/1725738‑1726035
                                                                                                                                                      |                                                                                                                                                    
aGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGATTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGa                                                                                                                                                      >  1:275223/1‑150 (MQ=255)
 gCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGaa                                                                                                                                                     >  1:1735838/1‑150 (MQ=255)
 gCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCATCGGCCAGCTAGACGGGCTGGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGaa                                                                                                                                                     <  1:1487998/150‑1 (MQ=255)
   ccGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGc                                                                                                                                                   >  1:1657644/1‑150 (MQ=255)
     gtcgtcACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCaa                                                                                                                                                 <  1:995449/150‑1 (MQ=255)
        gtcACGGTGTAGGGACCCCTGGGGACTCATTGGGCGCAGCTCGGCCCACAGTACTCGCTTGGCCCGAACGCAGTCCACCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGt                                                                                                                                              <  1:645118/150‑1 (MQ=255)
        gtcACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTTCAGCAAggc                                                                                                                                              <  1:1530024/150‑2 (MQ=255)
        gtcACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGt                                                                                                                                              >  2:354004/1‑150 (MQ=255)
         tcACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTTCAGCAAggcg                                                                                                                                             <  2:1002801/150‑3 (MQ=255)
             ggTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTGCTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACCCCGCTTGCCTTGCGCGTCGGCGTGCAGCACTCCTTGCTGAAGCAAGGTCtgaa                                                                                                                                         <  1:300681/150‑1 (MQ=255)
                  aGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAg                                                                                                                                    >  2:1280915/1‑150 (MQ=255)
                   gggCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGt                                                                                                                                   <  1:1746915/150‑1 (MQ=255)
                    ggCCACCTGGCGGATCGACAGGCGCAGCTCGGACGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTg                                                                                                                                  >  2:2345955/1‑150 (MQ=255)
                      ccACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGc                                                                                                                                <  2:2167579/150‑1 (MQ=255)
                               ggATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCg                                                                                                                       >  1:1608517/1‑150 (MQ=255)
                                gATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTTCAGCAAggc                                                                                                                                              <  1:2235146/126‑2 (MQ=255)
                                   cGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCtt                                                                                                                   <  1:1804061/150‑1 (MQ=255)
                                     aCAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAAc                                                                                                                           >  1:607696/1‑140 (MQ=255)
                                     aCAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAAc                                                                                                                           <  2:607696/140‑1 (MQ=255)
                                     aCAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGt                                                                                                                 >  1:763651/1‑150 (MQ=255)
                                     aCAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGACAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGt                                                                                                                 >  2:2424526/1‑150 (MQ=255)
                                      cAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATACCTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTc                                                                                                                <  1:2029265/150‑1 (MQ=255)
                                      cAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTc                                                                                                                <  1:793167/150‑1 (MQ=255)
                                        ggCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCGCGGCCAGCTCGGCGAGCTTGGCGTGATAACTCGGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACGGCGTCTTGTCga                                                                                                              >  2:1329364/1‑150 (MQ=255)
                                         gcgcAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCCGAATGAAGTGGCGGAACAGCGTCTTGTCgaa                                                                                                             <  1:752414/150‑1 (MQ=255)
                                          cgcAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCgaag                                                                                                            <  2:367901/150‑1 (MQ=255)
                                           gcAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGc                                                                                                                        <  2:637673/137‑1 (MQ=255)
                                           gcAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGc                                                                                                                        >  1:637673/1‑137 (MQ=255)
                                           gcAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCgaaga                                                                                                           >  2:1711150/1‑150 (MQ=255)
                                             aGCGCGGACCACAGGATTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAAc                                                                                                         <  2:735093/150‑1 (MQ=255)
                                                tCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGc                                                                                                                                >  1:1185881/1‑124 (MQ=255)
                                                tCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTCAAGCAAGGTCTGAATGAAGTGGc                                                                                                                                <  2:1185881/124‑1 (MQ=255)
                                                tCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTCGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTcc                                                                                                      <  2:1657644/150‑1 (MQ=255)
                                                  ggCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCtgaatg                                                                                                                                       <  1:1148184/115‑1 (MQ=255)
                                                  ggCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCtgaatg                                                                                                                                       >  2:1148184/1‑115 (MQ=255)
                                                  ggCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCtgaat                                                                                                                                        >  2:2211858/1‑114 (MQ=255)
                                                  ggCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCtgaat                                                                                                                                        <  1:2211858/114‑1 (MQ=255)
                                                   tgctaCAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGt                                                                                                   <  2:502139/146‑1 (MQ=255)
                                                   gCCGACAGTACTCGCTTGGCCACGCCCTCGGTCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGt                                                                                                   <  1:852755/150‑1 (MQ=255)
                                                   gCCGACAGTACTCGCTTGGCCACGCCCTCGGACAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGt                                                                                                   <  2:1248231/150‑1 (MQ=255)
                                                             cTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCga                                                                                                              <  1:270715/129‑1 (MQ=255)
                                                             cTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGc                                                                                         >  1:1997452/1‑150 (MQ=255)
                                                             cTCGCTTGGCCACGCCCTCGGCCAGCTCGCAGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGt                                                                                                   <  2:1238372/140‑1 (MQ=255)
                                                             cTCGCTTGGCCACGCCCTCGGCCAGCTCGCAGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGt                                                                                                   >  1:1238372/1‑140 (MQ=255)
                                                              ttgCTTGGCCACGCCCGCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGcc                                                                                        >  2:1215702/3‑150 (MQ=255)
                                                              tCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGt                                                                                                                 >  1:2194516/1‑125 (MQ=255)
                                                              tCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGt                                                                                                                 <  2:2194516/125‑1 (MQ=255)
                                                              tCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGcc                                                                                        >  1:2383681/1‑150 (MQ=255)
                                                                gCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCCGTCAGCCTTGCGCGTCGGCGTGCAGTACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCAt                                                                                      >  1:1226247/1‑150 (MQ=255)
                                                                gCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGc                                                                                                                        >  2:2408264/1‑116 (MQ=255)
                                                                gCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGc                                                                                                                        <  1:2408264/116‑1 (MQ=255)
                                                                  ttGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGc                                                                                    >  1:737918/1‑150 (MQ=255)
                                                                   tGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCa                                                                                   <  2:275223/150‑1 (MQ=255)
                                                                     gCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAgg                                                                                 <  1:1711150/150‑1 (MQ=255)
                                                                      ccACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGa                                                                                >  2:2284417/1‑150 (MQ=255)
                                                                      ccACGCCCTCGGCCAGCTCGCCGAGCTTGACGTGATAAGTCAGCTTGCCTTGCGCGTGGGGGTGCAGCACGCGTTGATGCTGCGCCGTACGAATTATGGGTCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGATCAGTCCATTCAGGa                                                                                >  1:542986/1‑150 (MQ=255)
                                                                        aCGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATc                                                                              <  2:737918/150‑1 (MQ=255)
                                                                          gCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGa                                                                            <  2:1608517/150‑1 (MQ=255)
                                                                             cTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAg                                                                         >  2:1961262/1‑150 (MQ=255)
                                                                              tCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAgg                                                                        >  2:2437806/1‑150 (MQ=255)
                                                                                 ggCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGgcg                                                                      >  1:2448723/1‑150 (MQ=255)
                                                                                   ccAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCt                                                                    >  2:59219/1‑150 (MQ=255)
                                                                                   ccAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCATTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCt                                                                    >  1:1993978/1‑150 (MQ=255)
                                                                                    cAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTg                                                                   >  1:862673/1‑150 (MQ=255)
                                                                                    cAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTg                                                                   <  2:542986/150‑1 (MQ=255)
                                                                                    cAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTg                                                                   >  2:2145721/1‑150 (MQ=255)
                                                                                     aGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTgg                                                                  <  1:2437806/150‑1 (MQ=255)
                                                                                          gCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCcat                                                             >  2:1490597/1‑150 (MQ=255)
                                                                                          gCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCcat                                                             >  2:2075515/1‑150 (MQ=255)
                                                                                             gAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCcatcat                                                          >  2:1074871/1‑150 (MQ=255)
                                                                                                cTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCAc                                                       >  2:234961/1‑150 (MQ=255)
                                                                                                cTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCAc                                                       <  1:975979/150‑1 (MQ=255)
                                                                                                 ttGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACg                                                      >  1:1304333/1‑150 (MQ=255)
                                                                                                   gTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTg                                                                                                                  >  2:519824/1‑88 (MQ=255)
                                                                                                   gTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTg                                                                                                                  <  1:519824/88‑1 (MQ=255)
                                                                                                    tCGTGATAACTCAGCTTGCCTTGCGCATCGGCGTGCAGCACGCCTTGCTGAAGCAAg                                                                                                                                                >  2:1834337/1‑57 (MQ=255)
                                                                                                     cGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCAAGcaca                                                  <  2:595227/150‑1 (MQ=255)
                                                                                                      gTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAg                                                 <  2:1801092/150‑1 (MQ=255)
                                                                                                      gTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAg                                                 <  1:2424526/150‑1 (MQ=255)
                                                                                                          tAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCt                                             >  1:1179970/1‑150 (MQ=255)
                                                                                                           aaCTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCtt                                            <  1:2345955/150‑1 (MQ=255)
                                                                                                            aCTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTc                                           <  1:1879051/150‑1 (MQ=255)
                                                                                                            aCTCAGCTTGCATTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTc                                           <  2:1993978/150‑1 (MQ=255)
                                                                                                              tCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTg                                                                                                                  <  2:2034884/77‑1 (MQ=255)
                                                                                                              tCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTg                                                                                                                  >  1:2034884/1‑77 (MQ=255)
                                                                                                              tCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGGCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCa                                         >  1:1402837/1‑150 (MQ=255)
                                                                                                                  cttgccttgcGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGt                                                                                                   <  2:207021/88‑1 (MQ=255)
                                                                                                                       cttgcGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCgg                                >  1:1177821/1‑150 (MQ=255)
                                                                                                                        ttgcGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTggg                                                                 >  1:1995178/1‑116 (MQ=255)
                                                                                                                        ttgcGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTggg                                                                 <  2:1995178/116‑1 (MQ=255)
                                                                                                                                 ggCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCt                      <  1:1215702/150‑1 (MQ=255)
                                                                                                                                    gTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGtt                   <  2:1179970/150‑1 (MQ=255)
                                                                                                                                        ttcagaccTTGCTTCAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTGc               >  1:2361860/7‑150 (MQ=255)
                                                                                                                                         gcaNGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTGcc              >  2:1498093/1‑150 (MQ=255)
                                                                                                                                         gcaCGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATc                                                                              <  1:1466677/86‑1 (MQ=255)
                                                                                                                                          caCGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTGccc             >  2:1765005/1‑150 (MQ=255)
                                                                                                                                                ttGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTGCCCGCTGtt       <  2:1304333/150‑1 (MQ=255)
                                                                                                                                                   cTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTGCCCGCTGTTGAt    <  1:1214557/150‑1 (MQ=255)
                                                                                                                                                   cTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGGCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTGCCCGCTGTTGAt    <  1:375674/150‑1 (MQ=255)
                                                                                                                                                     gAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTGCCCGCTGTTGATca  <  2:1818971/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
AGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTGCCCGCTGTTGATCA  >  NC_002947/1725738‑1726035

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