Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 394,474 A→C 5.1% intergenic (+56/‑32) yaiT → / → iprA putative autotransporter YaiT/putative DNA‑binding transcriptional regulator IprA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913394,4740AC 5.1% 488.8 / 12.7 154intergenic (+56/‑32)yaiT/iprAputative autotransporter YaiT/putative DNA‑binding transcriptional regulator IprA
Reads supporting (aligned to +/- strand):  ref base A (68/77);  new base C (4/4);  total (72/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 159 total aligned reads displayed.

GAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGATGTC  >  NC_000913/394325‑394622
                                                                                                                                                     |                                                                                                                                                    
caaCTTCAGAGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCAGCGAACTTGGGTGTTAAATATACCTGGTAATATTATTCACTCCGAAGAAATACTGGTAATTTAATATAAATAATGCCCGTCAAGGa                                                                                                                                                      >  2:1666998/2‑150 (MQ=255)
   attCAGCGCCTATACCGATGCTAACTACCTCGGTTGTTGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTTGTAATATTCTTCACTCCGAATAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGAttt                                                                                                                                                   <  1:221487/149‑1 (MQ=255)
     tCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGa                                                                                                                                                 >  2:859577/1‑150 (MQ=255)
       aGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACAGGTCCGCGAACGAGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTCGTAATTTAATCTAAATAATGCCCGTCAAGCATTTGACg                                                                                                                                               >  2:1276169/1‑150 (MQ=255)
          gcCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTCGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAATCCTTGACGGGc                                                                                                                                            <  1:2052569/150‑1 (MQ=255)
          gcCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAATCCTTGACGGGc                                                                                                                                            <  2:21385/150‑1 (MQ=255)
            cTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCAt                                                                                                                                          <  1:1773503/150‑1 (MQ=255)
            cTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAATCCTTGACGGGCAt                                                                                                                                          <  1:274029/150‑1 (MQ=255)
                  cGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGc                                                                                                                                    <  2:623208/150‑1 (MQ=255)
                  cGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCAAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGc                                                                                                                                    >  2:1112257/1‑150 (MQ=255)
                    aTGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCaa                                                                                                                                  <  2:1144451/150‑1 (MQ=255)
                      gCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAg                                                                                                                                >  1:1591597/1‑150 (MQ=255)
                      gCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAg                                                                                                                                <  1:79901/150‑1 (MQ=255)
                       cTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGa                                                                                                                                                      <  1:1162950/127‑1 (MQ=255)
                       cTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGa                                                                                                                                                      >  2:1162950/1‑127 (MQ=255)
                       cTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAgg                                                                                                                               <  1:515955/150‑1 (MQ=255)
                       cTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAgg                                                                                                                               <  1:2042706/150‑1 (MQ=255)
                        tAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCGGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAgaa                                                                                                                              >  2:762102/1‑148 (MQ=255)
                        tAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGa                                                                                                                              >  1:1671364/1‑150 (MQ=255)
                            tACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGAcgcg                                                                                                                          <  1:101419/150‑1 (MQ=255)
                                tCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACggg                                                                                                                                             >  2:1348234/1‑127 (MQ=255)
                                tCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACggg                                                                                                                                             <  1:1348234/127‑1 (MQ=255)
                                  ggtggtggtGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATg                                                                                                                    <  2:340319/150‑1 (MQ=255)
                                     ggtggtgACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTa                                                                                                                 <  1:1879624/150‑1 (MQ=255)
                                       tggtgACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATc                                                                                                               <  2:397891/150‑1 (MQ=255)
                                       tggtgACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATc                                                                                                               <  2:283150/150‑1 (MQ=255)
                                       tggtgACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATc                                                                                                               <  1:565154/150‑1 (MQ=255)
                                        ggtcACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTATTCACTCAGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGAATTTGACGGGCATTACTGAAAAGGACGCCCAAATGTTATCt                                                                                                              >  2:2479836/1‑150 (MQ=255)
                                           gACGTAGATCAAGACTGGTCCGCGAACGGGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAATTAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTa                                                                                                           <  2:1800734/150‑1 (MQ=255)
                                             cGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATAGCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGt                                                                                                         <  2:197203/150‑1 (MQ=255)
                                             cGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGt                                                                                                         >  1:91998/1‑150 (MQ=255)
                                             cGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGt                                                                                                         >  2:2374030/1‑150 (MQ=255)
                                             cGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTATTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGt                                                                                                         <  2:188560/150‑1 (MQ=255)
                                                 gATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTaaa                                                                                                     >  1:1865782/1‑150 (MQ=255)
                                                  aTCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGt                                                                                                         <  1:2307214/145‑1 (MQ=255)
                                                  aTCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGt                                                                                                         >  2:2307214/1‑145 (MQ=255)
                                                   tCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAAcc                                                                                                   >  1:1588483/1‑150 (MQ=255)
                                                      aGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACctct                                                                                                >  1:1701313/1‑150 (MQ=255)
                                                        aCTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACTGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTc                                                                                              <  1:2464736/150‑1 (MQ=255)
                                                        aCTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTc                                                                                              <  2:861754/150‑1 (MQ=255)
                                                         cTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCa                                                                                             <  1:603166/150‑1 (MQ=255)
                                                          tGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAg                                                                                            <  2:846142/150‑1 (MQ=255)
                                                          tGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAg                                                                                            <  1:1156225/150‑1 (MQ=255)
                                                           ggTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAgg                                                                                           <  1:2004409/150‑1 (MQ=255)
                                                              ccGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAAt                                                                                        <  2:2309627/150‑1 (MQ=255)
                                                                gcgAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAAttt                                                                                      >  2:2200819/1‑150 (MQ=255)
                                                                gcgAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAAttt                                                                                      <  1:23240/150‑1 (MQ=255)
                                                                  gAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAAttttg                                                                                    <  2:1973618/150‑3 (MQ=255)
                                                                   aaCGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTTTAGTTAAACCTCTTCAGGATTTTGGt                                                                                   <  2:1840057/150‑1 (MQ=255)
                                                                         ggTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTc                                                                             >  2:2236590/1‑150 (MQ=255)
                                                                         ggTGTTAAATATACCTGGTAAAATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTc                                                                             <  1:380781/150‑1 (MQ=255)
                                                                            gtTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGAt                                                                          <  2:1013978/150‑1 (MQ=255)
                                                                            gtTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGAt                                                                          >  1:464564/1‑150 (MQ=255)
                                                                             ttAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATa                                                                         <  2:2085207/150‑1 (MQ=255)
                                                                              tAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATaa                                                                        >  1:75878/1‑150 (MQ=255)
                                                                                aaTATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAAt                                                                      <  1:1648116/150‑1 (MQ=255)
                                                                                aaTATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAAt                                                                      >  1:886933/1‑150 (MQ=255)
                                                                                 atatACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATg                                                                     <  2:1518385/150‑1 (MQ=255)
                                                                                    taCCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAGCCTCTTCAGGAATTTGGTAAGCTCGATAAATGttt                                                                  >  2:2146822/1‑150 (MQ=255)
                                                                                    taCCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGttt                                                                  >  2:2037137/1‑150 (MQ=255)
                                                                                     aCCTGGTAATATTCTTCACTCCGAAGAAATACTGGTATTTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAATCGTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTg                                                                 <  1:832200/150‑1 (MQ=255)
                                                                                      ccTGGTAATATTCTTCACTCCGAAGAAATAGTGGTAATTTAATCTAAATAATGCCCGTCAATTATTTGACGGTCATTACTTCAAATGACGCGCATATGTTATCTTTAGTTAAACCTCTTCAGTAATTTGGTAAGCTCGATAAATGTTTGt                                                                <  2:1314512/150‑1 (MQ=255)
                                                                                        tGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTcc                                                              >  2:1169265/1‑150 (MQ=255)
                                                                                         ggTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCa                                                             >  1:1438036/1‑150 (MQ=255)
                                                                                         ggTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCa                                                             <  2:1458026/150‑1 (MQ=255)
                                                                                         ggTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCa                                                             <  2:1516983/150‑1 (MQ=255)
                                                                                         ggTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCa                                                             >  2:1771015/1‑150 (MQ=255)
                                                                                          gTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAg                                                            <  2:546145/150‑1 (MQ=255)
                                                                                                  cttcACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTa                                                    >  2:2282310/1‑150 (MQ=255)
                                                                                                       cTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACgcgc                                               >  1:268101/1‑150 (MQ=255)
                                                                                                         ccGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCtt                                             >  1:1088825/1‑150 (MQ=255)
                                                                                                          cGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTc                                            >  2:758530/1‑150 (MQ=255)
                                                                                                           gaagaaATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCg                                           <  1:6750/150‑1 (MQ=255)
                                                                                                                aaTACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGttt                                      >  1:1501414/1‑150 (MQ=255)
                                                                                                                 aTACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTta                                     <  1:211349/150‑1 (MQ=255)
                                                                                                                  tACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTtaa                                    <  1:1935584/150‑1 (MQ=255)
                                                                                                                     tGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCCCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTtaataa                                 >  2:495635/1‑150 (MQ=255)
                                                                                                                      ggTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTaataat                                >  1:1197878/1‑150 (MQ=255)
                                                                                                                       gTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGTGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCATGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATg                               <  2:976692/150‑1 (MQ=255)
                                                                                                                              aaTCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCa                        >  1:1838503/1‑150 (MQ=255)
                                                                                                                               aTCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCaa                       >  1:2090485/1‑150 (MQ=255)
                                                                                                                                 cTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGt                     <  1:859577/150‑1 (MQ=255)
                                                                                                                                  tAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGtt                    >  1:1848245/1‑150 (MQ=255)
                                                                                                                                     ataatGCCCGTCAATGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTtata                 <  1:1338259/150‑1 (MQ=255)
                                                                                                                                     ataatGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTtata                 <  1:1523645/150‑1 (MQ=255)
                                                                                                                                        atGCCCGTCAAATCCTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATAttt              >  1:1139686/1‑150 (MQ=255)
                                                                                                                                         tGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATAtttt             >  2:1539022/1‑150 (MQ=255)
                                                                                                                                         tGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATAtttt             >  2:1530091/1‑150 (MQ=255)
                                                                                                                                          gCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTc            >  1:935669/1‑150 (MQ=255)
                                                                                                                                           cccGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTcc           <  1:803905/150‑1 (MQ=255)
                                                                                                                                           cccGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTcc           >  2:831904/1‑150 (MQ=255)
                                                                                                                                           cccGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTcc           <  2:2199804/150‑1 (MQ=255)
                                                                                                                                             cGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAg         <  1:1276169/150‑1 (MQ=255)
                                                                                                                                                cAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGa      >  2:1941980/1‑150 (MQ=255)
                                                                                                                                                cAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGa      >  2:1427454/1‑150 (MQ=255)
                                                                                                                                                cAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGAGCTTCGAGTTTAATAATGAAAAGCAAGTTATTTTTTCCAGTGa      >  2:1615502/1‑150 (MQ=255)
                                                                                                                                                  aGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGATg    <  2:1701313/150‑1 (MQ=255)
                                                                                                                                                    taTTTGACGGGCATTACTGCAAAGTACGCGCAAATGTTTTATGTAGTTAAACCTCATCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACTCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGATGTc  <  2:733022/149‑1 (MQ=255)
                                                                                                                                                    gATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGATGTc  <  2:217996/150‑1 (MQ=255)
                                                                                                                                                    gATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGATGTc  <  2:2090485/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGATGTC  >  NC_000913/394325‑394622

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: