Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,674,178 A→C 8.4% E89A (GAG→GCG) ‡ tqsA → pheromone AI‑2 transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,674,1780AC 8.4% 242.2 / 11.5 83E89A (GAG→GCG) ‡tqsApheromone AI‑2 transporter
Reads supporting (aligned to +/- strand):  ref base A (42/34);  new base C (3/4);  total (45/38)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

CGCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTC  >  NC_000913/1674029‑1674323
                                                                                                                                                     |                                                                                                                                                  
cGCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGa                                                                                                                                                    >  2:1501154/1‑150 (MQ=255)
 gCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAg                                                                                                                                                   >  2:1045362/1‑150 (MQ=255)
 gCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAg                                                                                                                                                   <  1:1433058/150‑1 (MQ=255)
  cATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGt                                                                                                                                                  >  2:9344/1‑150 (MQ=255)
      atttattGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGAc                                                                                                                                             >  1:867492/1‑150 (MQ=255)
         tattGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGTCAACTCGTTGagcgcg                                                                                                                                          <  2:1414187/150‑6 (MQ=255)
          attGCTGTTATTCTTAACCCGCTGGGGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCgg                                                                                                                                         <  2:1606302/150‑1 (MQ=255)
           ttGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGTCAACTCGTTGAGCGCGGa                                                                                                                                        <  1:1508389/150‑1 (MQ=255)
            tGCTGTTATTCTTAACCCGCTGGTGCAACACTTGGGCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCTATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGTCAACTCGTTGAGCGCGGaa                                                                                                                                       <  1:394990/150‑2 (MQ=255)
               tGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGtt                                                                                                                                    >  1:210950/1‑150 (MQ=255)
                gTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTa                                                                                                                                   >  1:1275332/1‑150 (MQ=255)
                gTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTa                                                                                                                                   >  1:1702261/1‑150 (MQ=255)
                 ttATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGTCAACTCGTTGAGCGCGGaacccag                                                                                                                                  <  2:1390103/150‑7 (MQ=255)
                  tATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTAcc                                                                                                                                 <  1:963923/150‑1 (MQ=255)
                   aTTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCg                                                                                                                                <  1:703945/150‑1 (MQ=255)
                   aTTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCg                                                                                                                                >  1:90649/1‑150 (MQ=255)
                     tCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCa                                                                                                                              >  2:409952/1‑150 (MQ=255)
                             cGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGAAGTTACCGCAATATCGCa                                                                                                                      >  1:1220800/1‑150 (MQ=255)
                              gCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCaa                                                                                                                     >  2:75255/1‑150 (MQ=255)
                                   tGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCta                                                                                                                >  2:515605/1‑150 (MQ=255)
                                     cAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCACCTCtatt                                                                                                              >  2:835308/1‑150 (MQ=255)
                                     cAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCtatt                                                                                                              <  1:515605/150‑1 (MQ=255)
                                      aaCACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCGTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCtattt                                                                                                             >  1:300786/1‑149 (MQ=255)
                                            tGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGAcgc                                                                                                       >  1:1707088/1‑150 (MQ=255)
                                                  gCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGc                                                                                                 >  1:1179386/1‑150 (MQ=255)
                                                     ggCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAg                                                                                              <  2:500157/150‑1 (MQ=255)
                                                         tgtgCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCt                                                                                          <  1:75255/150‑1 (MQ=255)
                                                            gCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGa                                                                                                                                        <  1:720062/101‑1 (MQ=255)
                                                            gCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGa                                                                                                                                        >  2:720062/1‑101 (MQ=255)
                                                              cgcgTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAAc                                                                                     <  1:309792/150‑1 (MQ=255)
                                                              cgcgTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTGCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAAc                                                                                     <  1:502148/150‑1 (MQ=255)
                                                                  ggtACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgtt                                                                                 <  1:474093/149‑1 (MQ=255)
                                                                    tACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttgt                                                                               >  1:401078/1‑150 (MQ=255)
                                                                        ggCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCa                                                                           <  2:1417947/150‑1 (MQ=255)
                                                                        ggCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCa                                                                           <  1:409952/150‑1 (MQ=255)
                                                                           ggTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAAcg                                                                        >  2:1197835/1‑150 (MQ=255)
                                                                               tCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTa                                                                    >  1:275601/1‑150 (MQ=255)
                                                                                cGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAg                                                                   >  1:755130/1‑150 (MQ=255)
                                                                                 gATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAgg                                                                  <  2:867492/150‑1 (MQ=255)
                                                                                  aTTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAggg                                                                 >  1:1018171/1‑150 (MQ=255)
                                                                                   ttttGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGa                                                                >  1:291919/1‑150 (MQ=255)
                                                                                       gatgaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGa                                                            >  2:956532/1‑150 (MQ=255)
                                                                                       gatgaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGa                                                            <  1:1045362/150‑1 (MQ=255)
                                                                                        atgaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGAc                                                           <  2:291919/150‑1 (MQ=255)
                                                                                        atgaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGAc                                                           <  1:33962/150‑1 (MQ=255)
                                                                                           aCCATCATCGTGATGGCGATGGTGCTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGtc                                                        >  1:635171/1‑150 (MQ=255)
                                                                                            ccATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCACAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGtct                                                       >  2:1118486/1‑150 (MQ=255)
                                                                                             catcatCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGtctc                                                      <  2:755130/150‑1 (MQ=255)
                                                                                               tcatcGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAg                                                    <  1:928798/150‑1 (MQ=255)
                                                                                               tcatcGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAg                                                    <  1:934993/150‑1 (MQ=255)
                                                                                                catcGTGTTGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCTCGCTCAACGAGTTGACCCGGACGTTACCGCATTATCGCAACTCTATTATGTCGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGt                                                   <  1:881288/150‑1 (MQ=255)
                                                                                                    gTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGAc                                               >  1:765608/1‑150 (MQ=255)
                                                                                                    gTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGAc                                               <  2:1275332/150‑1 (MQ=255)
                                                                                                     tGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGAcc                                              <  1:1501154/150‑1 (MQ=255)
                                                                                                       aTGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCAGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGCCGTCTCAGTTGACCAg                                            >  2:207335/1‑150 (MQ=255)
                                                                                                           cGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTgg                                        <  1:9344/150‑1 (MQ=255)
                                                                                                           cGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTgg                                        >  1:766587/1‑150 (MQ=255)
                                                                                                            gATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGgc                                       <  2:90649/150‑1 (MQ=255)
                                                                                                              tGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGgcgc                                     <  1:1672647/150‑1 (MQ=255)
                                                                                                                gtgtTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCAt                                   <  2:1702261/150‑1 (MQ=255)
                                                                                                                 tgtTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCAtt                                  >  1:1697439/1‑150 (MQ=255)
                                                                                                                 tgtTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCAtt                                  >  2:177417/1‑150 (MQ=255)
                                                                                                                    tGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATtata                               >  2:1614907/1‑150 (MQ=255)
                                                                                                                     gCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCATCTATAATATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGATGGCGCATTatat                              >  1:297221/1‑150 (MQ=255)
                                                                                                                      cTATTAGCTTATCTGTGTTGCGCGCTCTACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTTTTATGCCGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACGAGCTGGCGCATTAtttt                             <  1:993341/150‑1 (MQ=255)
                                                                                                                       tATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTg                            >  1:1057415/1‑150 (MQ=255)
                                                                                                                           aGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATcc                        <  1:956532/150‑1 (MQ=255)
                                                                                                                            gCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCg                       >  2:947317/1‑150 (MQ=255)
                                                                                                                                 tCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAAcgc                  >  2:461283/1‑150 (MQ=255)
                                                                                                                                 ggtaacgtCCGCGTCAACTCGTTGAGCGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCACCGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAAcgc                  >  2:193347/8‑150 (MQ=255)
                                                                                                                                  gtaacgtCCGCGTCAACTCGTTGAGCGCGGACGTTACCGCAATATCGCCACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGCCGTCTCAGTTGCCCAGCTGGCGCATTATATTGATCCGCTcgcg                 >  1:1642687/7‑150 (MQ=255)
                                                                                                                                    gggTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACgcggc               <  2:1707088/150‑1 (MQ=255)
                                                                                                                                    gggTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACgcggc               <  2:1244514/150‑1 (MQ=255)
                                                                                                                                      cgtCCGCGTCAACTCGTTGAGCGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGa             >  1:541453/3‑150 (MQ=255)
                                                                                                                                       ttCCGCGCTCAACGAGTTGACGCGGACGTTACCACAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGAt            <  1:1118486/150‑1 (MQ=255)
                                                                                                                                        tCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATg           <  1:1614907/150‑1 (MQ=255)
                                                                                                                                        tCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATg           >  1:204870/1‑150 (MQ=255)
                                                                                                                                        tCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATg           <  1:1197835/150‑1 (MQ=255)
                                                                                                                                          cgcgCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGAc         >  1:1712949/1‑150 (MQ=255)
                                                                                                                                              cTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTg     >  1:686328/1‑150 (MQ=255)
                                                                                                                                                cAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCt   >  2:1374491/1‑150 (MQ=255)
                                                                                                                                                cAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCt   >  2:1200532/1‑150 (MQ=255)
                                                                                                                                                 aaCGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTc  >  2:711217/1‑150 (MQ=255)
                                                                                                                                                 aaCGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTAGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTc  >  1:1707256/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTC  >  NC_000913/1674029‑1674323

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