Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 847,214 T→A 6.0% 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 6.0% 388.2 / 11.5 117intergenic (‑95/+44)glnP/glnHglutamine transporter subunit/glutamine transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (58/52);  new base A (3/4);  total (61/56)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

Only 100 of 118 total aligned reads displayed.

ACCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGC  >  NC_000913/847066‑847361
                                                                                                                                                    |                                                                                                                                                    
aCCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGGCAAACTGCATATCTTGTTACTGTTACCGTGGTATCAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGCTTTGCGCCCCTCACCCGAATCCGCTCCCCAt                                                                                                                                                      >  3:797794/1‑149 (MQ=255)
  cTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATa                                                                                                                                                     <  1:64061/148‑1 (MQ=255)
   ttCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAgg                                                                                                                                                  >  4:885385/1‑149 (MQ=255)
   ttCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCATGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAgg                                                                                                                                                  >  3:551001/1‑149 (MQ=255)
    tCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAggg                                                                                                                                                 >  3:139186/1‑149 (MQ=255)
    tCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAgg                                                                                                                                                  <  1:553687/148‑1 (MQ=255)
      aaTCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAgggg                                                                                                                                                <  1:724793/148‑1 (MQ=255)
        tCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGgaga                                                                                                                                             <  3:765145/149‑1 (MQ=255)
        tCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGgag                                                                                                                                              <  1:343139/148‑1 (MQ=255)
         cAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTATGGGgaga                                                                                                                                             <  1:824863/148‑1 (MQ=255)
         cAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGgaga                                                                                                                                             >  1:157717/1‑148 (MQ=255)
         cAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGgaga                                                                                                                                             <  2:337956/148‑1 (MQ=255)
         cAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGgaga                                                                                                                                             >  1:26546/1‑148 (MQ=255)
          aGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGagag                                                                                                                                            >  1:447892/1‑148 (MQ=255)
             aGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGa                                                                                                                                         <  1:117435/148‑1 (MQ=255)
              gCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGaa                                                                                                                                        >  1:939408/1‑148 (MQ=255)
               cGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTATGGGGAGAggatta                                                                                                                                       <  2:360798/148‑5 (MQ=255)
               cGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGGTCGTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTATGGGGAGAggattag                                                                                                                                      <  3:302627/149‑6 (MQ=255)
                ggAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAACCCTCTCCCCTATGGGGAGAggattag                                                                                                                                      <  2:137370/148‑6 (MQ=255)
                  aaTGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGcc                                                                                                                                   >  3:697581/1‑149 (MQ=255)
                         ggCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGt                                                                                                                            <  3:198699/149‑1 (MQ=255)
                           ccAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGtt                                                                                                                           >  1:136260/1‑148 (MQ=255)
                              gATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTa                                                                                                                        >  2:575698/1‑148 (MQ=255)
                                  gCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTTTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAgg                                                                                                                    <  2:274836/148‑1 (MQ=255)
                                       ccAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAg                                                                                                               >  2:855423/1‑148 (MQ=255)
                                           tCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGtta                                                                                                          >  4:559147/1‑149 (MQ=255)
                                               aCTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGATAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGttatta                                                                                                       <  1:137568/148‑1 (MQ=255)
                                                  gCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTAttt                                                                                                    <  2:123112/148‑1 (MQ=255)
                                                    atatGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCgg                                                                                                 >  3:561431/1‑149 (MQ=255)
                                                     tatGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTGATCCTCTCACCATAGGGGAGAGGGAACTGCCAGTGCGTTTTGCAGGTGTAGCGTTATTATTTCgg                                                                                                 >  2:741458/1‑148 (MQ=255)
                                                     tatGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGGACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGt                                                                                                >  3:855958/1‑149 (MQ=255)
                                                           gttCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTa                                                                                          >  3:624110/1‑149 (MQ=255)
                                                            ttCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTAc                                                                                         <  3:548096/149‑1 (MQ=255)
                                                            ttCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTAc                                                                                         >  3:543683/1‑149 (MQ=255)
                                                               cTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGa                                                                                      >  4:568727/1‑149 (MQ=255)
                                                                tGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGaa                                                                                     >  3:549506/1‑149 (MQ=255)
                                                                tGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGa                                                                                      >  2:679129/1‑148 (MQ=255)
                                                                 gTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGaa                                                                                     >  2:544499/1‑148 (MQ=255)
                                                                     ccGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCa                                                                                                                                  <  4:132824/99‑1 (MQ=255)
                                                                     ccGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCa                                                                                                                                  >  3:132824/1‑99 (MQ=255)
                                                                      cGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCGGTACTGAAccgaaa                                                                               <  3:676850/149‑5 (MQ=255)
                                                                       gtggtgTAAATCATCAAATTCAGGGTAATTGATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAAGGGAGAGGGAACTGCCAGTGCGTTTTAAAGGTGTAGCGTTATTATTTCGGTGCTGAACCGaaataataa                                                                              <  4:797794/149‑8 (MQ=255)
                                                                         ggtgTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCAtttttt                                                                            <  3:347961/149‑1 (MQ=255)
                                                                          gtgtAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTg                                                                           <  3:544617/149‑1 (MQ=255)
                                                                           tgtAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGt                                                                          >  3:384963/1‑149 (MQ=255)
                                                                                   atCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCg                                                                  <  3:460084/149‑1 (MQ=255)
                                                                                    tCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGt                                                                 >  3:928758/1‑149 (MQ=255)
                                                                                     cAAATTCAGGGTAATTAATGGCCCCGGTGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGt                                                                 >  2:880548/1‑148 (MQ=255)
                                                                                     cAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGCGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTAATTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGtt                                                                <  4:231167/149‑1 (MQ=255)
                                                                                     cAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGtt                                                                <  3:126973/149‑1 (MQ=255)
                                                                                     cAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGtt                                                                <  4:697581/149‑1 (MQ=255)
                                                                                     cAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGtt                                                                <  3:722313/149‑1 (MQ=255)
                                                                                         ttCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTaa                                                            <  4:679945/149‑1 (MQ=255)
                                                                                         ttCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTa                                                             <  1:756255/148‑1 (MQ=255)
                                                                                             gggTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTCACAGGTGTAGCGTTATTATTTCGGATCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTc                                                        <  4:551001/149‑1 (MQ=255)
                                                                                              ggTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTcc                                                       <  3:358007/149‑1 (MQ=255)
                                                                                              ggTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTc                                                        <  2:437081/148‑1 (MQ=255)
                                                                                              ggTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTc                                                        >  1:349816/1‑148 (MQ=255)
                                                                                               gTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTcc                                                       >  1:255330/1‑148 (MQ=255)
                                                                                                        ggCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCgcgc                                             >  3:720017/1‑149 (MQ=255)
                                                                                                           cccGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAgg                                          >  4:351648/1‑149 (MQ=255)
                                                                                                            ccGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAgg                                          >  1:673253/1‑148 (MQ=255)
                                                                                                             cGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATCTCGTTGTAAGTTCCGTTCTCGCGCAggg                                         >  1:554178/1‑148 (MQ=255)
                                                                                                                aGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAgg                                                                                                                    >  1:23322/1‑70 (MQ=255)
                                                                                                                aGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAgg                                                                                                                    <  2:23322/70‑1 (MQ=255)
                                                                                                                    ggTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCaa                                  >  1:802237/1‑148 (MQ=255)
                                                                                                                    ggTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAAc                                 <  3:351648/149‑1 (MQ=255)
                                                                                                                               cTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCATTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGa                       >  2:605598/1‑148 (MQ=255)
                                                                                                                               cTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTg                                                               <  3:313410/108‑1 (MQ=255)
                                                                                                                               cTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTg                                                               >  4:313410/1‑108 (MQ=255)
                                                                                                                               cTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGa                       <  2:941085/148‑1 (MQ=255)
                                                                                                                                tCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGATCCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACt                     >  3:557337/1‑149 (MQ=255)
                                                                                                                                  aCCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACtt                    <  2:157717/148‑1 (MQ=255)
                                                                                                                                   cccTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCAtttttt                                                                            >  4:790341/1‑91 (MQ=255)
                                                                                                                                    ccTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGtttt                                     <  4:230973/129‑1 (MQ=255)
                                                                                                                                    ccTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAAcgc                               >  3:350572/1‑135 (MQ=255)
                                                                                                                                    ccTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAAcgc                               <  4:350572/135‑1 (MQ=255)
                                                                                                                                    ccTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTg                  <  1:127133/148‑1 (MQ=255)
                                                                                                                                     cTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGt                 >  1:149390/1‑148 (MQ=255)
                                                                                                                                     cTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGt                 >  2:868215/1‑148 (MQ=255)
                                                                                                                                     cTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTc                <  3:568727/149‑1 (MQ=255)
                                                                                                                                      tactCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCa               <  3:885385/146‑1 (MQ=255)
                                                                                                                                      tAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTc                <  2:759110/148‑1 (MQ=255)
                                                                                                                                      tAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCa               >  4:632941/1‑149 (MQ=255)
                                                                                                                                       aaTCCCCTACCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAATCATTTTTTGTAGATTACGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCa               <  2:421759/148‑1 (MQ=255)
                                                                                                                                          ccTCTCCCCTATGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGc            >  2:140113/1‑148 (MQ=255)
                                                                                                                                           ctctCCCCTATGGGGAGAGGGAACTGCCAGTGCGTTCTACAGGTGTAGCGTTATTATTTCGGTACAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGAGTTTTCAACGCGCCGTTGACTTTCTCACGCAg          >  3:596663/1‑149 (MQ=255)
                                                                                                                                            tctcCCCTATGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTc                                                        <  2:16993/102‑1 (MQ=255)
                                                                                                                                            tctcCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTGACGCAg          >  1:938255/1‑148 (MQ=255)
                                                                                                                                            tctcCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGc         >  3:277782/1‑149 (MQ=255)
                                                                                                                                                cccATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgt     >  3:349152/1‑149 (MQ=255)
                                                                                                                                                cccATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgt     >  3:393280/1‑149 (MQ=255)
                                                                                                                                                cccATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgt     <  3:559147/149‑1 (MQ=255)
                                                                                                                                                 ccATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgtc    >  3:155456/1‑149 (MQ=255)
                                                                                                                                                 ccATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgtc    <  4:519084/149‑1 (MQ=255)
                                                                                                                                                 ccATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgtc    >  3:334455/1‑149 (MQ=255)
                                                                                                                                                 ccATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgtc    <  4:651220/149‑1 (MQ=255)
                                                                                                                                                 ccATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgt     >  1:244724/1‑148 (MQ=255)
                                                                                                                                                   aTAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgtcg   >  2:500656/1‑148 (MQ=255)
                                                                                                                                                   aTAGGGGAGAGGGAACTGCCAGTGCGTTTTACAAGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGc  >  4:124034/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGC  >  NC_000913/847066‑847361

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: