Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,902,119 C→A 5.5% intergenic (‑271/+399) yobF ← / ← yebO DUF2527 family heat‑induced protein/putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,902,1190CA 5.5% 552.0 / 11.9 146intergenic (‑271/+399)yobF/yebODUF2527 family heat‑induced protein/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (44/94);  new base A (2/6);  total (46/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

Only 100 of 155 total aligned reads displayed.

TGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGAAAAATTGATGAATGAGCAA  >  CP009273/1901970‑1902268
                                                                                                                                                     |                                                                                                                                                     
tGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAgtgcgtgc                                                                                                                                                       >  1:4600086/1‑150 (MQ=255)
tGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAgtgcgtgc                                                                                                                                                       >  2:1297356/1‑150 (MQ=255)
tGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAgtgcgtgc                                                                                                                                                       >  1:2355537/1‑150 (MQ=255)
 gCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCa                                                                                                                                                      <  1:4442084/150‑1 (MQ=255)
 gCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCa                                                                                                                                                      >  2:3817796/1‑150 (MQ=255)
   acacTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAat                                                                                                                                                    <  1:423114/150‑1 (MQ=255)
   acacTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAat                                                                                                                                                    <  1:3712033/150‑1 (MQ=255)
   acacTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAat                                                                                                                                                    <  1:3711983/150‑1 (MQ=255)
           gCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACAtt                                                                                                                                            <  1:921576/150‑1 (MQ=255)
           gCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACAtt                                                                                                                                            <  1:1920788/150‑1 (MQ=255)
            cAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACAttt                                                                                                                                           <  2:603962/150‑1 (MQ=255)
            cAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACAttt                                                                                                                                           >  2:4270388/1‑150 (MQ=255)
            cAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACAttt                                                                                                                                           >  2:4680414/1‑150 (MQ=255)
            cAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACAttt                                                                                                                                           >  2:4680358/1‑150 (MQ=255)
              aCGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACAttttt                                                                                                                                         <  1:4541217/150‑1 (MQ=255)
                   aTACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTa                                                                                                                                    >  2:1202758/1‑150 (MQ=255)
                     aCGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATg                                                                                                                                  <  1:1286931/150‑1 (MQ=255)
                     aCGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGATCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGACGTTATTACATAGTGCGTGCCATATCACATTTTTTGCTATg                                                                                                                                  >  2:3886942/1‑150 (MQ=255)
                     aCGCGCAGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGAAATATCACATTTTTTGCTATg                                                                                                                                  <  2:3748530/150‑1 (MQ=255)
                      cgcgCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGc                                                                                                                                 >  2:2850290/1‑150 (MQ=255)
                      cgcgCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGc                                                                                                                                 >  2:458945/1‑150 (MQ=255)
                       gcgcCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCa                                                                                                                                <  2:407363/150‑1 (MQ=255)
                          ccGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATg                                                                                                                             <  2:458455/150‑1 (MQ=255)
                          ccGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATg                                                                                                                             <  2:1904867/150‑1 (MQ=255)
                          ccGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATg                                                                                                                             >  2:2371356/1‑150 (MQ=255)
                            gATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGaa                                                                                                                           <  1:3094511/150‑1 (MQ=255)
                                 atCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATaaaa                                                                                                                      <  1:1540825/150‑1 (MQ=255)
                                  tCATTGTGCTTAACCTTGCCAGTTCAGTCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTTTTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATaaaaa                                                                                                                     <  1:2489337/150‑1 (MQ=255)
                                  tCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATaaaaa                                                                                                                     <  1:3778136/150‑1 (MQ=255)
                                  tCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATaaaaa                                                                                                                     <  1:2489042/150‑1 (MQ=255)
                                    aTTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGt                                                                                                                   >  1:4016530/1‑150 (MQ=255)
                                     ttGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGtt                                                                                                                  <  2:2325845/150‑1 (MQ=255)
                                      tgtgCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTta                                                                                                                 <  2:4484433/150‑1 (MQ=255)
                                          cTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATc                                                                                                             >  2:1208959/1‑150 (MQ=255)
                                              aCCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACtt                                                                                                         <  1:1077080/150‑1 (MQ=255)
                                               ccTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACttt                                                                                                        <  2:14034/150‑1 (MQ=255)
                                                cTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGAAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACtttt                                                                                                       >  1:1063821/1‑150 (MQ=255)
                                                cTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGAAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACtttt                                                                                                       >  1:392461/1‑150 (MQ=255)
                                                  tGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTtc                                                                                                     <  1:4774881/150‑1 (MQ=255)
                                                  tGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTtc                                                                                                     <  1:1848030/150‑1 (MQ=255)
                                                  tGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTtc                                                                                                     >  2:1910646/1‑150 (MQ=255)
                                                     cAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCa                                                                                                  >  1:2605470/1‑150 (MQ=255)
                                                     cAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCa                                                                                                  >  2:3055579/1‑150 (MQ=255)
                                                               aGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAAcagtc                                                                                        <  2:4774297/150‑1 (MQ=255)
                                                                gATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAAcagtca                                                                                       <  2:1109652/150‑1 (MQ=255)
                                                                 aTACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGACTAAAAAGTTATATCACTTTTTCTCATAAAAcagtcag                                                                                      >  1:2750489/1‑150 (MQ=255)
                                                                 aTACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAAcagtcag                                                                                      <  2:4072204/150‑1 (MQ=255)
                                                                  tACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAAcagtcagt                                                                                     >  2:186665/1‑150 (MQ=255)
                                                                  tACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGAAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAAcagtcagt                                                                                     <  2:2674911/150‑1 (MQ=255)
                                                                    cTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTa                                                                                   >  1:3949409/1‑150 (MQ=255)
                                                                      tAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAAc                                                                                 <  1:3900881/150‑1 (MQ=255)
                                                                          cTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCt                                                                             >  1:4156761/1‑150 (MQ=255)
                                                                              cGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTatta                                                                         <  2:966048/150‑1 (MQ=255)
                                                                              cGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTatta                                                                         <  2:1238629/150‑1 (MQ=255)
                                                                               gTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTattaa                                                                        >  1:1768275/1‑150 (MQ=255)
                                                                               gTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTattaa                                                                        >  1:3281265/1‑150 (MQ=255)
                                                                               gTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTattaa                                                                        >  2:3666890/1‑150 (MQ=255)
                                                                               gTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTattaa                                                                        <  2:2846996/150‑1 (MQ=255)
                                                                                tATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTattaat                                                                       <  1:4326174/150‑1 (MQ=255)
                                                                                tATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTattaat                                                                       <  2:3709997/150‑1 (MQ=255)
                                                                                tATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTattaat                                                                       >  2:3599950/1‑150 (MQ=255)
                                                                                    ttCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTAcc                                                                   <  2:923998/150‑1 (MQ=255)
                                                                                     tCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTAccc                                                                  <  1:9225/150‑1 (MQ=255)
                                                                                      cTAACATAGTTCAATAAGTTCAGTTGTTCCATGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCt                                                                 <  2:4650510/150‑1 (MQ=255)
                                                                                      cTAACATAGTTCAATAAGTTCAGTTGCTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCt                                                                 <  2:2249972/150‑1 (MQ=255)
                                                                                            tAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAagag                                                           >  1:1786452/1‑150 (MQ=255)
                                                                                             aGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAgaga                                                          <  1:3158231/150‑1 (MQ=255)
                                                                                             aGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAgaga                                                          >  2:1458343/1‑150 (MQ=255)
                                                                                             aGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAgaga                                                          <  1:506397/150‑1 (MQ=255)
                                                                                             aGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAgaga                                                          >  1:1871076/1‑150 (MQ=255)
                                                                                             aGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAgaga                                                          <  1:2899731/150‑1 (MQ=255)
                                                                                               ttCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAg                                                        <  1:800270/150‑1 (MQ=255)
                                                                                               ttCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAg                                                        <  1:800114/150‑1 (MQ=255)
                                                                                                tCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGt                                                       <  2:391402/150‑1 (MQ=255)
                                                                                                   aTAAGTTCTGTTGTTCCAGGCCAGTTGGAGAAGTTATTACATAGTGAGTTCAATATCACATTTTTTGCTATTCAATGAATATAAATTTATATCACTTTTTCTCATAAAAAAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCaa                                                    <  1:1589269/150‑1 (MQ=255)
                                                                                                        ttCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATcccc                                               <  2:1768275/150‑1 (MQ=255)
                                                                                                          cAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGAAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCaa                                             <  2:392461/150‑1 (MQ=255)
                                                                                                           aGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCaaa                                            >  2:4474806/1‑150 (MQ=255)
                                                                                                           aGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCaaa                                            >  2:3103198/1‑150 (MQ=255)
                                                                                                            gttgttCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAg                                           <  1:2850290/150‑1 (MQ=255)
                                                                                                             ttgttCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAgg                                          <  1:4421083/150‑1 (MQ=255)
                                                                                                                ttCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGAt                                       <  1:1264756/150‑1 (MQ=255)
                                                                                                                ttCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGAt                                       <  2:2931621/150‑1 (MQ=255)
                                                                                                                ttCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGAt                                       <  1:2044843/150‑1 (MQ=255)
                                                                                                                 tCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGAtt                                      <  2:4409388/150‑1 (MQ=255)
                                                                                                                 tCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGAtt                                      <  1:3057827/150‑1 (MQ=255)
                                                                                                                     ggCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTaa                                  <  1:3486831/150‑1 (MQ=255)
                                                                                                                        cAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAAtt                               <  1:1910646/150‑1 (MQ=255)
                                                                                                                            gggAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTaaa                           <  1:4362288/150‑1 (MQ=255)
                                                                                                                              gagaAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAAt                         <  2:2614919/150‑1 (MQ=255)
                                                                                                                              gagaAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAAt                         <  2:2560297/150‑1 (MQ=255)
                                                                                                                               agaAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATa                        <  2:618870/150‑1 (MQ=255)
                                                                                                                                 aaGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAg                      <  2:3810351/150‑1 (MQ=255)
                                                                                                                                        taCATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGAAAAAtt               <  2:107202/150‑1 (MQ=255)
                                                                                                                                          cATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGAAAAATtga             <  2:2355537/150‑1 (MQ=255)
                                                                                                                                             aGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGAAAAATtgatga          <  2:2528932/150‑1 (MQ=255)
                                                                                                                                             aGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGAAAAATtgatga          <  2:4637750/150‑1 (MQ=255)
                                                                                                                                             aGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGAAAAATtgatga          >  1:691942/1‑150 (MQ=255)
                                                                                                                                                 cgtgcAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGAAAAATTGatgaatga      <  2:882484/150‑1 (MQ=255)
                                                                                                                                                     cAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGAAAAATTGATGAATGAGCaa  >  1:18646/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTGCAATATCACATTTTTTGCTATGCAATGAATAAAAAGTTATATCACTTTTTCTCATAAAACAGTCAGTTAACGGCTATTAATTACCCTAAAGAGAAGTCAATCCCCAAAGGGATTGTAAATTTAAAATAAGAAAAATTGATGAATGAGCAA  >  CP009273/1901970‑1902268

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: