Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 847,215 A→T 5.9% intergenic (‑96/+43) 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,2150AT 5.9% 383.2 / 10.0 115intergenic (‑96/+43)glnP/glnHglutamine transporter subunit/glutamine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (58/50);  new base T (3/4);  total (61/54)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 118 total aligned reads displayed.

CCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTA  >  NC_000913/847067‑847363
                                                                                                                                                    |                                                                                                                                                     
 cTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATa                                                                                                                                                       <  1:64061/148‑1 (MQ=255)
  ttCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑gg                                                                                                                                                    >  4:885385/1‑149 (MQ=255)
  ttCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCATGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑gg                                                                                                                                                    >  3:551001/1‑149 (MQ=255)
   tCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑ggg                                                                                                                                                   >  3:139186/1‑149 (MQ=255)
   tCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑gg                                                                                                                                                    <  1:553687/148‑1 (MQ=255)
     aaTCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑gggg                                                                                                                                                  <  1:724793/148‑1 (MQ=255)
       tCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGgaga                                                                                                                                               <  3:765145/149‑1 (MQ=255)
       tCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGgaga                                                                                                                                               >  4:23091/1‑149 (MQ=255)
       tCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGgag                                                                                                                                                <  1:343139/148‑1 (MQ=255)
        cAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTAT‑GGGgaga                                                                                                                                               <  1:824863/148‑1 (MQ=255)
        cAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGgaga                                                                                                                                               >  1:157717/1‑148 (MQ=255)
        cAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGgaga                                                                                                                                               <  2:337956/148‑1 (MQ=255)
        cAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGgaga                                                                                                                                               >  1:26546/1‑148 (MQ=255)
         aGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGagag                                                                                                                                              >  1:447892/1‑148 (MQ=255)
            aGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGa                                                                                                                                           <  1:117435/148‑1 (MQ=255)
             gCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGaa                                                                                                                                          >  1:939408/1‑148 (MQ=255)
              cGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCC‑TATGGGGAGAggatta                                                                                                                                         <  2:360798/148‑5 (MQ=255)
              cGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGGTCGTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTAT‑GGGGAGAggattag                                                                                                                                        <  3:302627/149‑6 (MQ=255)
               ggAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAACCCTCTCCCC‑TATGGGGAGAggattag                                                                                                                                        <  2:137370/148‑6 (MQ=255)
                 aaTGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGcc                                                                                                                                     >  3:697581/1‑149 (MQ=255)
                        ggCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGt                                                                                                                              <  3:198699/149‑1 (MQ=255)
                          ccAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGtt                                                                                                                             >  1:136260/1‑148 (MQ=255)
                             gATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTa                                                                                                                          >  2:575698/1‑148 (MQ=255)
                                 gCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTTTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAgg                                                                                                                      <  2:274836/148‑1 (MQ=255)
                                      ccAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAg                                                                                                                 >  2:855423/1‑148 (MQ=255)
                                          tCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGtta                                                                                                            >  4:559147/1‑149 (MQ=255)
                                              aCTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGATAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGttatta                                                                                                         <  1:137568/148‑1 (MQ=255)
                                                 gCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTAttt                                                                                                      <  2:123112/148‑1 (MQ=255)
                                                   atatGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCgg                                                                                                   >  3:561431/1‑149 (MQ=255)
                                                    tatGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTGATCCTCTCACCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTGCAGGTGTAGCGTTATTATTTCgg                                                                                                   >  2:741458/1‑148 (MQ=255)
                                                    tatGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGGACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGt                                                                                                  >  3:855958/1‑149 (MQ=255)
                                                          gttCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTa                                                                                            >  3:624110/1‑149 (MQ=255)
                                                           ttCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTAc                                                                                           <  3:548096/149‑1 (MQ=255)
                                                           ttCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTAc                                                                                           >  3:543683/1‑149 (MQ=255)
                                                              cTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGa                                                                                        >  4:568727/1‑149 (MQ=255)
                                                               tGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGaa                                                                                       >  3:549506/1‑149 (MQ=255)
                                                               tGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGa                                                                                        >  2:679129/1‑148 (MQ=255)
                                                                gTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGaa                                                                                       >  2:544499/1‑148 (MQ=255)
                                                                    ccGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCa                                                                                                                                    <  4:132824/99‑1 (MQ=255)
                                                                    ccGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCa                                                                                                                                    >  3:132824/1‑99 (MQ=255)
                                                                     cGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCGGTACTGAAccgaaa                                                                                 <  3:676850/149‑5 (MQ=255)
                                                                      gtggtgTAAATCATCAAATTCAGGGTAATTGATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑AGGGAGAGGGAACTGCCAGTGCGTTTTAAAGGTGTAGCGTTATTATTTCGGTGCTGAACCGaaataataa                                                                                <  4:797794/149‑8 (MQ=255)
                                                                        ggtgTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCAtttttt                                                                              <  3:347961/149‑1 (MQ=255)
                                                                         gtgtAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTg                                                                             <  3:544617/149‑1 (MQ=255)
                                                                          tgtAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGt                                                                            >  3:384963/1‑149 (MQ=255)
                                                                                  atCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCg                                                                    <  3:460084/149‑1 (MQ=255)
                                                                                   tCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTTGATTTCGt                                                                   <  4:450239/149‑1 (MQ=255)
                                                                                   tCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGt                                                                   >  3:928758/1‑149 (MQ=255)
                                                                                    cAAATTCAGGGTAATTAATGGCCCCGGTGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGt                                                                   >  2:880548/1‑148 (MQ=255)
                                                                                    cAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGCGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTAATTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGtt                                                                  <  4:231167/149‑1 (MQ=255)
                                                                                    cAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGtt                                                                  <  3:126973/149‑1 (MQ=255)
                                                                                    cAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGtt                                                                  <  4:697581/149‑1 (MQ=255)
                                                                                    cAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGtt                                                                  <  3:722313/149‑1 (MQ=255)
                                                                                        ttCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTaa                                                              <  4:679945/149‑1 (MQ=255)
                                                                                        ttCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTa                                                               <  1:756255/148‑1 (MQ=255)
                                                                                            gggTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTCACAGGTGTAGCGTTATTATTTCGGATCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTc                                                          <  4:551001/149‑1 (MQ=255)
                                                                                             ggTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTcc                                                         <  3:358007/149‑1 (MQ=255)
                                                                                             ggTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTc                                                          <  2:437081/148‑1 (MQ=255)
                                                                                             ggTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTc                                                          >  1:349816/1‑148 (MQ=255)
                                                                                              gTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTcc                                                         >  1:255330/1‑148 (MQ=255)
                                                                                                       ggCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCgcgc                                               >  3:720017/1‑149 (MQ=255)
                                                                                                          cccGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAgg                                            >  4:351648/1‑149 (MQ=255)
                                                                                                           ccGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAgg                                            >  1:673253/1‑148 (MQ=255)
                                                                                                            cGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATCTCGTTGTAAGTTCCGTTCTCGCGCAggg                                           >  1:554178/1‑148 (MQ=255)
                                                                                                               aGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAgg                                                                                                                      >  1:23322/1‑70 (MQ=255)
                                                                                                               aGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAgg                                                                                                                      <  2:23322/70‑1 (MQ=255)
                                                                                                                   ggTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCaa                                    >  1:802237/1‑148 (MQ=255)
                                                                                                                   ggTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAAc                                   <  3:351648/149‑1 (MQ=255)
                                                                                                                              cTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCATTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGa                         >  2:605598/1‑148 (MQ=255)
                                                                                                                              cTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTg                                                                 >  4:313410/1‑108 (MQ=255)
                                                                                                                              cTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTg                                                                 <  3:313410/108‑1 (MQ=255)
                                                                                                                              cTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGa                         <  2:941085/148‑1 (MQ=255)
                                                                                                                               tCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGATCCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACt                       >  3:557337/1‑149 (MQ=255)
                                                                                                                                 aCCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACtt                      <  2:157717/148‑1 (MQ=255)
                                                                                                                                  cccTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCAtttttt                                                                              >  4:790341/1‑91 (MQ=255)
                                                                                                                                   ccTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAAcgc                                 >  3:350572/1‑135 (MQ=255)
                                                                                                                                   ccTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAAcgc                                 <  4:350572/135‑1 (MQ=255)
                                                                                                                                   ccTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTg                    <  1:127133/148‑1 (MQ=255)
                                                                                                                                    cTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGt                   >  1:149390/1‑148 (MQ=255)
                                                                                                                                    cTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGt                   >  2:868215/1‑148 (MQ=255)
                                                                                                                                    cTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTc                  <  3:568727/149‑1 (MQ=255)
                                                                                                                                     tactCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCa                 <  3:885385/146‑1 (MQ=255)
                                                                                                                                     tAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTc                  <  2:759110/148‑1 (MQ=255)
                                                                                                                                     tAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCa                 >  4:632941/1‑149 (MQ=255)
                                                                                                                                      aaTCCCCTACCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAATCATTTTTTGTAGATTACGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCa                 <  2:421759/148‑1 (MQ=255)
                                                                                                                                         ccTCTCCCCTAT‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGc              >  2:140113/1‑148 (MQ=255)
                                                                                                                                          ctctCCCCTAT‑GGGGAGAGGGAACTGCCAGTGCGTTCTACAGGTGTAGCGTTATTATTTCGGTACAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGAGTTTTCAACGCGCCGTTGACTTTCTCACGCAg            >  3:596663/1‑149 (MQ=255)
                                                                                                                                           tctcCCCTAT‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTc                                                          <  2:16993/102‑1 (MQ=255)
                                                                                                                                           tctcCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTGACGCAg            >  1:938255/1‑148 (MQ=255)
                                                                                                                                           tctcCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGc           >  3:277782/1‑149 (MQ=255)
                                                                                                                                               cccATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgt       >  3:349152/1‑149 (MQ=255)
                                                                                                                                               cccATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgt       >  3:393280/1‑149 (MQ=255)
                                                                                                                                               cccATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgt       <  3:559147/149‑1 (MQ=255)
                                                                                                                                                ccATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgtc      >  3:155456/1‑149 (MQ=255)
                                                                                                                                                ccATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgtc      >  3:334455/1‑149 (MQ=255)
                                                                                                                                                ccATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgtc      <  4:651220/149‑1 (MQ=255)
                                                                                                                                                ccATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgt       >  1:244724/1‑148 (MQ=255)
                                                                                                                                                  aTA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgtcg     >  2:500656/1‑148 (MQ=255)
                                                                                                                                                  aTA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAAGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGc    >  4:124034/1‑149 (MQ=255)
                                                                                                                                                    a‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTa  >  4:636981/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATA‑GGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTA  >  NC_000913/847067‑847363

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