Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,228,2750TC50.7% 53.2 / 82.4 67pseudogene (335/1287 nt)mdtQpseudogene; putative channel/filament proteins
Reads supporting (aligned to +/- strand):  ref base T (16/17);  new base C (10/24);  total (26/41)
Fisher's exact test for biased strand distribution p-value = 1.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.04e-07

TCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCGTTCCG‑CCGCCCGTG‑C‑TTTAACCGTACCC‑AGGCGGGCA‑GTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTT  >  NC_000913/2228128‑2228415
                                                                                                                                                         |                                                                                                                                            
tctTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTTTTCGCGTTTCGCCCCCCCGGC‑C‑TTTAACCGCCCCC‑CGG‑GGGCC‑GTAACCTcc                                                                                                                                            >  2:480434/1‑150 (MQ=255)
 ctTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTGCTCGCGGTTCC‑CCCCCCCCG‑G‑TTTAAACGTACCC‑CGGCGGGGA‑GGAACcccc                                                                                                                                            >  1:1085478/1‑149 (MQ=255)
   ttttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTGCTCCCGGTCCG‑CCGCCCCCG‑C‑TTTTACCCCACCC‑CGGCGGGGA‑GGAACcccccccc                                                                                                                                        >  2:406211/1‑151 (MQ=255)
     ttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCGTTCCG‑CCGCCCGTG‑C‑TTTAACCGTACCC‑AGGCGGGCA‑GTAACCTCCGCCCg                                                                                                                                       <  1:607729/150‑1 (MQ=255)
     ttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCGTTCCG‑CCGCCCGTG‑C‑TTTAACCGTACCC‑AGGCGGGCA‑GTAACCTCCGCCCg                                                                                                                                       <  2:1278638/150‑1 (MQ=255)
     ttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTTCTCGCGTTCCG‑CCCCCCGGG‑CTTTTAACCG‑ACCC‑CGGCGGGCA‑GGAACCTcccccc                                                                                                                                        >  1:685465/1‑149 (MQ=255)
      tCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTTCTCGCGGTTCG‑CCCCCCGCG‑CTTTTAACCCAACCC‑AGGCGGGGG‑GAAACCCCCCCCCgt                                                                                                                                      >  2:1242831/1‑150 (MQ=255)
         aTTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCGTTCCG‑CCGCCCGTG‑C‑TTTAACCGTACCC‑AGGCGGGCA‑GTAACCTCCGCCCGAtt                                                                                                                                    <  2:630909/149‑1 (MQ=255)
         aTTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCGTTCCG‑CCGCCCGTG‑C‑TTTAACCGTACCC‑AGGCGGGCA‑GTAACCTCCGCCCGAt                                                                                                                                     <  2:404564/148‑1 (MQ=255)
              ctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCCGCAATTGGCGGGGTTTGC‑CGCGGTTCC‑CCGCCCCCGCC‑TTTAACC‑TACCC‑CAGCGGGGAGGAAACCCCCCCCCGCTTCTTccc                                                                                                                              >  2:140088/1‑151 (MQ=255)
                gcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCGTTCCG‑CCGCCCGTG‑C‑TTTAACCGTACCC‑AGGCGGGCA‑GTAACCTCCGCCCGATTCTTTcccc                                                                                                                            <  1:969332/150‑1 (MQ=255)
                  aaGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTAAGCTGCCAGCCAGCAATTGGC‑GGGCTTTTGCTCGTTTCG‑CCGCCCCCG‑C‑CTTTACCCGACCC‑CCGGGGGCA‑GTTACCCCAGCCCGATTCtttttcc                                                                                                                            >  2:628036/1‑144 (MQ=255)
                       cGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCGTTCCG‑CCGCCCGTG‑C‑TTTAACCGTACCC‑AGGCGGGCA‑GTAACCTCCGCCCGATTCTTTCCCCAGATATCg                                                                                                                    <  2:226844/151‑1 (MQ=255)
                       cGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTTTTCCCGTTTCG‑CCCCCCCTG‑C‑TTTTACCCTACCC‑CGGGGGGGC‑GTAACCCCCCCCCCCTTTTTTTCCCACatat                                                                                                                      >  2:916588/1‑149 (MQ=255)
                        gtgtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCGTTCCG‑CCGCCCGTG‑C‑TTTAACCGTACCC‑AGGCGGGCA‑GTAACCTCCGCCCGATTCTTTCCCCAGATATCga                                                                                                                   <  2:374222/151‑1 (MQ=255)
                              aaCGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCGTTCCG‑CCGCCCGTG‑C‑TTTAACCGTACCC‑AGGCGGGCA‑GTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGc                                                                                                              <  1:1210028/150‑1 (MQ=255)
                                cgccgcCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCGTTCCG‑CCGCCCGTG‑C‑TTTAACCGTACCC‑AGGCGGGCA‑GTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAg                                                                                                           <  1:691275/151‑1 (MQ=255)
                                cgccgcCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTTCCCGCGTTCCC‑CCGCCCGCG‑G‑CTTTACCGCACCC‑CCGCGCGGC‑GTAAACCCCGCCCCCTTTTTTCCCCCGATATCGGGATGcc                                                                                                             >  2:135789/1‑149 (MQ=255)
                                  ccgccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTTCTCGCGGTTCG‑CCGCCCCTG‑C‑TTTTACCGGACCC‑CGGCGGGGA‑GGAAACCCCCCCCGCTTTTTTTCCCCGATTTCGAGATGgccgcc                                                                                                         >  2:405183/1‑151 (MQ=255)
                                   cgccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCGTTCCG‑CCGCCCGTG‑C‑TTTAACCGTACCC‑AGGCGGGCA‑GTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGccc                                                                                                        <  1:971808/151‑1 (MQ=255)
                                   cgccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTGCTCGGGTTCCG‑CCGCCCGTG‑C‑TTTATCCGTACCC‑AGGCGGGGA‑GTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGccc                                                                                                        <  1:1102931/151‑1 (MQ=255)
                                      cTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTTCTCGCGGTCCG‑CCGCCCGGG‑C‑TTTTACCGTACCC‑AGGCGGGGA‑GTAACCCCCCCCCGATTCTTTTCCCCGATATCGAGATGCCAGcccccc                                                                                                     >  1:952141/1‑151 (MQ=255)
                                             tGCCACTCCCAGTACAGGCGGGGTACGCTGCCAGCCAGCAATTTGC‑CGGTTTGCTCCCGCTCCC‑CCGCCCCGG‑C‑GTTTACACTGCCC‑CGGGGGGCA‑GGAACCCGCCCCCCCTTCTTTCTCCCCATATCGCGGTGCCACCccgccccgaag                                                                                               >  2:457709/1‑144 (MQ=255)
                                                           cAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTTCTCGCGTTCCG‑CCGCCCGTG‑C‑CTTAACCCTACCC‑CGGCGGGCC‑GGAACCCCCCCCCCGCTCTTTTCCCAGATATTGCGATGGCACGCCCCCCTTAAGGCCAAAggaac                                                                                    >  2:426519/1‑143 (MQ=255)
                                                           cAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCGTTCCG‑CCGCCCGTG‑C‑TTTAAACGTACCC‑AGGCGGGCA‑GTAACCTCCCCCCGATTCTTTCCCCAGATATCGAGATGGCAGCCCGCCGTTAAGCCCAAAGAACCGtt                                                                                 >  2:653400/1‑150 (MQ=255)
                                                              gCGGGCTCCGCTGCCGGCCAGCAATTGGC‑GGGTTTGCTCGCGTCCCG‑CCCCCCGTT‑C‑TTTACCCGTCCCC‑AGGCGGGCA‑GTACCCTCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGgtg                                                                             <  2:11938/151‑1 (MQ=255)
                                                                  gCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCGTTCCG‑CCGCCCGTG‑C‑TTTAACCGTACCC‑AGGCGGGCA‑GTAACCTCCCCCCGCTTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGCTGGTGTAcc                                                                         >  2:200727/1‑151 (MQ=255)
                                                                    tACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCGTTCCG‑CCGCCCGTG‑C‑TTTAACCGTACCC‑AGGCGGGCA‑GTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGAACCGTTGGTGTAc                                                                          >  1:948955/1‑148 (MQ=255)
                                                                      ccctccccgccaccAATTGGC‑GGGTTTGCTCGCTTCCCC‑CCGCCCGTG‑C‑TTTAACCGTCCCC‑AGGGGGGCA‑GTCACCCCCGCCCGATTTTTCCCCCATATCTCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTAcc                                                                         <  2:170177/139‑1 (MQ=255)
                                                                      cGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCGGTCCG‑CCGCCCGTG‑C‑TTTAAACGTACCC‑AGGCGGGCA‑GTAACCTCCGCCCGATTCTTTCCCCCGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGGACCGTTGGTGTACCACgg                                                                     >  1:58106/1‑151 (MQ=255)
                                                                      cGCTGCCACCCAGCAATGGGG‑GGGTTTGCTCGCTTCCCC‑CCGCCCGTT‑T‑TTTACCCGTCCCC‑AGGCGGGCA‑GTACCCCCCGCCCGTTTCTTCCCCCAGATTTCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCa                                                                        <  2:808163/148‑1 (MQ=255)
                                                                       gCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCGTTCCG‑CCGCCCGCG‑C‑TTTAACCGTACCC‑AGGCGGGCA‑GTAACCCCCCCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCCTTGGTGTACCACggg                                                                    >  2:934896/1‑151 (MQ=255)
                                                                        cTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCGTTCCC‑CCGCCCGTG‑C‑TTTAACCCTAACC‑AGGCGGGCA‑GTAAACTCCGCCCGAATCTTTTCCCCGATATCCAGATGCCAGCCCGCCCTTAAGCCCAAAGTACCGTTGGTGTACCACcgg                                                                    >  2:782351/1‑150 (MQ=255)
                                                                           ccagccagCAATTGGC‑GGGTTTGCTCGCGTTCCG‑CCGCCCGTG‑C‑TTTACCCGTCCCC‑GGGCGGGCA‑GTACCCCCCGCCCTATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGcc                                                                  <  1:344561/149‑1 (MQ=255)
                                                                            cggccgacAATTGGC‑GGTTTTCTTCGCGTTCC‑‑CCGGCCATGTC‑TTTAACCTTCCCC‑TGGAGTGCG‑ATTACTCCGGCTCGATTATCGCGCAATTTCTCGAGACGCAGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCg                                                                 <  2:911859/142‑1 (MQ=255)
                                                                             agcccccAATTGGC‑GGGTTTGCTCGCGTCCCG‑CCGCCGGTG‑G‑TTTAACGGCACCC‑CGCCGGTAG‑GTACCACCCGCCGAATTCTTGCCCCAAATACGGAGATGCCCGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGc                                                                   <  1:311383/140‑1 (MQ=255)
                                                                               caagcAATTGGC‑GGTTTTGCGCGCTTTCCC‑CCGCCCGTG‑C‑TTAACCCGTACCC‑GGGCGGGTA‑TTACCCCCCGCCCGTTTTTTTCCCCAGATATCGAGATCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGc                                                                   <  2:937582/142‑1 (MQ=255)
                                                                                      tgggC‑GGTTTTGCTCGCGTTCCG‑CCGCCCGTG‑T‑TTTACCCGTCCCC‑AGGCGGGCA‑GTAACCCCCCCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGc                                                      <  1:128426/148‑1 (MQ=255)
                                                                                       tGGC‑GGGTTTGCTCGCGTTCCG‑CCGCCCGTG‑C‑TTTAACCGTACCC‑AGGCGGGCA‑GTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCgg                                                    >  1:995675/1‑151 (MQ=255)
                                                                                        ggC‑GGGTTGCCTCGGGTTCCC‑CGGCCGGTC‑C‑TTTACCCTTACCC‑AGGGGGGCG‑GTACCCCCCGCCGGTTTTTTCCCCCGGATACCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCg                                                     <  1:492526/149‑1 (MQ=255)
                                                                                         gggggttttGCTCGCGTCCCG‑CCGCCCGTG‑C‑TTTACCCGTCCCC‑GGGCGGGCA‑GTACCCCCCGCCCGTTTCTTCCCCCAGATATCGAGATGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCgg                                                    <  2:165186/144‑1 (MQ=255)
                                                                                         gC‑GGGTTCGCCCGCGTCCCC‑CCGCCCGTG‑T‑TTTACCCGTCCCA‑GGGGGGGCA‑GTACCCCCCGCCCGTTTCTTTCCCCAGATAACGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGcc                                                  <  2:989269/151‑1 (MQ=255)
                                                                                          c‑GGGTTTGCTCGCGTTCCG‑CCGCCCGTG‑C‑TTTAACCGTACCC‑AGGCGGGCA‑GTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCg                                                 >  2:711665/1‑151 (MQ=255)
                                                                                            gggTTTGCCCGGTTCCCG‑CCGCCCGTG‑C‑TTTACCCGTCCCC‑GGGCGGGCA‑GAACCCTCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCgg                                                <  1:197255/151‑1 (MQ=255)
                                                                                             gttttGCTCGCGTCCCC‑CCGCCCGTG‑T‑TTTACCCGTCCCA‑GGGCGGACA‑GAACCCCCCGCCCTTTTTTTCCCCCAGATATCGGGGTCCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCgg                                                <  2:191536/148‑1 (MQ=255)
                                                                                               tttGCTCGCTTCCCG‑CCCCCCGTG‑T‑TTAACCCGTACCC‑GGGCGGGAA‑GTAACCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGAt                                              <  1:753510/150‑1 (MQ=255)
                                                                                               tttGCTCGCGTTCCG‑CCGCCCGTG‑C‑TTTATCCGTACCC‑AGGCGGGCA‑GTAACCTCCGCCCGATTCTTTCCCCAGATATCGCGATGCCAGCCCGCCGTTATCCCAAAAGTACCGATGACGTCCCTCGCGCCGGTCGTACCTGCGGCCGGa                                               >  2:975839/1‑149 (MQ=255)
                                                                                                   cccgcgTCCCC‑CCGCCTTTC‑T‑TTTAACCGTCCCA‑GGGCGGGCA‑GCACCCTCCCCCCGATTCTTGCCCAAGATATGGATATCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATc                                             <  2:790102/145‑1 (MQ=255)
                                                                                                      gctttCCC‑CCCCCTGTG‑T‑TTTACCCGTCCCC‑GGGCGGGTA‑GTACCCCCCGCCCGATTTTTCCCCCAGATATGGAGATCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGt                                           <  1:369751/143‑1 (MQ=255)
                                                                                                      gcgTTCCG‑CCG‑CCGGG‑G‑TGTAACCCTTCCCGCGGCGGGAG‑GTACCCCCCGCCGGATTTTTTCCCCAGATAACGAGACGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGtt                                          <  1:1248511/147‑1 (MQ=255)
                                                                                                        gTTCCG‑CCGCCCGTG‑C‑TTTAACCGTACCC‑AGGCGGGCA‑GTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGc                                    >  2:1267427/1‑151 (MQ=255)
                                                                                                           ccg‑ccgCCCGTG‑T‑TTTAACCGTCCCC‑AGGCGGGCA‑TTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGc                                    <  2:1251360/148‑1 (MQ=255)
                                                                                                           ccg‑ccgCCCGTG‑C‑TTTAACCGTCCCC‑GGGCGGGCA‑GTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCa                                   <  2:203447/149‑1 (MQ=255)
                                                                                                           ccg‑ccgCCCGTG‑C‑TTTAACCGTACCC‑AGGCGGGCA‑GTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCaa                                  >  1:344259/1‑150 (MQ=255)
                                                                                                               ccgccCGTG‑C‑TTTAACCGTACCC‑AGGCGGGCA‑GTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACgg                              >  2:960347/1‑151 (MQ=255)
                                                                                                                  cccGTG‑T‑TTTACCCGTCCCC‑GGGCGGGCA‑GTACCCCCCCCCCGTTTTTTCCCCCGGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAAc                                <  2:683883/146‑1 (MQ=255)
                                                                                                                  cccGTC‑C‑TTTAGCCGTCCCG‑GGGCGGAAA‑GTGACCCCCGCCCGTTTTTTCCCCCAAATGTCGAAAGCCCAGGCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATc                                             <  1:1165899/133‑1 (MQ=255)
                                                                                                                      tG‑C‑TTAAACCGTACCG‑GGGCGGGGA‑GTAACCTCCGCCCGATTTTTCCCCCAAATACAAAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCTCCAt                        <  1:879550/150‑1 (MQ=255)
                                                                                                                            ttAACCGTCCCG‑AGGAGGGCA‑TTACCCCCCGCCAGATTTTTTCCCAAGATATCGAGGTGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAg                    <  1:1226809/150‑1 (MQ=255)
                                                                                                                              accctaaCCC‑GGGGGGGGA‑TTACCCCCCGCCCGATTCTTTCCCAAAATAGCGAGATGCCGCCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTa                      <  1:295523/140‑1 (MQ=255)
                                                                                                                              acccGTCCCC‑AGGGGGGCA‑GTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGcc                  <  1:1257026/148‑1 (MQ=255)
                                                                                                                                cggtcccc‑AGGGGGGAA‑GTACCCCCCGCCCGTTTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCtt                <  1:533625/145‑1 (MQ=255)
                                                                                                                                cctcaCCC‑GGGCGGGCG‑ATACCCCCCCCCCGATTCTTCCCCCAGATAGCGGGAGGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCtt                <  2:304077/146‑1 (MQ=255)
                                                                                                                                ccGTCCCC‑AGGCGGGCA‑GAACCCTCCCCCCGTTTCTTTCCCCAGATACCGAGATGCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGcc                  <  1:480454/148‑1 (MQ=255)
                                                                                                                                 cGTACCC‑AGGCGGGCA‑GTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTc               >  1:579096/1‑150 (MQ=255)
                                                                                                                                  gTACCC‑AGGCGGGCA‑GTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTg             >  1:829407/1‑151 (MQ=255)
                                                                                                                                  gTACCC‑AGGCGGGCA‑GTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTg             >  2:977701/1‑151 (MQ=255)
                                                                                                                                             gggCA‑TTACCCCCCCCCCGTTTTTTTCCCCGGATATGGAGATGCCCGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAttt    <  1:1012241/150‑1 (MQ=255)
                                                                                                                                              ggta‑GTACCCCCCGCCCGTTTCTTTCCCGACATATGGGAAGCCCATCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAttttt  <  1:628060/148‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                            
TCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCGTTCCG‑CCGCCCGTG‑C‑TTTAACCGTACCC‑AGGCGGGCA‑GTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTT  >  NC_000913/2228128‑2228415

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