Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,228,2220TG41.8% 87.1 / 10.9 55pseudogene (388/1287 nt)mdtQpseudogene; putative channel/filament proteins
Reads supporting (aligned to +/- strand):  ref base T (9/23);  new base G (23/0);  total (32/23)
Fisher's exact test for biased strand distribution p-value = 1.98e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGG  >  NC_000913/2228074‑2228369
                                                                                                                                                    |                                                                                                                                                       
gaagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGttt                                                                                                                                                       <  1:224662/151‑1 (MQ=255)
     aaTCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTTGCGGGt                                                                                                                                                         >  1:797433/1‑144 (MQ=255)
     aaTCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTcg                                                                                                                                                  <  1:882950/151‑1 (MQ=255)
             tCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGtttc                                                                                                                                            >  1:928189/1‑147 (MQ=255)
                 atatAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCgccgcc                                                                                                                                      <  1:199298/151‑1 (MQ=255)
                    taGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTCCGCCGccc                                                                                                                                     >  1:608992/1‑149 (MQ=255)
                      gCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCCCGTTTCGCCGcccct                                                                                                                                   >  1:343924/1‑147 (MQ=255)
                         ggCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCttt                                                                                                                             <  2:355348/151‑1 (MQ=255)
                           cGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTCCGCCCCCCGTGCttt                                                                                                                             >  1:98053/1‑149 (MQ=255)
                             aTCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCCCGTTCCGCCCCCCCTGCttt                                                                                                                             >  1:78513/1‑147 (MQ=255)
                                   cgcgATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGGTTCCCCGCCCGTGCTTTTAACGTAccc                                                                                                                  >  2:49110/1‑151 (MQ=255)
                                     cgATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGGTCCCCCCCCCGTGCTTTTACCCGAcccc                                                                                                                 >  1:556667/1‑149 (MQ=255)
                                      gATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAg                                                                                                                <  1:803782/150‑1 (MQ=255)
                                      gATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCCGCCCGCCAGCAATTGGCGGGTGTTGTCGCGGTCCTCCCCCCCTGCTTTTACCCGACCCAgg                                                                                                               >  2:403823/1‑151 (MQ=255)
                                         aaTGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGGTCCCGGTCCGCCGCCCGCGCTTTTAACCCACCCcggcgg                                                                                                           >  2:934206/1‑151 (MQ=255)
                                              tgtTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTGTCGCGGTTCGCCCCCCGGGCTTTTAACGTACCCAGGCggg                                                                                                          >  2:319098/1‑147 (MQ=255)
                                              tgtTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTTCCGTTTCCCCGCCCCGGCTTTTAACCCACCCACGCggg                                                                                                          >  1:579520/1‑147 (MQ=255)
                                               gtTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGt                                                                                                      <  1:322748/150‑1 (MQ=255)
                                                 tCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCCCGTTCCCCCGCCCGGGCTTTTACCGGACCCAGGCGGGCAGTaa                                                                                                    >  1:235076/1‑150 (MQ=255)
                                                    gCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAAc                                                                                                  <  1:486084/148‑1 (MQ=255)
                                                    gCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCt                                                                                                <  2:312238/150‑1 (MQ=255)
                                                     ctctTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTCCCCCGCCCGTGCTTTAACCGTACCCCGGCGGGCAGTAACCTcc                                                                                              >  1:958438/1‑151 (MQ=255)
                                                      tctTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGGTCGCGTTCCCCCGCCCGCGCTTTAACCGGACCCCGGGGGGGAGTAACcccc                                                                                              >  1:880392/1‑150 (MQ=255)
                                                       ctTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGGCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCTAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTcc                                                                                              <  1:998932/149‑1 (MQ=255)
                                                       ctTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTcc                                                                                              <  1:669927/149‑1 (MQ=255)
                                                        tttttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTTACCGTACCCAGGCGGGCAGTAACCTCCGcc                                                                                           >  1:633967/1‑151 (MQ=255)
                                                        tttttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCTCGTTCCCCCCCCCCTGCTTTTACCGTACCCAGGCGGGGAGTAAACTCcccc                                                                                           >  2:290545/1‑151 (MQ=255)
                                                         ttttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCg                                                                                             <  2:979993/148‑1 (MQ=255)
                                                           ttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGGTCCGCCCCCCGGGGTTTAACCGGAACCCAGGGGGGCGTAACCTcccccc                                                                                          >  1:859085/1‑149 (MQ=255)
                                                            tCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGa                                                                                        <  1:983910/150‑1 (MQ=255)
                                                             cTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCCCGTTCCCCCGCCCGCGCCTTTACCCGCCCCCGGCGGGGCGGAAACTccccccccat                                                                                      >  2:640494/1‑147 (MQ=255)
                                                               aTTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGAt                                                                                       <  1:510619/148‑1 (MQ=255)
                                                               aTTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTTCGCCGCCCGTGCTTTAAACGTACCCCAGCGGGGAGGAAACTCCCCCCGATTCt                                                                                    >  1:108227/1‑151 (MQ=255)
                                                                  tgctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCCCCCCCCCGGCTTTTACCGGACCCCGGCGGGCAGTAAACTcccccccgttcctt                                                                                  >  1:261028/1‑142 (MQ=255)
                                                                      gcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTGTCGCGGTCCGGCGCCCGGGGTTTTAACCGTCCCCGGGGGGGAGGAAACCCCCCCCCGTTTTTTTcccc                                                                              >  2:780161/1‑151 (MQ=255)
                                                                       cAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAg                                                                            <  2:330409/151‑1 (MQ=255)
                                                                             cGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATc                                                                       <  1:951919/150‑1 (MQ=255)
                                                                               tgtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGGTTCGCCGCCCCTGCTTTAAACGTACCCAGGCGGGCAGTAAACTCCCCCCGATTCTTTTCCCAGata                                                                         >  1:837245/1‑146 (MQ=255)
                                                                                gtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCga                                                                     <  1:64841/149‑1 (MQ=255)
                                                                                 tttAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGGTCCGCCGCCCCCGCCTTTACCCGACCCAGGCGGGGCGGAAACCCCCCCCCGTTCTTTCCCCCGa                                                                           >  2:93430/1‑142 (MQ=255)
                                                                                   tAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCGAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCTCGTGCTTTAACCGTCCCCAGGCGGGCAGTACCCTCCGCCCGATTCTTTCCCCAGATATCTAGATg                                                                 <  2:518018/150‑1 (MQ=255)
                                                                                   tAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATg                                                                 <  1:668121/150‑1 (MQ=255)
                                                                                      cgccgcCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAg                                                             <  2:637400/151‑1 (MQ=255)
                                                                                         cgccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGccc                                                          <  2:295120/151‑1 (MQ=255)
                                                                                              gggTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTCCCCCGCCCGTGCTTTTACCGCAACCCGGCGGGCAGTAACCCCCCCCCCCATCTTTTCCCCGATATCGAGATGCCcgcccgcc                                                       >  2:840498/1‑149 (MQ=255)
                                                                                                  ttGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCGCCGCCCGGGCCTTAAACGGAACCAGGCGGGCAGGAAACTCCGCCCGATTTTTTTCCCAGATATCGAGATGCCCGCCCGCCGttta                                                  >  1:894380/1‑148 (MQ=255)
                                                                                                     ccACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTCCCCAGGCGGGCAGTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAg                                                 <  1:216904/148‑1 (MQ=255)
                                                                                                      cACTCCCATTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTCCCCAGGCGGGCAGTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCaa                                             <  1:797077/151‑1 (MQ=255)
                                                                                                      cACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCGCCGCCCGTGCTTTAAACGTACCCAGGCGGGCAGTAAACTCCGCCCGATTCTTTCCCCCGATAACGAGATGCCAACCCCCCGTTAAAcc                                               >  1:38403/1‑149 (MQ=255)
                                                                                                        cTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCCTTCCGCCGCCCCTGCTTTTAACCTACCCCGGCGGGGCGGAACCTCCCCCCCCATCTTTTCCCCGATATCGAGATGCCCGCCCGCCGGTTAACCaaa                                            >  2:438345/1‑150 (MQ=255)
                                                                                                            cAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGGTTCCCCCCCCCTGCTTTTAACCTACCCCGGGGGGCAGTAACCCCCCCCCCCCTTTTTTCCCCAATTTCCAGAGGCCAACCCCCCCCTAAACCcaaaaga                                        >  2:774858/1‑147 (MQ=255)
                                                                                                                aCAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTTCGCCGCCCGTGCTTTTACCCTACCCAGGCGGGGAGGAAACTCCCCCCCGTTCTTTTCCCCGATATCGAGATGGCAGCCCCCCCGTAAGCCaaaaagaccctt                                   >  1:632268/1‑143 (MQ=255)
                                                                                                                   gggggcttccGCTGCCAGCCAGCAATTGGGGGGTTTGCGCGCGTCCCCCCCCCCGTGTTTAACCCGTACCCAGGCGGGCGGTACCCTCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTg                                  <  2:678983/140‑1 (MQ=255)
                                                                                                                          tACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTTACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTAACGTTGGTGTAcccc                         >  2:40944/1‑149 (MQ=255)
                                                                                                                              cTGCGAGCCAGCAATCGGCGGATCTGCCCGCGTTCCGCGGCCCGTGATTAAACCGCACCCAGGGGGATACAAACCTCCGCCCGTTTTTTGCCACACATATCGAGAGCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACgg                       <  2:73216/149‑1 (MQ=255)
                                                                                                                                       agcAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTAc            >  1:892657/1‑151 (MQ=255)
                                                                                                                                                 ggttttGCTCGCGTCCCCCCGCCCGTGCTTTACCCGTCCCCGGGCGGGCAGAACCCCCCGCCCGATTCTTTCCCAAGATATCGAGATGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCg   <  2:337924/147‑1 (MQ=255)
                                                                                                                                                 gggTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCgg  >  1:59387/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
GAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGG  >  NC_000913/2228074‑2228369

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