Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,054,977 T→G 6.4% intergenic (+18/‑84) 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,9770TG 6.4% 386.3 / 11.4 124intergenic (+18/‑84)shiA/amnshikimate transporter/AMP nucleosidase
Reads supporting (aligned to +/- strand):  ref base T (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 = 8.74e-01

Only 100 of 127 total aligned reads displayed.

CACCTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATG  >  NC_000913/2054829‑2055124
                                                                                                                                                    |                                                                                                                                                   
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)
        aTTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTc                                                                                                                                            >  2:1471998/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)
            ccGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCGGATGATCTTCGcca                                                                                                                                        <  1:1234859/150‑2 (MQ=255)
            ccGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCt                                                                                                                                        <  1:1674429/150‑1 (MQ=255)
             cGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCtt                                                                                                                                       >  1:1460845/1‑150 (MQ=255)
                 gCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTacac                                                                                                                                   >  1:705919/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)
                          acttacttTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCGTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTc                                                                                                                          <  2:1392364/150‑1 (MQ=255)
                            ttccttTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTcac                                                                                                                        <  2:1312306/147‑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)
                                       ggAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGAc                                                                                                             <  2:559748/150‑1 (MQ=255)
                                           cTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATAc                                                                                                         <  1:148302/150‑1 (MQ=255)
                                           cTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATAc                                                                                                         >  2:1712807/1‑150 (MQ=255)
                                           cTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATAc                                                                                                         >  1:1100734/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)
                                                        gCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCg                                                                                            >  1:615821/1‑150 (MQ=255)
                                                         ccATTTATTTGCTTGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCgg                                                                                           <  2:1545636/150‑1 (MQ=255)
                                                         ccATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCgg                                                                                           <  1:152879/150‑1 (MQ=255)
                                                         ccATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCgg                                                                                           <  2:28151/150‑1 (MQ=255)
                                                           atttatttGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTa                                                                                         <  2:1288395/150‑1 (MQ=255)
                                                              tatttGCTGGCTGGATGCCTGATTTCCGCAATGCCCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAAtt                                                                                      >  1:1328607/1‑150 (MQ=255)
                                                               atttGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACATTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTAACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTg                                                                                     <  2:1190249/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)
                                                                   gctggctgGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATATGCGGTAATTGTATg                                                                                 <  1:751031/150‑1 (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:1721031/1‑150 (MQ=255)
                                                                           gATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGga                                                                         <  2:52782/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)
                                                                                cTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGacac                                                                    >  1:409416/1‑150 (MQ=255)
                                                                                  gATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGacacac                                                                  >  2:179959/1‑150 (MQ=255)
                                                                                      tCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGa                                                              >  2:1536545/1‑150 (MQ=255)
                                                                                       ccGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGaa                                                             <  1:215237/150‑1 (MQ=255)
                                                                                        cGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGaat                                                            <  2:181256/150‑1 (MQ=255)
                                                                                        cGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGaat                                                            <  1:1788712/150‑1 (MQ=255)
                                                                                              gACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAg                                                      >  2:899523/1‑150 (MQ=255)
                                                                                               aCCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAgg                                                     <  1:1544617/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)
                                                                                                         tgatgaAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCt                                           >  1:871265/1‑150 (MQ=255)
                                                                                                          tatgaAAGACAGTCAACGCGCTTGATAGCCTGTCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAACGGCTCCGGTCTg                                          <  2:505559/149‑1 (MQ=255)
                                                                                                            tgaAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGAc                                        >  2:1157784/1‑150 (MQ=255)
                                                                                                            tgaAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGAc                                        >  1:222521/1‑150 (MQ=255)
                                                                                                             gaAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTTCACTTTTGTTTCACATTTATGTGACATACTATCGGATGTGCGGTTATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGAcc                                       >  2:518268/1‑150 (MQ=255)
                                                                                                             gaAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGAcc                                       <  2:1328607/150‑1 (MQ=255)
                                                                                                             gaAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGAcc                                       >  1:788065/1‑150 (MQ=255)
                                                                                                             gaAAGACAGTAAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTTCACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATATTAAGGGCTACGGTCTGAcc                                       >  1:1799019/1‑150 (MQ=255)
                                                                                                             caaaGACAGTCCACGCGCTTGCTACCCTGGCGAAGATCATCCGATCCTCGCCGTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATCAATAATAAGGGCTCCGGGCTGAcc                                       <  2:1005275/149‑1 (MQ=255)
                                                                                                               aaGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGAcccc                                     <  1:1626882/150‑1 (MQ=255)
                                                                                                                 gACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAg                                   >  2:1186247/1‑150 (MQ=255)
                                                                                                                 gACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAg                                   <  1:1131255/150‑1 (MQ=255)
                                                                                                                  aCAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACAATCGGATGTGCGGTAAATGTAGGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGc                                  <  2:332009/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)
                                                                                                                       cAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAgg                             >  2:1080285/1‑150 (MQ=255)
                                                                                                                       cAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGGAATTGTATGGAACAGGAGACACACATGAATAAGAAGGGCTCCGGTCTGACCCCAGCTCAgg                             >  2:35237/1‑150 (MQ=255)
                                                                                                                            cgcTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTg                        <  1:1569996/150‑1 (MQ=255)
                                                                                                                                tGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATa                    >  2:1292051/1‑150 (MQ=255)
                                                                                                                                tGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATa                    >  1:755257/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)
                                                                                                                                              aGATCATCCTATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTg      <  2:310201/150‑1 (MQ=255)
                                                                                                                                               gATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGt     >  1:10140/1‑150 (MQ=255)
                                                                                                                                                  catcCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGATTAATAAGGGCTCCGGTCTGACCCCCGCTCAGTCACTGGATAAACTCGACGCGCTGTATg  >  2:367402/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CACCTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATG  >  NC_000913/2054829‑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: