Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 847,214 T→A 5.2% intergenic (‑95/+44) glnP ← / ← glnH glutamine transporter subunit/glutamine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913847,2140TA 5.2% 553.5 / 13.6 182intergenic (‑95/+44)glnP/glnHglutamine transporter subunit/glutamine transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (89/76);  new base A (6/4);  total (98/84)
Fisher's exact test for biased strand distribution p-value = 7.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 187 total aligned reads displayed.

CACCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTA  >  NC_000913/847065‑847363
                                                                                                                                                     |                                                                                                                                                     
  ccTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCGGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAg                                                                                                                                                     <  2:3210239/150‑1 (MQ=255)
    ttCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAggg                                                                                                                                                   >  1:3460139/1‑150 (MQ=255)
    ttCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAggg                                                                                                                                                   >  2:3289721/1‑150 (MQ=255)
    ttCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAggg                                                                                                                                                   >  1:1240259/1‑150 (MQ=255)
    ttCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAggg                                                                                                                                                   >  1:422054/1‑150 (MQ=255)
    ttCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAggg                                                                                                                                                   <  1:1520503/150‑1 (MQ=255)
        aTCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTATGGGgaga                                                                                                                                               <  2:3258383/150‑1 (MQ=255)
         tCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTATGGGGagag                                                                                                                                              <  1:3349012/150‑1 (MQ=255)
         tCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGagag                                                                                                                                              <  2:2674707/150‑1 (MQ=255)
           aGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTATGGGGAGAgga                                                                                                                                            <  1:137742/150‑2 (MQ=255)
           aGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAggg                                                                                                                                            <  2:589653/150‑1 (MQ=255)
             aaGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTATATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGaa                                                                                                                                          <  2:2022008/150‑1 (MQ=255)
             aaGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGaa                                                                                                                                          <  1:3472502/150‑1 (MQ=255)
               gCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTCATCCTCTCCCCATAGGGGAGAGGGACCt                                                                                                                                        >  1:697656/1‑150 (MQ=255)
               gCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACt                                                                                                                                        >  2:4060147/1‑150 (MQ=255)
               gCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACt                                                                                                                                        >  1:2022594/1‑150 (MQ=255)
                 ggAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAACCCTCTCCCCTATGGGGAGAggattaggg                                                                                                                                      <  1:1087530/150‑8 (MQ=255)
                  gAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGcc                                                                                                                                     >  2:3420525/1‑150 (MQ=255)
                    aTGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAg                                                                                                                                   >  1:300771/1‑150 (MQ=255)
                         aGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGt                                                                                                                              >  1:107579/1‑150 (MQ=255)
                          ggCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGtt                                                                                                                             <  1:3109832/150‑1 (MQ=255)
                             cAGATGGCACTCCAGTCAAACTGCATATTTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTa                                                                                                                          <  1:4047039/150‑1 (MQ=255)
                              aGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTAc                                                                                                                         <  2:2365932/150‑1 (MQ=255)
                              aGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTAc                                                                                                                         <  1:1163088/150‑1 (MQ=255)
                                aTGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAg                                                                                                                       <  2:703419/150‑1 (MQ=255)
                                 tGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAgg                                                                                                                      <  2:1887472/150‑1 (MQ=255)
                                 tGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAgg                                                                                                                      >  2:3389265/1‑150 (MQ=255)
                                 tGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAgg                                                                                                                      <  2:52828/150‑1 (MQ=255)
                                  ggCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTCGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGgt                                                                                                                     <  2:453231/150‑1 (MQ=255)
                                      cTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAg                                                                                                                 <  1:3012893/150‑1 (MQ=255)
                                       tCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGc                                                                                                                >  2:1830141/1‑150 (MQ=255)
                                       tCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGc                                                                                                                >  2:2802698/1‑150 (MQ=255)
                                       tCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCAGAGGGGAGTGGGACCTGCCAGTGCGGTTTACAGGTGTAGc                                                                                                                >  2:229482/1‑150 (MQ=255)
                                            tCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGttat                                                                                                           >  1:2717694/1‑150 (MQ=255)
                                            tCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGttat                                                                                                           >  2:1932639/1‑150 (MQ=255)
                                              aaaCTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGttatta                                                                                                         <  2:89955/150‑1 (MQ=255)
                                               aaCTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTtattat                                                                                                        >  1:358199/1‑150 (MQ=255)
                                                    cATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCgg                                                                                                   <  1:2245441/150‑1 (MQ=255)
                                                      tatGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGtt                                                                                                 >  1:3592194/1‑150 (MQ=255)
                                                             ttCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTAcc                                                                                          <  1:2945436/150‑1 (MQ=255)
                                                               ccTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGACCTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGa                                                                                        >  1:1724405/1‑150 (MQ=255)
                                                                cTGTTACCGTGGTGTAAATAATCAAATTCAGGGTTATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGaa                                                                                       <  2:1841230/150‑1 (MQ=255)
                                                                  gTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAAcc                                                                                     >  2:3590253/1‑150 (MQ=255)
                                                                    tACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCAt                                                                                   >  2:3670412/1‑150 (MQ=255)
                                                                       cGTGGTGTAAATAATCATATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCACATAGGGGAGTGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCAtttt                                                                                <  2:1394639/150‑1 (MQ=255)
                                                                         tggtgTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCAtttttt                                                                              >  1:3267846/1‑150 (MQ=255)
                                                                         tggtgTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCAtttttt                                                                              <  2:1344560/150‑1 (MQ=255)
                                                                         cggtgTAAATAATCAAAGTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCAtttttt                                                                              <  2:353528/149‑1 (MQ=255)
                                                                          ggtgTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCACTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTg                                                                             <  1:2124470/150‑1 (MQ=255)
                                                                            tgtAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTa                                                                           <  1:2475929/150‑1 (MQ=255)
                                                                             gtAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAg                                                                          >  2:2406514/1‑150 (MQ=255)
                                                                                aataatCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGAtt                                                                       >  1:483462/1‑150 (MQ=255)
                                                                                   aatCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCg                                                                    <  2:2088645/150‑1 (MQ=255)
                                                                                   aatCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCg                                                                    <  1:591094/150‑1 (MQ=255)
                                                                                    atCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTAGGTTCAGTACCGAACCATTTTTTGTAGATTTCGt                                                                   <  1:2697803/150‑1 (MQ=255)
                                                                                     tCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGtt                                                                  >  1:1417927/1‑150 (MQ=255)
                                                                                     tCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGtt                                                                  >  1:1417960/1‑150 (MQ=255)
                                                                                     tCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTGTTACAGGGGTAGCGTTATTATTTCGGTTCAGTACCGAACCATGTTTTGTAGATGTCGtt                                                                  >  1:2359723/1‑150 (MQ=255)
                                                                                     tCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCATCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTAAAGGGGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGtt                                                                  >  2:414524/1‑150 (MQ=255)
                                                                                      cAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTg                                                                 >  1:179808/1‑150 (MQ=255)
                                                                                      cAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTg                                                                 >  1:2988319/1‑150 (MQ=255)
                                                                                       aaaTTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGt                                                                >  1:87569/1‑150 (MQ=255)
                                                                                        aaTTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTa                                                               >  2:1300199/1‑150 (MQ=255)
                                                                                            cAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCAGAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGtt                                                           >  2:2662678/1‑150 (MQ=255)
                                                                                              gggTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTcc                                                         <  2:1621573/150‑1 (MQ=255)
                                                                                              gggTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTcc                                                         <  2:1029560/150‑1 (MQ=255)
                                                                                               ggTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCg                                                        >  2:4045225/1‑150 (MQ=255)
                                                                                                gTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGt                                                       >  2:251248/1‑150 (MQ=255)
                                                                                                 taattaatGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGtt                                                      >  2:246437/1‑150 (MQ=255)
                                                                                                   attaatGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTtct                                                    >  2:1611563/1‑150 (MQ=255)
                                                                                                       atGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTcgcg                                                >  1:1477327/1‑150 (MQ=255)
                                                                                                         ggCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCa                                              >  1:2260605/1‑150 (MQ=255)
                                                                                                         ggCCCCGGAGCGGGGTTGCGCCCCTCACCCTAATCCTCTCCCCCTAGGGGGGCGGGCCCTGCCCGGGCGGTTTACAGGGGGGGCGGGCTTAGTTCGGGGCCGTCCCGACCCATGTGTTGGAGAGTTCGGTGTGAGTGCCGTTCTCgcgcg                                              >  1:1241832/1‑149 (MQ=255)
                                                                                                          gCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAg                                             >  2:2996824/1‑150 (MQ=255)
                                                                                                          gCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCAGAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATGTCGTTGTAAGTTCCGTTCTCGCGCAg                                             <  1:581712/150‑1 (MQ=255)
                                                                                                            cccGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCAACCCCAAGGGGAGAGGGAACAGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTATATTTCGTTGTAAGTTACGTTCTCGCGCAcgg                                           >  1:382797/1‑150 (MQ=255)
                                                                                                             ccGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGt                                          >  2:1217424/1‑150 (MQ=255)
                                                                                                             ccGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGt                                          >  2:3505231/1‑150 (MQ=255)
                                                                                                                gAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGtttt                                       <  1:1285944/150‑1 (MQ=255)
                                                                                                                 aGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTc                                      <  1:1932639/150‑1 (MQ=255)
                                                                                                                     ggTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAAcg                                  <  1:251258/150‑1 (MQ=255)
                                                                                                                              cccTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGGAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGa                         >  2:1194529/1‑150 (MQ=255)
                                                                                                                               ccTCACCCTAATCCTCTCCCCAGAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGAc                        <  2:431759/150‑1 (MQ=255)
                                                                                                                                  cACCCTAAGCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACttt                     <  1:229482/150‑1 (MQ=255)
                                                                                                                                   aCCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTg                    <  2:382797/150‑1 (MQ=255)
                                                                                                                                   aCCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTg                    <  1:2406514/150‑1 (MQ=255)
                                                                                                                                    cccTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGt                   <  2:300771/150‑1 (MQ=255)
                                                                                                                                        aaTCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACg               <  1:241188/150‑1 (MQ=255)
                                                                                                                                          tCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCa             >  2:3174426/1‑150 (MQ=255)
                                                                                                                                          tCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGATTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCa             <  1:3423103/150‑1 (MQ=255)
                                                                                                                                          cccTCTCCCCTATGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCa             >  2:1722102/2‑150 (MQ=255)
                                                                                                                                            ctctCCCCTATGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGc           >  2:147588/1‑150 (MQ=255)
                                                                                                                                            ctctCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGc           >  1:2886695/1‑150 (MQ=255)
                                                                                                                                            ctctCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGc           >  2:1605490/1‑150 (MQ=255)
                                                                                                                                             tctcCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCt          >  2:2266725/1‑150 (MQ=255)
                                                                                                                                               tcCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcg        <  2:345040/150‑1 (MQ=255)
                                                                                                                                               tcCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcg        >  2:854213/1‑150 (MQ=255)
                                                                                                                                                 cccATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgtc      >  1:2633453/1‑150 (MQ=255)
                                                                                                                                                  ccATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgtcg     >  1:4073600/1‑150 (MQ=255)
                                                                                                                                                  ccATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgtcg     >  1:3971752/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTA  >  NC_000913/847065‑847363

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: