Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,054,976 A→G 6.7% intergenic (+17/‑85) shiA → / → amn shikimate transporter/AMP nucleosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,054,9760AG 6.7% 341.2 / 12.3 119intergenic (+17/‑85)shiA/amnshikimate transporter/AMP nucleosidase
Reads supporting (aligned to +/- strand):  ref base A (61/50);  new base G (2/6);  total (63/56)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.12e-01

Only 100 of 125 total aligned reads displayed.

TACACCTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTA  >  NC_000913/2054827‑2055122
                                                                                                                                                     |                                                                                                                                                  
tACACCTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGatca                                                                                                                                                    >  2:1569347/1‑150 (MQ=255)
tACACCTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGatca                                                                                                                                                    >  2:1372163/1‑150 (MQ=255)
tACACCTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGatca                                                                                                                                                    >  1:911617/1‑150 (MQ=255)
   acCTTTTATTGCCGCTGCCCTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCAGTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGCCCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATcc                                                                                                                                                 >  1:226987/1‑150 (MQ=255)
   acCTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTAATTTGCTGGCTGGATGCCTGATTACCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATcc                                                                                                                                                 >  1:825968/1‑150 (MQ=255)
    ccTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCg                                                                                                                                                >  1:190961/1‑150 (MQ=255)
     cTTTTATTGCCGCTGCACTAATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTATGTTGATGAAAGAAAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGa                                                                                                                                               >  2:312188/1‑150 (MQ=255)
      ttttATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGAt                                                                                                                                              >  1:1526038/1‑150 (MQ=255)
       tttATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCAAAGATCATCCGATc                                                                                                                                             <  2:647377/150‑1 (MQ=255)
           ttGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCGGATGATCTTCg                                                                                                                                         <  2:1107016/150‑1 (MQ=255)
             gCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCGGATGATCTTCGcc                                                                                                                                       <  1:1050525/150‑1 (MQ=255)
             gCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCGGATGATCTTCGcc                                                                                                                                       <  2:1119883/150‑1 (MQ=255)
             gCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGcc                                                                                                                                       >  2:952735/1‑150 (MQ=255)
             gCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGcc                                                                                                                                       >  2:1003477/1‑150 (MQ=255)
             gCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGcc                                                                                                                                       >  1:501035/1‑150 (MQ=255)
              ccGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCtt                                                                                                                                     <  2:468716/150‑1 (MQ=255)
               cGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCGGATGATCTTCGccag                                                                                                                                     <  1:1424059/150‑3 (MQ=255)
                gCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCGGATGATCTTCGccagg                                                                                                                                    <  2:12057/150‑4 (MQ=255)
                gCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGCTCTTCGCCTTa                                                                                                                                    >  1:333239/1‑150 (MQ=255)
                 cTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGCCAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTac                                                                                                                                   >  1:158102/1‑150 (MQ=255)
                 cTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCttcc                                                                                                                                   >  1:125986/1‑148 (MQ=255)
                   gCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTacac                                                                                                                                 >  1:1314228/1‑150 (MQ=255)
                    cACTCATCACTTACTTTGCCGGGATCTGGCATAGCGTCGCCATTTATTTGCTGGCTCGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACt                                                                                                                                >  1:104095/1‑150 (MQ=255)
                    cACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACt                                                                                                                                >  1:99481/1‑150 (MQ=255)
                        catcaCTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTg                                                                                                                            >  2:540166/1‑150 (MQ=255)
                            acttacttTGCCTGGAACTGGCATAGCGTCGCCATTTATTTTCTGTCTGGATGCCTGATTTACGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTc                                                                                                                        <  1:29339/150‑1 (MQ=255)
                            acttacttTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTc                                                                                                                        <  2:812285/150‑1 (MQ=255)
                              ttacttTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTcac                                                                                                                      <  2:914261/150‑1 (MQ=255)
                               tacttTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTcaca                                                                                                                     >  1:1406604/1‑150 (MQ=255)
                                    tGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCt                                                                                                                >  2:1333295/1‑150 (MQ=255)
                                     gCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCtg                                                                                                               >  2:877201/1‑150 (MQ=255)
                                      ccGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCtgt                                                                                                              >  2:1682410/1‑150 (MQ=255)
                                      ccGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCtgt                                                                                                              >  2:40127/1‑150 (MQ=255)
                                       cGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCtgtg                                                                                                             >  1:493669/1‑150 (MQ=255)
                                       cGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCtgtg                                                                                                             >  2:156529/1‑150 (MQ=255)
                                       cGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCtgtg                                                                                                             <  2:916752/150‑1 (MQ=255)
                                                   tAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCgg                                                                                                 <  2:501035/150‑1 (MQ=255)
                                                    aGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGa                                                                                                <  1:1232165/150‑1 (MQ=255)
                                                       gTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGAtgt                                                                                             >  1:621827/1‑150 (MQ=255)
                                                          gCCATTTATTTGCTGGCTGGATGCCTGGTTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACCTTTCGGTGACATACTATCGGATGTGCg                                                                                          >  1:952003/1‑150 (MQ=255)
                                                          gCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCg                                                                                          <  1:540166/150‑1 (MQ=255)
                                                           ccATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCgg                                                                                         >  1:131266/1‑150 (MQ=255)
                                                            cATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGt                                                                                        >  2:231114/1‑150 (MQ=255)
                                                            cATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGt                                                                                        >  1:917676/1‑150 (MQ=255)
                                                             atttatttGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTa                                                                                       <  2:601858/150‑1 (MQ=255)
                                                                      ctggctggATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATgg                                                                              <  1:231114/150‑1 (MQ=255)
                                                                           tggATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAg                                                                         <  1:1682410/150‑1 (MQ=255)
                                                                            ggATGCCTGATTTCCGCAATGACCTCTTTGTTGATGAAAGACAGTCAACGCTCTTGATAGTCTGGCGAAGATCATCCGATCTTCGACTTACACTTTTGTTTCACATTTCTGTGACATTCTATCGGATGTGCGGTAATTGTATGGAACAgg                                                                        <  1:97448/150‑1 (MQ=255)
                                                                             gATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGga                                                                       <  1:1569347/150‑1 (MQ=255)
                                                                                gCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGac                                                                    >  2:1320896/1‑150 (MQ=255)
                                                                                gCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGac                                                                    <  1:167752/150‑1 (MQ=255)
                                                                                  cTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGacac                                                                  >  2:1343300/1‑150 (MQ=255)
                                                                                   tGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGacaca                                                                 <  1:844854/150‑1 (MQ=255)
                                                                                    gATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGacacac                                                                <  1:376127/150‑1 (MQ=255)
                                                                                     aTTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGAcacaca                                                               <  1:1372163/150‑1 (MQ=255)
                                                                                     aTTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGAcacaca                                                               <  1:1442096/150‑1 (MQ=255)
                                                                                      tttCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACAt                                                              >  2:1008594/1‑150 (MQ=255)
                                                                                       ttCCGCAATGACCTCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATg                                                             <  2:1360455/150‑1 (MQ=255)
                                                                                       ttCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATg                                                             <  1:952735/150‑1 (MQ=255)
                                                                                       ttCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATg                                                             <  1:172959/150‑1 (MQ=255)
                                                                                         ttggaaTGACATCTTTGTTGATGTAAGTAATTTAACTCGGTTGATTGCCTGGCGAATATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGaa                                                           <  1:1003477/146‑1 (MQ=255)
                                                                                          cGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGaat                                                          >  2:1279328/1‑150 (MQ=255)
                                                                                            cAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGaataa                                                        <  2:1406604/150‑1 (MQ=255)
                                                                                               tGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAAtaataa                                                     >  1:1694145/1‑150 (MQ=255)
                                                                                                  ccGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGTGACACACATGAATAATAAggg                                                  >  1:900219/1‑150 (MQ=255)
                                                                                                      tttGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTcc                                              >  1:501359/1‑150 (MQ=255)
                                                                                                         gttgATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACCCTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGt                                           >  2:1203394/1‑150 (MQ=255)
                                                                                                         gttgATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCAGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGt                                           >  1:377983/1‑150 (MQ=255)
                                                                                                           tgatgaAAGACAGTGAACGCGGTTGATAGCCTGGCGAAGATCATCCGATCTTTGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCTGTAATTGTATGGAACAGTAGACACACATGAATAATATGGTCTCCCGTCt                                         <  2:88017/150‑1 (MQ=255)
                                                                                                           tgatgaAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCt                                         >  1:1073205/1‑150 (MQ=255)
                                                                                                              tgaAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGAc                                      >  1:779136/1‑150 (MQ=255)
                                                                                                              tgaAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGAc                                      <  1:1320896/150‑1 (MQ=255)
                                                                                                                aaaGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATTAATAATACGGGCTCCGGTCTGAccc                                    <  2:952003/150‑1 (MQ=255)
                                                                                                                aaaGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGAccc                                    <  1:177287/150‑1 (MQ=255)
                                                                                                                    aCAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGc                                >  1:1610676/1‑150 (MQ=255)
                                                                                                                    aCAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGc                                <  1:40127/150‑1 (MQ=255)
                                                                                                                     cAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCt                               >  1:417254/1‑150 (MQ=255)
                                                                                                                      aGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTc                              <  2:1526038/150‑1 (MQ=255)
                                                                                                                               gcTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTgg                     >  2:461412/1‑150 (MQ=255)
                                                                                                                                cTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGa                    >  1:188630/1‑150 (MQ=255)
                                                                                                                                cTTGATAGCCTGGCGAAGATCATCCGATATTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTAACGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCAAGCTCATGCACTGGa                    >  1:850301/1‑150 (MQ=255)
                                                                                                                                 ttGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGAt                   >  1:564912/1‑150 (MQ=255)
                                                                                                                                  tGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATg                                                                               <  1:515901/89‑1 (MQ=255)
                                                                                                                                  tGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGTGCTCCGGTCTGACCCCAGCTCAGGCACTGGATa                  >  1:885160/1‑150 (MQ=255)
                                                                                                                                  tGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATa                  >  2:231167/1‑150 (MQ=255)
                                                                                                                                  tGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAAGTGTATGGAACAGGAGACACACATGAAGAATAAGGGCTCCGGTATGACCCCCGCTCAGGCACTGGATa                  >  1:1131047/1‑150 (MQ=255)
                                                                                                                                  tGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTAAACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATg                                                                               >  2:515901/1‑89 (MQ=255)
                                                                                                                                      aGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACt              <  2:621827/150‑1 (MQ=255)
                                                                                                                                       gCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAtctc             >  2:791802/1‑150 (MQ=255)
                                                                                                                                       gCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTc             <  1:1279328/150‑1 (MQ=255)
                                                                                                                                         cTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGa           >  1:951538/1‑150 (MQ=255)
                                                                                                                                           ggCGAAGATCGGATGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGAcg         >  1:451678/1‑150 (MQ=255)
                                                                                                                                           ggCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGAcg         <  1:1412469/150‑1 (MQ=255)
                                                                                                                                           ggCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGAcg         <  2:1314228/150‑1 (MQ=255)
                                                                                                                                             cGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGAcgcg       >  2:1475304/1‑150 (MQ=255)
                                                                                                                                              gAAGATCATCAGATCTTCGCCTTACACTTTTGTTTCACTTTTCTGTGACATACTATCGGATGTCCGTTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCTGTATGACCCCACCTCAGTCACAGGAAAAACTCCAcgctc      >  2:38332/1‑148 (MQ=255)
                                                                                                                                               aaGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTAGCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCt     <  1:947976/150‑1 (MQ=255)
                                                                                                                                                aGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCTGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTg    <  2:190961/150‑1 (MQ=255)
                                                                                                                                                  atcatcCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTa  <  1:156529/150‑1 (MQ=255)
                                                                                                                                                  atcatcCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTa  <  2:333239/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TACACCTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTA  >  NC_000913/2054827‑2055122

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