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

GCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGC  >  NC_000913/2228126‑2228421
                                                                                                                                                          |                                                                                                                                                  
gCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTc                                                                                                                                                   <  1:691775/151‑1 (MQ=255)
        tCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGTTCCCCCGCCCGGGCTTTAACCGTACCCAGGCGGGCAGTAACCCCCCCCCCAt                                                                                                                                           >  1:1670381/1‑151 (MQ=255)
         cTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGGTCCCCCGCCCCTGGTTTTACCGTACCCAGGCGGGCCGTAACCTCCCCCCCAtt                                                                                                                                          >  1:1174545/1‑151 (MQ=255)
           aTTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGGTTGCCACTCCCAGTACAGGCGGGCTACGCGCCCAACCAGCAATTGGCGGGGTTTCTTGCCGTTCGCCGCCCCTCCCTTTAACGTAACCAGCCCGGCCGGAAACTcccgcc                                                                                                                                              >  2:55434/1‑141 (MQ=255)
             ttGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCtt                                                                                                                                       <  2:1748257/150‑1 (MQ=255)
             ttGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTTCGCCGCCCCGGCTTTTACCCTACCCCGGCGGGGAGTAACCCCCCCCCCATTCttt                                                                                                                                      >  2:1219182/1‑151 (MQ=255)
               gctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGGTCGCGGTCCGCCGCCCCTGCTTTAACCGTACCCCGGGGGGCAGGAACCCCCCCCCGGTTTTTTcc                                                                                                                                    >  1:1110392/1‑151 (MQ=255)
                   cAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTTCGCCGCCCGTGCTTTTAACCTACCCAGGCGGGCAGTAAACTCCCCCCGATTCTTTtccc                                                                                                                                  >  1:579478/1‑149 (MQ=255)
                    aaGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGa                                                                                                                               <  1:828674/151‑1 (MQ=255)
                    aaGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCGCCGCCCGTGGTTTAACCGTACCCAAGCGGGCAGTAACCCCCCCCCCGTTCTTTtccc                                                                                                                                  >  1:1037228/1‑148 (MQ=255)
                     aGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGat                                                                                                                              <  1:593464/151‑1 (MQ=255)
                      gACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGata                                                                                                                             <  1:1703866/151‑1 (MQ=255)
                      gACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCGCCGCCCGTGCTTTTACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTTCCCAGa                                                                                                                               >  1:1417360/1‑149 (MQ=255)
                       aCCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGatat                                                                                                                            <  1:421437/151‑1 (MQ=255)
                        ccGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGGTCCGCCGCCCCGGCTTTTACCCTAACCCAGCGGGGCGGAAACCCCCCCCCGCTTTTTTccc                                                                                                                                   >  2:535143/1‑144 (MQ=255)
                         cGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCCCGGTCCCCCGCCCCGGCTTTTAACCCACCCCAGGGGGGAGGAAACTCCCCCCCCTTCTTTtcccccata                                                                                                                             >  2:111397/1‑143 (MQ=255)
                            gtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCCCGGTTCCCCGCCCCGGCTTTAACCGTACCCCCGGGGGGCGGAAACCCCCCCCCAATTTTTTccc                                                                                                                                   >  2:326432/1‑140 (MQ=255)
                            gtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTTCGCCCCCCGGGCTTTTACCCTACCCAGGCGGGCAGTAAACTCCCCCCGCTTCTTTTCCCCGata                                                                                                                             >  1:760218/1‑145 (MQ=255)
                             tttAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGGTTCGCCGCCCCCGCTTTTAACGTACCCCGGCGGGGAGGAACCTCCGCCCGGTTCTTTTCCCCGATATCgag                                                                                                                        >  2:96434/1‑149 (MQ=255)
                              ttAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGTTCCCCCCCCCGTGCTTTAAACGGACCCAGGCGGGCAGTAACCTCCCCCCCATTCTTTCCCCCGATATCgaga                                                                                                                       >  1:1009611/1‑149 (MQ=255)
                               tAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCCCGGTCCGCCGCCCGTGCTTTAACCGGACCCAGGCGGGCAGTAAACTCCCCCCCATTCTTTCCCCAGATATCGAGATGc                                                                                                                    >  1:378340/1‑151 (MQ=255)
                                  cgccgcCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCCCGCGTTCCCCCCCCCCTGCTTTTAACGTACCCCGGCGGGCAGTAAACCCCCCCCGATTCTTTCCCCCGATATTGAGATGcc                                                                                                                   >  2:534004/1‑149 (MQ=255)
                                  cgccgcCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGGTCCGCCCCCCCTGCTTTTAACCTACCCCAGCGGGGAGGAAACTCCCCCCGATTCTTTTCCCCGATAACGAGAAGcc                                                                                                                   >  1:408558/1‑149 (MQ=255)
                                      gccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCCCCCGATTCTTTTCCCCGATATCGAGATGCCAGccc                                                                                                              >  1:874329/1‑150 (MQ=255)
                                       ccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCCCGCTTCCCCGCCCGTGCTTTTACCGCACCCCGGCGGGCAGGAAACCCCCCCCGATTCTTTCCCCCAATATCGAGATTCCAGccc                                                                                                              >  1:312294/1‑149 (MQ=255)
                                        cTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCCCGCGTTCCGCCCCCCGGGGTTTTAACCTAACCCGGGGGGGCGTGACCCCCCCCCCCTTTTTTTCCCCCATATCCAGATGGCAGcccccc                                                                                                           >  2:1411038/1‑151 (MQ=255)
                                        cTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCCCGTTCCCCCCCCCGTGCTTTAACCCTACCCAGGCGGGCCGGAACCTCCCCCCCCTTCTTTCCCCCGATATCGAGATGCCCGcccccc                                                                                                           >  1:255067/1‑151 (MQ=255)
                                         tGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGGTTCCCCCCCCCTGCTTTTAACCTAACCCAGCGGGGCGGAAACCCCCCCCCATTCTTTTCCCCCAAATCCAGGTGcccccccccc                                                                                                           >  2:731110/1‑142 (MQ=255)
                                          gggTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTACCCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGt                                                                                                         <  1:774169/151‑1 (MQ=255)
                                          gggTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCccgccg                                                                                                          <  2:392142/150‑1 (MQ=255)
                                            gTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTCCGCCCCCCCTGCCTTTACCCTACCCCGGCGGGCAGTAAACCCCGCCCCCTTTTTTTCCCCGATATCGAGATGCCCACcccccgct                                                                                                        >  2:172762/1‑148 (MQ=255)
                                             tttGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTCCCCCGCCCCTGCTTTTAACGGACCCAGGGGGGGAGGAAACCCCCCCCCGTTTTTTCCCCCGATATAGAGATGGCAACccgcccct                                                                                                        >  2:64523/1‑146 (MQ=255)
                                                 ccccTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTACCCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGc                                                                                                    <  1:927249/146‑1 (MQ=255)
                                                  cACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCCGGCGGGCAGTAACCTCCCCCCGATTCTTTTCCCAGATATCGAGATGCCAGCCCGCCGTTAAGcc                                                                                                   >  1:692790/1‑149 (MQ=255)
                                                    cTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCCCGTTTCCCCGCCCCTGCTTTTAACCCACCCCGGCGGGCCGTAACCCCCCCCCCATTCTTTCCCCAGATATCGAGATGCCAGCCCCCCCTTAAGCCcaa                                                                                                >  2:1321981/1‑150 (MQ=255)
                                                       ccAGTACAGGCGGGTCCCGCTCCCAGCCACAAATGGGGGGTTTTGCTGGCTTCCCGCCCCCCGTGCTTAACCCGTACCCAGGCGGGAAGAACCCCCCGCCCGATTCTTTCCCCAGATACCGAGAGCCCACCCCGCCGTTAAGCCAAAAGTa                                                                                            <  2:408933/151‑1 (MQ=255)
                                                        cAGTACAGGCGGGCTACGCTGCCAGCCAGCATTTGGCGGTTTTGCTCGCGTTCCGCCGCCCGTGCTTTACCCGTCCCCAGGCGGGCAGTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTAc                                                                                           <  2:375179/151‑1 (MQ=255)
                                                         aGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTCCCCAGGCGGGCAGTACCCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTAc                                                                                           <  1:1507232/150‑1 (MQ=255)
                                                          gtaggggggggCGACGCTGCCAGCCAGCAATTGGGGGTTTTGCTCGCTTCCCGCCGCCCGTGCTTTAACCGTCCCCGGGGGGGCAATACCCCCCGCCCGATTCTTCCCCCAGATATCGAGTGGCCAGCCCGCCGTTAAGCCAAAAGTACCg                                                                                         <  2:968162/143‑1 (MQ=255)
                                                          gTACGGGGGGGCTCCGCGCCCAGCCAGCATTTGGCGGTTTTGCCCGCGTCCCCCCGCCCGTGCTTAACCCGTCCCCGGGGGGGCAGTACCCCCCCCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTAcc                                                                                          <  2:1234723/150‑1 (MQ=255)
                                                          gTACCGGCGGTCTACTGTGCGAGCGCGAATTTGGGGGTTTCGCGCGCGTCCGCCCGCCCGTGTTTTAGCGGGACCTGGCTGGTCACTGCCCGCCGCTCGATCCTTGCCCGAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTAc                                                                                           <  2:172109/149‑1 (MQ=255)
                                                          gTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGTTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCg                                                                                         <  1:1613468/151‑1 (MQ=255)
                                                            aCAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTCCGCCGCCCGTGCTTTAAACCTACCCAGGCGGGCAGTAACCTCCGCCCCATTTTTTCCCCCGATATCGAGATGCCAGCCCGCCCTTAAACCAAAAATAcc                                                                                          >  1:1204880/1‑148 (MQ=255)
                                                              aGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCCATTCTTTCCCCCGATATCGAGATGCCAGCCCGCCGTTAAGCCCAAAGTACCGtt                                                                                       >  2:685190/1‑149 (MQ=255)
                                                               ggggggCTACGCTGCCACCCAGCATTGGGCGGGTTGGCTCGCGTCCCGCCGCCCGTGTTTTACCCGTCCCCAGGCGGGCAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTgg                                                                                     <  2:1134041/147‑1 (MQ=255)
                                                               ggCGGGCTCCGCTCCCAGCCAGCAATTGGCGGTTTTGCCCGCTTCCCGCCCCCCTTCTTTTACCCTTCCCCAGGGGGGCAGTACCCCCCCCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTAcc                                                                                          <  2:703705/145‑1 (MQ=255)
                                                               ggCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTgg                                                                                     >  1:466829/1‑150 (MQ=255)
                                                                  ggtctccgcgGCCAGCCAGCAATGGGCGGTTTTGCTCGCTTCCCGCCCCCCGTGCTTAACCCGTCCCCAGGGGGGCAGAACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTa                                                                                 <  2:1271104/145‑1 (MQ=255)
                                                                  gggCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGTTTTAACCGTCCCCAGGCGGGCAGTACCCTCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTa                                                                                 <  1:508386/151‑1 (MQ=255)
                                                                  gggCTACGCGCCCAGCCAGCAATGGGCGGTTTTGCCCGCGTCCCGCCCCCGGTGCTTTACCCGTCCCCGGGGGGGCAGTACCCCCCGCCCGTTTCTTCCCCCAGATATCGAGATCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTa                                                                                 <  2:83641/151‑1 (MQ=255)
                                                                     cTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGGACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGGTAAGCCAAAAGTACCGTTGGGGTAc                                                                                >  2:211775/1‑149 (MQ=255)
                                                                      tACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGGTCCGCCGCCCGTTCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTTCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCACAAGGACCGTTGGTGGACCAc                                                                             >  2:1384428/1‑151 (MQ=255)
                                                                       ccGCTGCCACCCAGCAATGGGCGGGTTTCCTCGCGTTCCGCCCCCCGTGCTTTACCCGTCCCCGGGGGGGCAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACg                                                                            <  2:1738873/150‑1 (MQ=255)
                                                                           tGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGcc                                                                        >  1:699715/1‑151 (MQ=255)
                                                                             ccagccagCAATTGGCGGTTTTGCTCGCGTCCCGCCGCCCGTGCTTTACCCGTACCCAGCCGGGCAGTAACCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGcc                                                                        <  1:806108/149‑1 (MQ=255)
                                                                             ccagccagCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCCGGCGGGGAGTAACCTCCCCCCGATTCTTTCCCCAGAGATCGAGATTCCAGCCCGCCGTTAAGCCAAAAGGACCGTTTGTGGACCACGGGCCgg                                                                      >  2:356373/1‑151 (MQ=255)
                                                                             ccagccagCAATTGGCGGGTTTGCCCGCGTCCCGCCCCCCTTGCTTTACCCGTACCCAGGCGGGCAGTACCCCCCCCCCGATTCTTTCCCCAGATATCGAGATGCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGcc                                                                        <  2:1692013/149‑1 (MQ=255)
                                                                              cagccagcAATTGGCGGGTTTATTCCCCTTTCTGCCGCCGCGCTTTTACCCGACCCAGCCGGGCGTTAACCTCTGACCCCTTCTTTACCCCGGTGTCGAGATGCCCGCCCGGCAGTTCGCCAAAAGTaatgtaggggtaaca                                                                              >  2:1477657/1‑128 (MQ=255)
                                                                               agccagcAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCgg                                                                      >  1:431706/1‑149 (MQ=255)
                                                                                 ccagcAATTGGCGGGTTTGCTCGCGTCCCGCCGCCCGTGCTTTACCCGTACCCAGGCGGGCAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGt                                                                  <  1:1622762/151‑1 (MQ=255)
                                                                                  cagcAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGt                                                                  >  1:225337/1‑150 (MQ=255)
                                                                                    gcAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCCCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCCCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTAcc                                                               >  1:954221/1‑151 (MQ=255)
                                                                                     cAATTGGCGGGTTTGCTCGCGTTCCGCCCCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAATACCGTTGGTGTACCACGGGCCGGTCGGAcc                                                               >  2:1726401/1‑150 (MQ=255)
                                                                                      atttGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTACCCGTCCCCAGGCGGGCAGTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTg                                                             <  1:211781/149‑1 (MQ=255)
                                                                                      aaTTGGGGGTTTTGCTCGCGTTCCGCCCCCCGTGCTTTAACCGTCCCCAGGCGGGCAGTACCCCCCGCCCGATTCTTCCCCAAGATATCGAGATCCCGCCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTc                                                                    <  2:149505/144‑1 (MQ=255)
                                                                                      aaTTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTACCCGTACCCAGGCGGGCAGTACCCTCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTg                                                             <  1:702737/151‑1 (MQ=255)
                                                                                      aaTTGGCGGGTTTGCTCGCGTCCCGCCGCCCGTGCTTTACCCGTCCCCAGGCGGGCAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTg                                                             <  2:1417247/151‑1 (MQ=255)
                                                                                        ttGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCg                                                           >  2:1523255/1‑151 (MQ=255)
                                                                                                          gTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGc                                          >  2:561294/1‑151 (MQ=255)
                                                                                                           tcccgccgCCCGTGTTTTAACCGTCCCCGGGCGGGCAGTACCCCCCCCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGc                                          <  1:1327547/148‑1 (MQ=255)
                                                                                                           tcccgccgCCCGTGCTTTAACCGCCCCCAGGCGGGCAGTACCCCCCCCCCGATTCTTCCCCCAGATAGCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGc                                          <  2:1012237/148‑1 (MQ=255)
                                                                                                           tcccgcccCCCGTGCTTTACCCGTCCCCAGGGGGGAAGAACCCCCCCCCCGATTTTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCaga                                            <  2:135913/146‑1 (MQ=255)
                                                                                                            tCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCa                                         >  1:731110/1‑150 (MQ=255)
                                                                                                             ccgccgCCCGTCATTAACCCGTACCCAGGCGGGCAGAACCCCCCCCCCAATTCTTCCCCCAGATATCGATAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCa                                         <  2:1760417/149‑1 (MQ=255)
                                                                                                              cgccgcCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAAc                                      >  2:304634/1‑151 (MQ=255)
                                                                                                                ccccccGTGCTTTACCCGTCCCCGGGGGGGCAGTACCCCCCGCCCGTTTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACg                                     <  2:1605843/147‑1 (MQ=255)
                                                                                                                        tGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGCTTCTTTCCCCCGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCAtt                             >  1:872518/1‑151 (MQ=255)
                                                                                                                            tttACCCGTACCCAGGCGGGCAGTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTaa                           <  1:809651/150‑1 (MQ=255)
                                                                                                                            tttAACCGTCCCCAGGCGGGAAGTACCCTCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTa                            <  1:1021385/149‑1 (MQ=255)
                                                                                                                             tttcccGTCCCCAGGCGGGCAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAg                          <  1:1166234/146‑1 (MQ=255)
                                                                                                                              tAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGcc                        >  1:106294/1‑151 (MQ=255)
                                                                                                                               acccGTACCCGGGCGGGCAGTAACCCCCGCCCGATTCTTCCCCCAGAAAGCGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGcc                        <  1:491621/148‑1 (MQ=255)
                                                                                                                                 ccGTCCCCAGGCGGGCAGTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTc                     <  1:6177/151‑1 (MQ=255)
                                                                                                                                 ccGTACCCAGGCGGGCAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTc                     <  1:957567/151‑1 (MQ=255)
                                                                                                                                 ccGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTc                     >  2:751904/1‑151 (MQ=255)
                                                                                                                                   gTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTg                   >  1:836264/1‑151 (MQ=255)
                                                                                                                                   gTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGAGCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTg                   >  2:1343121/1‑151 (MQ=255)
                                                                                                                                    tACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGc                  >  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)
                                                                                                                                                          |                                                                                                                                                  
GCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGC  >  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: