Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,228,2750TC43.3% 101.4 / 16.9 99pseudogene (335/1287 nt)mdtQpseudogene; putative channel/filament proteins
Reads supporting (aligned to +/- strand):  ref base T (40/15);  new base C (12/30);  total (53/46)
Fisher's exact test for biased strand distribution p-value = 2.98e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGC  >  NC_000913/2228126‑2228421
                                                                                                                                                          |                                                                                                                                                  
gCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCTc                                                                                                                                                   <  1:691775/151‑1 (MQ=255)
        tCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCT‑CGCGTTCCCCCGCCC‑GGGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCCCCCCCCCAt                                                                                                                                           >  1:1670381/1‑151 (MQ=255)
         cTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCT‑CGCGGTCCCCCGCCC‑CTGG‑TTTTACCGT‑A‑CCCAGGCGGGCCGTAACCTCCCCCCCAtt                                                                                                                                          >  1:1174545/1‑151 (MQ=255)
           aTTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGGTTGCCACTCCCAGTACAGGCGGGCTACGCGCCCAACCAGCAATTGGCGGGGTTTCT‑TGCCGTTCGCCGCCC‑CTCC‑CTTTAACGT‑A‑ACCAGCCCGGCCGGAAACTcccgcc                                                                                                                                              >  2:55434/1‑141 (MQ=255)
             ttGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCtt                                                                                                                                       <  2:1748257/150‑1 (MQ=255)
             ttGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCT‑CGCGGTTCGCCGCCC‑CGGC‑TTTTACCCT‑A‑CCCCGGCGGGGAGTAACCCCCCCCCCATTCttt                                                                                                                                      >  2:1219182/1‑151 (MQ=255)
               gctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGGT‑CGCGGTCCGCCGCCC‑CTGC‑TTTAACCGT‑A‑CCCCGGGGGGCAGGAACCCCCCCCCGGTTTTTTcc                                                                                                                                    >  1:1110392/1‑151 (MQ=255)
                   cAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCT‑CGCGTTTCGCCGCCC‑GTGC‑TTTTAACCT‑A‑CCCAGGCGGGCAGTAAACTCCCCCCGATTCTTTtccc                                                                                                                                  >  1:579478/1‑149 (MQ=255)
                    aaGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGa                                                                                                                               <  1:828674/151‑1 (MQ=255)
                    aaGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCT‑CGCGTTCCGCCGCCC‑GTGG‑TTTAACCGT‑A‑CCCAAGCGGGCAGTAACCCCCCCCCCGTTCTTTtccc                                                                                                                                  >  1:1037228/1‑148 (MQ=255)
                     aGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGat                                                                                                                              <  1:593464/151‑1 (MQ=255)
                      gACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGata                                                                                                                             <  1:1703866/151‑1 (MQ=255)
                      gACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTTACCGT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTTCCCAGa                                                                                                                               >  1:1417360/1‑149 (MQ=255)
                       aCCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGatat                                                                                                                            <  1:421437/151‑1 (MQ=255)
                        ccGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCT‑CGCGGTCCGCCGCCC‑CGGC‑TTTTACCCTAA‑CCCAGCGGGGCGGAAACCCCCCCCCGCTTTTTTccc                                                                                                                                   >  2:535143/1‑144 (MQ=255)
                         cGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCT‑CCCGGTCCCCCGCCC‑CGGCTTTTAACCCA‑C‑CCCAGG‑GGGGAGGAAACTCCCCCCCCTTCTTTtcccccata                                                                                                                             >  2:111397/1‑143 (MQ=255)
                            gtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCT‑CCCGGTTCCCCGCCC‑CGGC‑TTTAACCGT‑ACCCCCGGGGGGCGGAAACCCCCCCCCAATTTTTTccc                                                                                                                                   >  2:326432/1‑140 (MQ=255)
                            gtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCT‑CGCGTTTCGCCCCCC‑GGGC‑TTTTACCCT‑A‑CCCAGGCGGGCAGTAAACTCCCCCCGCTTCTTTTCCCCGata                                                                                                                             >  1:760218/1‑145 (MQ=255)
                             tttAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCT‑CGCGGTTCGCCGCCC‑CCGC‑TTTTAACGT‑A‑CCCCGGCGGGGAGGAACCTCCGCCCGGTTCTTTTCCCCGATATCgag                                                                                                                        >  2:96434/1‑149 (MQ=255)
                              ttAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCT‑CGCGTTCCCCCCCCC‑GTGC‑TTTAAACGG‑A‑CCCAGGCGGGCAGTAACCTCCCCCCCATTCTTTCCCCCGATATCgaga                                                                                                                       >  1:1009611/1‑149 (MQ=255)
                               tAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCT‑CCCGGTCCGCCGCCC‑GTGC‑TTTAACCGG‑A‑CCCAGGCGGGCAGTAAACTCCCCCCCATTCTTTCCCCAGATATCGAGATGc                                                                                                                    >  1:378340/1‑151 (MQ=255)
                                  cgccgcCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCC‑CGCGTTCCCCCCCCC‑CTGC‑TTTTAACGT‑A‑CCCCGGCGGGCAGTAAACCCCCCCCGATTCTTTCCCCCGATATTGAGATGcc                                                                                                                   >  2:534004/1‑149 (MQ=255)
                                  cgccgcCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCT‑CGCGGTCCGCCCCCC‑CTGC‑TTTTAACCT‑A‑CCCCAGCGGGGAGGAAACTCCCCCCGATTCTTTTCCCCGATAACGAGAAGcc                                                                                                                   >  1:408558/1‑149 (MQ=255)
                                      gccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCTCCCCCCGATTCTTTTCCCCGATATCGAGATGCCAGccc                                                                                                              >  1:874329/1‑150 (MQ=255)
                                       ccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCT‑CCCGCTTCCCCGCCC‑GTGC‑TTTTACCGC‑A‑CCCCGGCGGGCAGGAAACCCCCCCCGATTCTTTCCCCCAATATCGAGATTCCAGccc                                                                                                              >  1:312294/1‑149 (MQ=255)
                                        cTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCC‑CGCGTTCCGCCCCCCGGGGT‑TTTAACCTA‑A‑CCCGGGGGGGC‑GTGACCCCCCCCCCCTTTTTTTCCCCCATATCCAGATGGCAGcccccc                                                                                                           >  2:1411038/1‑151 (MQ=255)
                                        cTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCT‑CCCGTTCCCCCCCCC‑GTGC‑TTTAACCCT‑A‑CCCAGGCGGGCCGGAACCTCCCCCCCCTTCTTTCCCCCGATATCGAGATGCCCGcccccc                                                                                                           >  1:255067/1‑151 (MQ=255)
                                         tGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCT‑CGCGGTTCCCCCCCC‑CTGCTTTTAACCTA‑A‑CCCAGCGGGGCGGAAACC‑CCCCCCCATTCTTTTCCCCCAAATCCAGGTGcccccccccc                                                                                                           >  2:731110/1‑142 (MQ=255)
                                          gggTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTACCCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGt                                                                                                         <  1:774169/151‑1 (MQ=255)
                                          gggTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCccgccg                                                                                                          <  2:392142/150‑1 (MQ=255)
                                            gTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCT‑CGCGGTCCGCCCCCC‑CTGC‑CTTTACCCT‑A‑CCCCGGCGGGCAGTAAACCCCGCCCCCTTTTTTTCCCCGATATCGAGATGCCCACcccccgct                                                                                                        >  2:172762/1‑148 (MQ=255)
                                             tttGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCT‑CGCGGTCCCCCGCCC‑CTGC‑TTTTAACGG‑A‑CCCAGGGGGGGAGGAAACCCCCCCCCGTTTTTTCCCCCGATATAGAGATGGCAACccgcccct                                                                                                        >  2:64523/1‑146 (MQ=255)
                                                 ccccTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTACCCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGc                                                                                                    <  1:927249/146‑1 (MQ=255)
                                                  cACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTAACCGT‑A‑CCCCGGCGGGCAGTAACCTCCCCCCGATTCTTTTCCCAGATATCGAGATGCCAGCCCGCCGTTAAGcc                                                                                                   >  1:692790/1‑149 (MQ=255)
                                                    cTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCT‑CCCGTTTCCCCGCCC‑CTGC‑TTTTAACCC‑A‑CCCCGGCGGGCCGTAACCCCCCCCCCATTCTTTCCCCAGATATCGAGATGCCAGCCCCCCCTTAAGCCcaa                                                                                                >  2:1321981/1‑150 (MQ=255)
                                                       ccAGTACAGGCGGGTCCCGCTCCCAGCCACAAATGGGGGGTTTTGCT‑GGCTTCCCGCCCCCC‑GTGC‑TTAACCCGT‑A‑CCCAGGCGGGAAGAACCCCCCGCCCGATTCTTTCCCCAGATACCGAGAGCCCACCCCGCCGTTAAGCCAAAAGTa                                                                                            <  2:408933/151‑1 (MQ=255)
                                                        cAGTACAGGCGGGCTACGCTGCCAGCCAGCATTTGGCGGTTTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTACCCGT‑C‑CCCAGGCGGGCAGTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTAc                                                                                           <  2:375179/151‑1 (MQ=255)
                                                         aGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTAACCGT‑C‑CCCAGGCGGGCAGTACCCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTAc                                                                                           <  1:1507232/150‑1 (MQ=255)
                                                          gtaggggggggCGACGCTGCCAGCCAGCAATTGGGGGTTTTGCT‑CGCTTCCCGCCGCCC‑GTGC‑TTTAACCGT‑C‑CCCGGGGGGGCAATACCCCCCGCCCGATTCTTCCCCCAGATATCGAGTGGCCAGCCCGCCGTTAAGCCAAAAGTACCg                                                                                         <  2:968162/143‑1 (MQ=255)
                                                          gTACGGGGGGGCTCCGCGCCCAGCCAGCATTTGGCGGTTTTGCC‑CGCGTCCCCCCGCCC‑GTGC‑TTAACCCGT‑C‑CCCGGGGGGGCAGTACCCCCCCCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTAcc                                                                                          <  2:1234723/150‑1 (MQ=255)
                                                          gTACCGGCGGTCTACTGTGCGAGCGCGAATTTGGGGGTTTCGCG‑CGCGTCCGCCCGCCC‑GTGT‑TTTAGCGGG‑A‑CCTGGCTGGTCACTGCCCGCCGCTCGATCCTTGCCCGAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTAc                                                                                           <  2:172109/149‑1 (MQ=255)
                                                          gTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCC‑GTGT‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCg                                                                                         <  1:1613468/151‑1 (MQ=255)
                                                            aCAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTAAACCT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCCATTTTTTCCCCCGATATCGAGATGCCAGCCCGCCCTTAAACCAAAAATAcc                                                                                          >  1:1204880/1‑148 (MQ=255)
                                                              aGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCCATTCTTTCCCCCGATATCGAGATGCCAGCCCGCCGTTAAGCCCAAAGTACCGtt                                                                                       >  2:685190/1‑149 (MQ=255)
                                                               ggggggCTACGCTGCCACCCAGCATTGGGCGGGTTGGCT‑CGCGTCCCGCCGCCC‑GTGT‑TTTACCCGT‑C‑CCCAGGCGGGCAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTgg                                                                                     <  2:1134041/147‑1 (MQ=255)
                                                               ggCGGGCTCCGCTCCCAGCCAGCAATTGGCGGTTTTGCC‑CGCTTCCCGCCCCCC‑TTCT‑TTTACCCTT‑C‑CCCAGGGGGGCAGTACCCCCCCCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTAcc                                                                                          <  2:703705/145‑1 (MQ=255)
                                                               ggCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTgg                                                                                     >  1:466829/1‑150 (MQ=255)
                                                                  ggtctccgcgGCCAGCCAGCAATGGGCGGTTTTGCT‑CGCTTCCCGCCCCCC‑GTGC‑TTAACCCGT‑C‑CCCAGGGGGGCAGAACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTa                                                                                 <  2:1271104/145‑1 (MQ=255)
                                                                  gggCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCC‑GTGT‑TTTAACCGT‑C‑CCCAGGCGGGCAGTACCCTCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTa                                                                                 <  1:508386/151‑1 (MQ=255)
                                                                  gggCTACGCGCCCAGCCAGCAATGGGCGGTTTTGCC‑CGCGTCCCGCCCCCG‑GTGC‑TTTACCCGT‑C‑CCCGGGGGGGCAGTACCCCCCGCCCGTTTCTTCCCCCAGATATCGAGATCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTa                                                                                 <  2:83641/151‑1 (MQ=255)
                                                                     cTACGCTGCCAGCCAGCAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTAACCGG‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGGTAAGCCAAAAGTACCGTTGGGGTAc                                                                                >  2:211775/1‑149 (MQ=255)
                                                                      tACGCTGCCAGCCAGCAATTGGCGGGTTTGCT‑CGCGGTCCGCCGCCC‑GTTC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTTCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCACAAGGACCGTTGGTGGACCAc                                                                             >  2:1384428/1‑151 (MQ=255)
                                                                       ccGCTGCCACCCAGCAATGGGCGGGTTTCCT‑CGCGTTCCGCCCCCC‑GTGC‑TTTACCCGT‑C‑CCCGGGGGGGCAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACg                                                                            <  2:1738873/150‑1 (MQ=255)
                                                                           tGCCAGCCAGCAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGcc                                                                        >  1:699715/1‑151 (MQ=255)
                                                                             ccagccagCAATTGGCGGTTTTGCT‑CGCGTCCCGCCGCCC‑GTGC‑TTTACCCGT‑A‑CCCAGCCGGGCAGTAACCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGcc                                                                        <  1:806108/149‑1 (MQ=255)
                                                                             ccagccagCAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTAACCGT‑A‑CCCCGGCGGGGAGTAACCTCCCCCCGATTCTTTCCCCAGAGATCGAGATTCCAGCCCGCCGTTAAGCCAAAAGGACCGTTTGTGGACCACGGGCCgg                                                                      >  2:356373/1‑151 (MQ=255)
                                                                             ccagccagCAATTGGCGGGTTTGCC‑CGCGTCCCGCCCCCC‑TTGC‑TTTACCCGT‑A‑CCCAGGCGGGCAGTACCCCCCCCCCGATTCTTTCCCCAGATATCGAGATGCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGcc                                                                        <  2:1692013/149‑1 (MQ=255)
                                                                              cagccagcAATTGGCGGGTTTATTCCCCTTTCTGCCG‑CC‑GCGC‑TTTTACCCG‑A‑CCCAGCCGGGCGTTAACCTCTGACCCCTTCTTTACCCCGGTGTCGAGATGCCCGCCCGGCAGTTCGCCAAAAGTaatgtaggggtaaca                                                                              >  2:1477657/1‑128 (MQ=255)
                                                                               agccagcAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCgg                                                                      >  1:431706/1‑149 (MQ=255)
                                                                                 ccagcAATTGGCGGGTTTGCT‑CGCGTCCCGCCGCCC‑GTGC‑TTTACCCGT‑A‑CCCAGGCGGGCAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGt                                                                  <  1:1622762/151‑1 (MQ=255)
                                                                                  cagcAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGt                                                                  >  1:225337/1‑150 (MQ=255)
                                                                                    gcAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCTCCCCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCCCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTAcc                                                               >  1:954221/1‑151 (MQ=255)
                                                                                     cAATTGGCGGGTTTGCT‑CGCGTTCCGCCCCCC‑GTGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAATACCGTTGGTGTACCACGGGCCGGTCGGAcc                                                               >  2:1726401/1‑150 (MQ=255)
                                                                                      atttGGCGGGTTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTACCCGT‑C‑CCCAGGCGGGCAGTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTg                                                             <  1:211781/149‑1 (MQ=255)
                                                                                      aaTTGGGGGTTTTGCT‑CGCGTTCCGCCCCCC‑GTGC‑TTTAACCGT‑C‑CCCAGGCGGGCAGTACCCCCCGCCCGATTCTTCCCCAAGATATCGAGATCCCGCCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTc                                                                    <  2:149505/144‑1 (MQ=255)
                                                                                      aaTTGGCGGGTTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTACCCGT‑A‑CCCAGGCGGGCAGTACCCTCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTg                                                             <  1:702737/151‑1 (MQ=255)
                                                                                      aaTTGGCGGGTTTGCT‑CGCGTCCCGCCGCCC‑GTGC‑TTTACCCGT‑C‑CCCAGGCGGGCAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTg                                                             <  2:1417247/151‑1 (MQ=255)
                                                                                        ttGGCGGGTTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCg                                                           >  2:1523255/1‑151 (MQ=255)
                                                                                                          gTTCCGCCGCCC‑GTGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGc                                          >  2:561294/1‑151 (MQ=255)
                                                                                                           tcccgccgCCC‑GTGT‑TTTAACCGT‑C‑CCCGGGCGGGCAGTACCCCCCCCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGc                                          <  1:1327547/148‑1 (MQ=255)
                                                                                                           tcccgccgCCC‑GTGC‑TTTAACCGC‑C‑CCCAGGCGGGCAGTACCCCCCCCCCGATTCTTCCCCCAGATAGCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGc                                          <  2:1012237/148‑1 (MQ=255)
                                                                                                           tcccgcccCCC‑GTGC‑TTTACCCGT‑C‑CCCAGGGGGGAAGAACCCCCCCCCCGATTTTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCaga                                            <  2:135913/146‑1 (MQ=255)
                                                                                                            tCCGCCGCCC‑GTGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCa                                         >  1:731110/1‑150 (MQ=255)
                                                                                                             ccgccgCCC‑GTCA‑TTAACCCGT‑A‑CCCAGGCGGGCAGAACCCCCCCCCCAATTCTTCCCCCAGATATCGATAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCa                                         <  2:1760417/149‑1 (MQ=255)
                                                                                                              cgccgcCC‑GTGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAAc                                      >  2:304634/1‑151 (MQ=255)
                                                                                                                cccccc‑GTGC‑TTTACCCGT‑C‑CCCGGGGGGGCAGTACCCCCCGCCCGTTTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACg                                     <  2:1605843/147‑1 (MQ=255)
                                                                                                                        tGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGCTTCTTTCCCCCGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCAtt                             >  1:872518/1‑151 (MQ=255)
                                                                                                                            tttACCCGT‑A‑CCCAGGCGGGCAGTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTaa                           <  1:809651/150‑1 (MQ=255)
                                                                                                                            tttAACCGT‑C‑CCCAGGCGGGAAGTACCCTCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTa                            <  1:1021385/149‑1 (MQ=255)
                                                                                                                             tttcccGT‑C‑CCCAGGCGGGCAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAg                          <  1:1166234/146‑1 (MQ=255)
                                                                                                                              tAACCGT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGcc                        >  1:106294/1‑151 (MQ=255)
                                                                                                                               acccGT‑A‑CCCGGGCGGGCAGTAACCCCCGCCCGATTCTTCCCCCAGAAAGCGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGcc                        <  1:491621/148‑1 (MQ=255)
                                                                                                                                 ccGT‑C‑CCCAGGCGGGCAGTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTc                     <  1:6177/151‑1 (MQ=255)
                                                                                                                                 ccGT‑A‑CCCAGGCGGGCAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTc                     <  1:957567/151‑1 (MQ=255)
                                                                                                                                 ccGT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTc                     >  2:751904/1‑151 (MQ=255)
                                                                                                                                   gT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTg                   >  1:836264/1‑151 (MQ=255)
                                                                                                                                   gT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGAGCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTg                   >  2:1343121/1‑151 (MQ=255)
                                                                                                                                    t‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGc                  >  1:847290/1‑151 (MQ=255)
                                                                                                                                        cccGGGGGGGCAGTACCCCCCGCCCGTTTCTTCCCCCAGATACCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTg                   <  2:582656/148‑1 (MQ=255)
                                                                                                                                        cccAGGCGGGCAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCt                    <  1:1219261/147‑1 (MQ=255)
                                                                                                                                         ccAGGCGGGCAGTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGa               <  1:411923/151‑1 (MQ=255)
                                                                                                                                          cAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGAc              >  1:411075/1‑151 (MQ=255)
                                                                                                                                             gCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAtt           >  2:1778597/1‑151 (MQ=255)
                                                                                                                                              cGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAttt          >  1:529410/1‑151 (MQ=255)
                                                                                                                                               gggCAGTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAttt          <  2:760183/150‑1 (MQ=255)
                                                                                                                                               gggCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAtttt         >  1:1186682/1‑151 (MQ=255)
                                                                                                                                               gggAAGTACCCCCCGCCCGTTTCTCCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAttt          <  2:1037175/150‑1 (MQ=255)
                                                                                                                                                      aaCCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGc  >  2:1592862/1‑151 (MQ=255)
                                                                                                                                                        ccTCCCCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGc  <  1:1386256/149‑1 (MQ=255)
                                                                                                                                                          |                                                                                                                                                  
GCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCC‑GTGC‑TTTAACCGT‑A‑CCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGC  >  NC_000913/2228126‑2228421

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