Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,969,7940AC48.3% 85.1 / 14.5 61G397G (GGT→GGG) tapmethyl‑accepting protein IV
Reads supporting (aligned to +/- strand):  ref base A (27/4);  new base C (0/29);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 3.62e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCG‑CGCTACGGCTGGCAAGATTGCG‑CACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCT  >  NC_000913/1969647‑1969942
                                                                                                                                                      |                                                                                                                                                    
ggTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCG‑CGCTACGGCTGGCAAGATTGCG‑CACTTCACCt                                                                                                                                                   >  1:1088159/1‑151 (MQ=255)
     ccGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCG‑CGCTACGGCTGGCAAGATTGCG‑CACTTCACCTGcc                                                                                                                                                <  2:481308/149‑1 (MQ=255)
       gAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCG‑CGCTACGGCTGGCAAGATTGCG‑CACTTCACCTGCCACt                                                                                                                                             >  2:894168/1‑150 (MQ=255)
        aaCTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCG‑CGCTACGGCTGGCAAGATTGCG‑CACTTCACCTGCCACTAc                                                                                                                                           >  1:585344/1‑151 (MQ=255)
                        cATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCG‑CGCTACGGCTGGCAAGATTGCG‑CACTTCACCTGCCACTACCGCAAAACCACGcccc                                                                                                                           >  1:713586/1‑151 (MQ=255)
                          tGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCG‑CGCTACGGCTGGCAAGATTGCG‑CACTTCACCTGCCACTACCGCAAAACCACGCCCCTg                                                                                                                         >  2:360208/1‑151 (MQ=255)
                           ggcggcggcggcgTTATTCCCCATTTTCGACCCCGGCTGGACACGATTGACTGCCTCTTCGACCAGCCTTTTG‑ATCTCTTTTGCTGCCTGCG‑CGCTCGGGTTGGCAAGATTGCG‑CACTCCCCCTGCCACTACCGCAAAACCACGCCCCt                                                                                                                          <  2:826951/145‑1 (MQ=255)
                           ggTCGGGGCGGGGTTTTTCACCATTTCCGACCCCTGTGGGACACAATTGACTAACTCTTCGATCCCCCCTTTG‑ATCTTTTTTCCTCCCTGCC‑CGCTCCGGCTGGCAAGATTGCG‑CACTCCCCCTGCCACTACCGCAAAACCACGCCCCt                                                                                                                          <  2:605678/149‑1 (MQ=255)
                               gcggcggcgTTATTCCCCAGTTCCGACCCCTGTTGGACCCGATTGCCTGACTCTTCGACCACCCCTTTG‑ATCTTTTTGGCTGCCCGCG‑CGCTACGCCGGGCAATTTTGCG‑CCCTTCGCCTGCCACTACCGCAAAACCACGCCCCTGctct                                                                                                                     <  1:485061/150‑1 (MQ=255)
                               gcggcggcgTTATTCACCAGTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTTTTTTGCTGCCTGCG‑CGCTACGGCTGGCAAGATTGCG‑CACTTCACCTGCCACTACCGCAAAACCACGCCCCTGctct                                                                                                                     <  2:549180/150‑1 (MQ=255)
                                    ggggTTATTCCCCATTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTTTTTTGCTGCCTGCG‑CGCTCCGGCTGGCAAGATTGCG‑CACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGgc                                                                                                                <  2:788742/147‑1 (MQ=255)
                                       gTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCG‑CGCTACGGCTGGCAAGATTGCG‑CACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGcgcg                                                                                                             >  2:459370/1‑150 (MQ=255)
                                        ttATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCG‑CGCTACGGCTGGCAAGATTGCG‑CACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGcgcg                                                                                                             >  2:562140/1‑149 (MQ=255)
                                           ttCACCAGTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCTTTTG‑ATCTCTTTGGCTGCCTGCG‑CGCTACGGCTGGAAAGATTGCG‑CACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGcc                                                                                                         <  2:830005/150‑1 (MQ=255)
                                            tCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCG‑CGCTACGGCTGGCAAGATTGCG‑CACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGc                                                                                                       >  1:559210/1‑151 (MQ=255)
                                             cACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCG‑CGCTACGGCTGGCAAGATTGCG‑CACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCt                                                                                                      >  1:1182016/1‑151 (MQ=255)
                                                cAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCG‑CGCTACGGCTGGCAAGATTGCG‑CACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTcc                                                                                                   >  2:618384/1‑151 (MQ=255)
                                                   tttCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCG‑CGCTACGGCTGGCAAGATTGCG‑CACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACt                                                                                                >  1:1124992/1‑151 (MQ=255)
                                                    tccaaccccTGCTGGCCACTTTTGCCTGCCTCTTCGACCCCCCTTTTG‑ATCTTTTTTGCTGCCTGCG‑CGCTCCGGCGGGAAAGTTTGCC‑CCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACt                                                                                                <  1:28476/144‑1 (MQ=255)
                                                        accccTGCTGGACACGATTGACTGCCTCTTCGATCACCCTTTTG‑TTCTTTTTTGCTGCCTGCG‑CGCTCGGGCGGGCAAGATGGCG‑CATTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTg                                                                                               <  2:874652/145‑1 (MQ=255)
                                                           ccTGCTGGCCACGATTGACTGACTCTTCGATCCGCCCTTTG‑CCCTTTTTGGCTCCCCGCC‑CGCTGGGGCGGGCAAGATCGCC‑CCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCg                                                                                            <  1:66502/147‑1 (MQ=255)
                                                              gctgGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCG‑CGCTACGGCTGGCAAGATTGCG‑CACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTc                                                                                       >  2:1198011/1‑149 (MQ=255)
                                                                   acacGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCG‑CGCTACGGCTGGCAAGATTGCG‑CACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGcc                                                                                 >  1:424456/1‑150 (MQ=255)
                                                                    cacGTTTGACTGCCTCTCCGATCCCCCCTTTG‑ATCTTTTTGGCTGCCGGCG‑CGCTACGGCGGGAAAGTTGGCG‑CCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAg                                                                                     <  1:929262/145‑1 (MQ=255)
                                                                      cGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCG‑CGCTACGGCTGGCAAGATTGCG‑CACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGa                                                                              >  2:389298/1‑150 (MQ=255)
                                                                        aTTGACTGACTCTGCGTCCAGCCTTTTA‑ACCTCTTTTGCGGCCGCCG‑CGCTCGGGCGGGAAAGATCGCC‑CCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAg                                                                               <  2:756706/147‑1 (MQ=255)
                                                                                  ttctCGATC‑GCCCTTTGAATCTTTTTGCCTGCCGGGG‑CGCTGCGGGTGGAAAGTTTGCC‑CTCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTg                                                                        <  2:528097/141‑1 (MQ=18)
                                                                                    tccaacccgCCTTTTG‑ATCTC‑TTTGCCGCCCGCGCCACTACGGCGGGCAGGTTGCCG‑CCCGCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCt                                                                    <  1:52477/138‑1 (MQ=16)
                                                                                      cgctCACCCTTTTC‑ATCTTTTTTGCTCCCGGCG‑CCCGCGGGCTGGAAAGATTCCC‑CCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTg                                                                   <  1:274234/142‑1 (MQ=18)
                                                                                        aTCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCG‑CGCTACGGCTGGCAAGCTTGCG‑CACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAg                                                              >  1:991211/1‑148 (MQ=255)
                                                                                               ctttggtcctctTTTGCCCCCTCCG‑CCCTACGGC‑CGCAAATTTGCGCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAg                                                              <  1:869155/134‑1 (MQ=18)
                                                                                                tttG‑ATCTTTTTTGCTCCCCGCG‑CGCTCCGGCTGCCAAGATTGCG‑CCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGt                                                     <  2:218579/149‑1 (MQ=25)
                                                                                                tttG‑ATCTCTTTTGCTGCCTGCG‑CGCTACGGCTGGCAAGATTGCG‑CACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCg                                                   >  2:896199/1‑151 (MQ=255)
                                                                                                 tttgcttttttttGTTGCCGCCG‑CCCTGCGGCGGGAAGGTTTGCC‑CTCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTc                                                    <  2:429743/141‑1 (MQ=21)
                                                                                                 ttG‑ATCTCTTTTGCTGCCTGCG‑CGCTACGGCTGGCAAGATTGCG‑CACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTc                                                    >  2:1059137/1‑149 (MQ=255)
                                                                                                  tG‑ATCTTTTTTGCTCCCGGCG‑CGCTCGGGCTGGCAAGATTGCG‑CATTCCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGa                                                  <  2:1088195/150‑1 (MQ=37)
                                                                                                   g‑ATCTCTTTTGCTGCCTGCG‑CGCTACGGCTGGCAAGATTGCG‑CACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGAt                                                 >  2:222477/1‑150 (MQ=255)
                                                                                                        tctTTTGCTGCCTGCG‑CGCTACGGCTGGCAAGATTGCG‑CACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGc                                            >  2:947722/1‑151 (MQ=255)
                                                                                                          ttttGCTCCCGGCC‑CGCGGCGGCGGCAAAGATTCCC‑CCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAAc                                              <  1:188734/147‑1 (MQ=21)
                                                                                                          ttttGCGCCCTGCG‑CGCAACGCCTGGAAAGTTTGCG‑CCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGa                                                  <  2:1190724/143‑1 (MQ=25)
                                                                                                           tttGCTGCCTGCG‑CGCTACGGCTGGCAAGATTGCG‑CACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTg                                          >  2:196792/1‑150 (MQ=255)
                                                                                                               ctgcctgcG‑CGCTACGGCTGGCAAGATTGCG‑CACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAAt                                     >  2:337466/1‑151 (MQ=255)
                                                                                                               ctgccggcG‑CGCTACGGCGGGCAAGTTTGCC‑CCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATa                                       <  2:1159353/149‑1 (MQ=37)
                                                                                                                  ccggcg‑cgcTACGCCGGCCAAGTTGCCC‑CCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTc                                  <  1:833259/148‑1 (MQ=25)
                                                                                                                   cggcg‑cgcTCCGGCGGGAAATATTGCC‑CCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCg                                 <  2:1114870/149‑1 (MQ=25)
                                                                                                                     gcg‑cgcTACGGCTGGCAAGATTGCG‑CACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGcc                               >  1:990677/1‑151 (MQ=255)
                                                                                                                            tACGGCTGGCAAGATTGCG‑CACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtttt                         >  2:121231/1‑151 (MQ=255)
                                                                                                                              cGGCTGGCAAGATTGCG‑CACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTTTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtttttt                       >  1:1093874/1‑151 (MQ=255)
                                                                                                                               ggCTGGAAAGTTTGCG‑CCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTc                                  <  2:650078/139‑1 (MQ=25)
                                                                                                                               ggCGGGCAAGTTTCCC‑CCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAttt                          <  2:1118305/147‑1 (MQ=25)
                                                                                                                               ggCGGGAAATATTGCG‑CACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAttttt                        <  1:603216/149‑1 (MQ=255)
                                                                                                                               ggCGGGAAAGTTTGCG‑CACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGc                                <  2:1073048/141‑1 (MQ=255)
                                                                                                                                gCTGGAAAGATTGCG‑CCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGc                     <  2:871569/151‑1 (MQ=255)
                                                                                                                                   ggaaaGTTTGCC‑CCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTg                      <  2:1128885/144‑1 (MQ=255)
                                                                                                                                   ggCAAGATTGCC‑CCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCg                    <  2:989266/149‑1 (MQ=255)
                                                                                                                                       aagtttgcccctttCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGc                 <  2:430970/137‑1 (MQ=255)
                                                                                                                                        gATTGCC‑CTCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCg                    <  2:855829/144‑1 (MQ=255)
                                                                                                                                             cg‑cACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCg         >  1:801942/1‑150 (MQ=255)
                                                                                                                                             cg‑cACTCCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCtggtgg           <  2:598095/148‑1 (MQ=255)
                                                                                                                                                cACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGAtc      >  2:793381/1‑151 (MQ=255)
                                                                                                                                                    tCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGCCGCGAGCCGCTTCCACTGCAGCATTCAGGGCCAGAATCTTCGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCt  >  2:825666/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
GGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCG‑CGCTACGGCTGGCAAGATTGCG‑CACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCT  >  NC_000913/1969647‑1969942

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