Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,054,976 A→G 6.3% 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.3% 375.0 / 8.7 124intergenic (+17/‑85)shiA/amnshikimate transporter/AMP nucleosidase
Reads supporting (aligned to +/- strand):  ref base A (58/58);  new base G (2/6);  total (60/64)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 128 total aligned reads displayed.

TACACCTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATG  >  NC_000913/2054827‑2055124
                                                                                                                                                     |                                                                                                                                                    
tACACCTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGatca                                                                                                                                                      >  2:1569996/1‑150 (MQ=255)
  cacCTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGatcatc                                                                                                                                                    >  2:704315/1‑150 (MQ=255)
   acCTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATcc                                                                                                                                                   <  2:662260/150‑1 (MQ=255)
    ccTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGCCAGTCAACGCGCTTGCTAGCCTGGCGAAGCTCATCCg                                                                                                                                                  >  2:537145/1‑150 (MQ=255)
    ccTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCg                                                                                                                                                  >  1:1633316/1‑150 (MQ=255)
      ttttATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTTGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGATTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGAt                                                                                                                                                >  1:1394465/1‑150 (MQ=255)
          aTTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTc                                                                                                                                            >  2:1471998/1‑150 (MQ=255)
            tGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGc                                                                                                                                          >  2:1131255/1‑150 (MQ=255)
             gCCGCTTCACTCATCACTTACTTTGCCGTGAACTTGCATAGCGTCGCAATTTATTTGCTGGCTGGATGCGTGATTTCCGCAATGACCGCTTTGTTGTTGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCGGATGATCTTCGcc                                                                                                                                         <  2:215384/150‑1 (MQ=255)
             gCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCGGATGATCTTCGcc                                                                                                                                         <  1:1638414/150‑1 (MQ=255)
             gCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCGGATGATCTTCGcc                                                                                                                                         <  1:1435959/150‑1 (MQ=255)
             gCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCACTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCGGATGATCTTCGcc                                                                                                                                         <  2:1529115/150‑1 (MQ=255)
              ccGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCGGATGATCTTCGcca                                                                                                                                        <  1:1234859/150‑2 (MQ=255)
               cGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCtt                                                                                                                                       >  1:1460845/1‑150 (MQ=255)
                   gCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTacac                                                                                                                                   >  1:705919/1‑150 (MQ=255)
                     aCTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACtt                                                                                                                                 >  2:1788712/1‑150 (MQ=255)
                      cTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACttt                                                                                                                                >  2:215237/1‑150 (MQ=255)
                           caCTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGttt                                                                                                                           <  1:1281169/150‑1 (MQ=255)
                            acttacttTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTc                                                                                                                          <  1:1732285/150‑1 (MQ=255)
                                 cttTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGAGTGCCTGGCGAAGATCATCCGATCTTCGCCTTACACGTTTGTTTCAGAtt                                                                                                                     >  2:312448/1‑150 (MQ=255)
                                   ttGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTc                                                                                                                   <  1:854492/150‑1 (MQ=255)
                                      ccGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTAGCTGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCtgt                                                                                                                >  1:1499348/1‑150 (MQ=255)
                                       cGGGAACTGGGATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCGTTGTTGAGGAAAGACAGTCCACGCGCTTGATAGCCTGGCGAAGAGGATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCtgtg                                                                                                               <  2:1123311/150‑1 (MQ=255)
                                        gggAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGa                                                                                                              <  2:515854/150‑1 (MQ=255)
                                         ggAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGAc                                                                                                             <  2:559748/150‑1 (MQ=255)
                                             cTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATAc                                                                                                         <  1:148302/150‑1 (MQ=255)
                                             cTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATAc                                                                                                         >  1:1100734/1‑150 (MQ=255)
                                             cTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATAc                                                                                                         >  1:52782/1‑150 (MQ=255)
                                               ggCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTa                                                                                                       >  1:975244/1‑150 (MQ=255)
                                                 cATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATc                                                                                                     >  1:1870730/1‑150 (MQ=255)
                                                 cATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATc                                                                                                     >  2:649193/1‑150 (MQ=255)
                                                     gCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGAt                                                                                                 <  1:1611264/150‑1 (MQ=255)
                                                      cGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGAtg                                                                                                >  1:1006535/1‑150 (MQ=255)
                                                         cGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGc                                                                                             <  2:1867548/150‑1 (MQ=255)
                                                          gCCATTTATTTGCTTGCTGGATGCCTGATTTCAGCAATGACCGCTTTTTTGATGAAAGACAGTCAACTCGCTTGATATCCTTTCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCg                                                                                            <  1:1262043/150‑1 (MQ=255)
                                                           ccATTTATTTGCTTGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCgg                                                                                           <  2:1545636/150‑1 (MQ=255)
                                                           ccATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCgg                                                                                           <  1:152879/150‑1 (MQ=255)
                                                           ccATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCgg                                                                                           <  2:28151/150‑1 (MQ=255)
                                                            cATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGt                                                                                          <  1:649193/150‑1 (MQ=255)
                                                             atttatttGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTa                                                                                         <  2:1288395/150‑1 (MQ=255)
                                                             atttatttGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTa                                                                                         <  1:277221/150‑1 (MQ=255)
                                                                 atttGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTg                                                                                     <  2:975244/150‑1 (MQ=255)
                                                                   ttGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTa                                                                                   <  2:1534485/150‑1 (MQ=255)
                                                                   ttGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCGGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTa                                                                                   <  1:1031342/150‑1 (MQ=255)
                                                                    tGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTAt                                                                                  >  1:950079/1‑150 (MQ=255)
                                                                     gctggctgGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGGCAACGCGCTTGATAGGCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGTTAATTGTATg                                                                                 <  1:1626050/150‑1 (MQ=255)
                                                                     gctgactgGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATg                                                                                 >  1:1154801/1‑150 (MQ=255)
                                                                      cttgcttgATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGAGCGTGATAGCCTGGCGAATTTCATCCGATCTTCGACTTAAACTTTTGTTTCACATTTCTGTGAAATACTATCGGATGTGCGGTAATTGTATgg                                                                                <  1:497054/147‑1 (MQ=255)
                                                                           tggATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAg                                                                           <  2:210169/150‑1 (MQ=255)
                                                                            ggATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAgg                                                                          <  2:500659/150‑1 (MQ=255)
                                                                             gATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGga                                                                         <  2:52782/150‑1 (MQ=255)
                                                                             gATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGga                                                                         >  2:1721031/1‑150 (MQ=255)
                                                                               tGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGgaga                                                                       <  1:787257/150‑1 (MQ=255)
                                                                                gCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGac                                                                      >  1:164334/1‑150 (MQ=255)
                                                                                 ccTGATTTCCGCAATGACCGCTTAGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCtgtg                                                                                                               >  1:721987/1‑108 (MQ=255)
                                                                                 ccTGATTTCCGCAATGACCGCTTAGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCtgtg                                                                                                               <  2:721987/108‑1 (MQ=255)
                                                                                  cggATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGacac                                                                    <  2:503209/148‑1 (MQ=255)
                                                                                         ccGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGaa                                                             <  1:215237/150‑1 (MQ=255)
                                                                                          cGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGaat                                                            <  2:181256/150‑1 (MQ=255)
                                                                                          cGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGaat                                                            <  1:1788712/150‑1 (MQ=255)
                                                                                              aTGACCGCTTTGTGGATGAAAGACAGTCAACGCGCTCGATAGCCGGGCGAAGAGCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAataata                                                        <  2:239052/150‑1 (MQ=255)
                                                                                                gACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAg                                                      >  2:899523/1‑150 (MQ=255)
                                                                                                 aCCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAgg                                                     <  2:1335413/150‑1 (MQ=255)
                                                                                                 aCCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAgg                                                     <  2:1534049/150‑1 (MQ=255)
                                                                                                  ccGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAggg                                                    >  2:174820/1‑150 (MQ=255)
                                                                                                      tttGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTcc                                                <  1:408616/150‑1 (MQ=255)
                                                                                                        tgttgATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTGACATTTCTGTGACATAGTATCTGATGTGCGGTAATTTTATGGAACAGGAGAAACACATGAATAATAAGGGCTCCgg                                              >  2:857636/1‑150 (MQ=255)
                                                                                                        tgttgATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCgg                                              >  1:1530702/1‑150 (MQ=255)
                                                                                                          ttgATGAAAGACAGTCAACGCGCGTGATAGCCTGGCGAAGATCATCCGATCTTCGACTTACAATTTTTTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATTATAAGGGCTACGGTc                                            <  2:1284504/150‑1 (MQ=255)
                                                                                                           tgatgaAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCt                                           <  1:1700116/150‑1 (MQ=255)
                                                                                                            tatgaAAGACAGTCAACGCGCTTGATAGCCTGTCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAACGGCTCCGGTCTg                                          <  2:505559/149‑1 (MQ=255)
                                                                                                              tgaAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGAc                                        >  2:1157784/1‑150 (MQ=255)
                                                                                                               gaAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGAcc                                       >  1:788065/1‑150 (MQ=255)
                                                                                                               gaAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGAcc                                       <  2:1328607/150‑1 (MQ=255)
                                                                                                               caaaGACAGTCCACGCGCTTGCTACCCTGGCGAAGATCATCCGATCCTCGCCGTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATCAATAATAAGGGCTCCGGGCTGAcc                                       <  2:1005275/149‑1 (MQ=255)
                                                                                                                   gACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAg                                   >  2:1186247/1‑150 (MQ=255)
                                                                                                                   gACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAg                                   <  1:1131255/150‑1 (MQ=255)
                                                                                                                    aCAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCGTACACTTTTGTTTGACATTTCTGTGACATACTATCGGATGTGCGGTATTTGTATGGAACAGGAGACACACATGAATATTAAGGGCTCCGGTCTGACCCCAGc                                  <  2:1633316/150‑1 (MQ=255)
                                                                                                                       gTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGCCACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCa                               >  1:568706/1‑150 (MQ=255)
                                                                                                                       gTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCCGCTcc                               >  2:1831368/1‑149 (MQ=255)
                                                                                                                       gTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCa                               >  2:1728638/1‑150 (MQ=255)
                                                                                                                        tCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAg                              >  1:1062532/1‑150 (MQ=255)
                                                                                                                         cAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGGTGTGCGGTAATTGTATGGTACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAgg                             >  1:56784/1‑150 (MQ=255)
                                                                                                                         cAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGGAATTGTATGGAACAGGAGACACACATGAATAAGAAGGGCTCCGGTCTGACCCCAGCTCAgg                             >  2:35237/1‑150 (MQ=255)
                                                                                                                             gcgcTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAATGGCTCCGGTCTGACCACAGCTCAGGCAcg                         >  1:1555361/1‑149 (MQ=255)
                                                                                                                              cgcTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTg                        <  1:1569996/150‑1 (MQ=255)
                                                                                                                                ctggATAGCCTGGCGAAGATCATCCTATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGa                      <  2:313649/147‑1 (MQ=255)
                                                                                                                                  tGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATa                    >  1:755257/1‑150 (MQ=255)
                                                                                                                                  tGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATa                    >  2:1292051/1‑150 (MQ=255)
                                                                                                                                  tGATAGCCTGGAGAACATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATa                    <  2:1345217/150‑1 (MQ=255)
                                                                                                                                    agaggcTGGCGAAGATAATCCGATCTTCGCCTTACAATTTTGTTTCACATTTCTGTGACATACTATCTGATGTGCGGTAATTGTATGCATCAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATaaa                  <  2:127018/145‑1 (MQ=255)
                                                                                                                                      aGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACt                >  1:1830898/1‑150 (MQ=255)
                                                                                                                                      aGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACt                >  1:1053701/1‑150 (MQ=255)
                                                                                                                                          aggCGAAGATCGGATGATCTTCGACTTACACTATTGTTACACATTTCTGTGACATACTATAGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTCACCCCAGCTCAGGCACTGGATAAACTCGAc            >  2:393208/2‑150 (MQ=255)
                                                                                                                                           ggCGAAGATCGGATGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGTACCAGGAGACACACATGAATAATAAGGGATCCGGTCTGACCCCAGATCATGCACTGGAGAAACTCGAcg           >  1:296845/1‑150 (MQ=255)
                                                                                                                                             cGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGAcgcg         <  1:1790577/150‑1 (MQ=255)
                                                                                                                                             cGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGAcgcg         >  1:1822456/1‑150 (MQ=255)
                                                                                                                                                aGATCATCCTATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTg      <  2:310201/150‑1 (MQ=255)
                                                                                                                                                 gATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGt     >  1:10140/1‑150 (MQ=255)
                                                                                                                                                    catcCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGATTAATAAGGGCTCCGGTCTGACCCCCGCTCAGTCACTGGATAAACTCGACGCGCTGTATg  >  2:367402/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TACACCTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATG  >  NC_000913/2054827‑2055124

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