Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,054,978 C→A 6.2% intergenic (+19/‑83) 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,9780CA 6.2% 332.9 / 10.9 122intergenic (+19/‑83)shiA/amnshikimate transporter/AMP nucleosidase
Reads supporting (aligned to +/- strand):  ref base C (57/57);  new base A (2/6);  total (59/63)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

Only 100 of 128 total aligned reads displayed.

CACCTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGC  >  NC_000913/2054829‑2055127
                                                                                                                                                     |                                                                                                                                                     
cacCTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGatcatc                                                                                                                                                       >  2:704315/1‑150 (MQ=255)
 acCTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATcc                                                                                                                                                      <  2:662260/150‑1 (MQ=255)
 acCTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATcc                                                                                                                                                      >  2:277221/1‑150 (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                                                                                                                                            <  2:437741/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)
              gCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTa                                                                                                                                         <  2:209922/150‑1 (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:52782/1‑150 (MQ=255)
                                           cTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATAc                                                                                                            >  2:1712807/1‑150 (MQ=255)
                                           cTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATAc                                                                                                            >  1:1100734/1‑150 (MQ=255)
                                             ggCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTa                                                                                                          >  1:975244/1‑150 (MQ=255)
                                               cATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATc                                                                                                        >  2:649193/1‑150 (MQ=255)
                                               cATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATc                                                                                                        >  1:1870730/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)
                                                          cATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGt                                                                                             <  1:649193/150‑1 (MQ=255)
                                                           atttatttGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTa                                                                                            <  1:277221/150‑1 (MQ=255)
                                                           atttatttGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTa                                                                                            <  2:1288395/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)
                                                                             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)
                                                                                  gATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGacacac                                                                     >  2:179959/1‑150 (MQ=255)
                                                                                       ccGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGaa                                                                <  1:215237/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)
                                                                                                             gaAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTTCACTTTTGTTTCACATTTATGTGACATACTATCGGATGTGCGGTTATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGAcc                                          >  2:518268/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                                      <  1:1131255/150‑1 (MQ=255)
                                                                                                                 gACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAg                                      >  2:1186247/1‑150 (MQ=255)
                                                                                                                  aCAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACAATCGGATGTGCGGTAAATGTAGGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGc                                     <  2:332009/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)
                                                                                                                           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:1822456/1‑150 (MQ=255)
                                                                                                                                           cGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGAcgcg            <  1:1790577/150‑1 (MQ=255)
                                                                                                                                              aGATCATCCTATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTg         <  2:310201/150‑1 (MQ=255)
                                                                                                                                               gATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGt        >  1:10140/1‑150 (MQ=255)
                                                                                                                                                  catcCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGATTAATAAGGGCTCCGGTCTGACCCCCGCTCAGTCACTGGATAAACTCGACGCGCTGTATg     >  2:367402/1‑150 (MQ=255)
                                                                                                                                                     ccGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGc  <  2:566440/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACCTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGAGC  >  NC_000913/2054829‑2055127

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