Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,969,7940AC45.0% 80.8 / 28.2 60G397G (GGT→GGG) tapmethyl‑accepting protein IV
Reads supporting (aligned to +/- strand):  ref base A (28/5);  new base C (0/27);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 2.43e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCT  >  NC_000913/1969650‑1969928
                                                                                                                                                 |                                                                                                                                      
cACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGcc                                                                                                                                  >  1:1061944/1‑151 (MQ=255)
    gAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACt                                                                                                                               >  1:638991/1‑150 (MQ=255)
     aaCTGACAATATCGATCATGGTCGCGGCGGCGTTATTCCCCAGTTTCGACCCCTGCTGGACACGATTGACTGCCTCTTCGATCAGCCCTTTGACCTCTTTGGCTGCCTGCGCGTTACGGCTGGAAAGATTGCGCAC‑TTCCCCTGCCACTa                                                                                                                              <  2:89379/150‑1 (MQ=255)
     aaCTGACAATATCGATCATGGTCGCGGCGGCGTTATTCCCCAGTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCCCCTGCCACTa                                                                                                                              <  2:369863/150‑1 (MQ=255)
     aaCTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTa                                                                                                                              <  1:3675/150‑1 (MQ=255)
           cAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCaa                                                                                                                        >  1:22428/1‑150 (MQ=255)
            aaTATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCaaa                                                                                                                       >  2:687470/1‑150 (MQ=255)
                     cATGGTCGCGGCGGCGTTATTCCCCAGTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTGCGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCAc                                                                                                                  <  1:1008976/146‑1 (MQ=255)
                             cggcggcgTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTcc                                                                                                     >  1:791729/1‑151 (MQ=255)
                                cggcgTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGACTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGgc                                                                                                  >  2:1258226/1‑151 (MQ=255)
                                 ggctTTATTCCCCAGTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTTTTTTGCTGCCTGCGCGCTCCGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGgcg                                                                                                 <  1:1370729/151‑1 (MQ=255)
                                         tCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGc                                                                                         >  1:294663/1‑151 (MQ=255)
                                              atttccGACCCCTCCGGGACCCGATTGCCGTCCTCTTCCCCCAGCCCTTTGACCTCTTTCGCTCCCCGCCCGTGTCGGCTGAAAAGTTCGCTCAC‑TTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTcc                                                                                     <  1:1101455/145‑1 (MQ=255)
                                                ttccGACCCCTGCTGGACCCGATTGCCTGACTCTTCGATCAGCCTTTTGATCTTTTTTGCTGCCTGCGCGCTACGGCTGGCAAGTTTGCGCAC‑TTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCa                                                                                    <  1:1049540/146‑1 (MQ=255)
                                                   cGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGcc                                                                               >  2:922201/1‑151 (MQ=255)
                                                     aaCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGc                                                                             >  1:253456/1‑151 (MQ=255)
                                                     aaCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGc                                                                             >  2:537826/1‑151 (MQ=255)
                                                     aaCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TACACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGc                                                                             >  1:1145015/1‑151 (MQ=255)
                                                       cccTGCGGGACACTTTTGACTCACTCTCCGATCACCCTTTTGATCTTTTTTGCTGCCGGCGCGCTACGGGGGGAAAGTTGGCCCCC‑TCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAt                                                                           <  2:781886/151‑1 (MQ=255)
                                                        ccTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAtt                                                                          >  1:1337963/1‑151 (MQ=255)
                                                          tgctgGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAtt                                                                          >  1:595286/1‑149 (MQ=255)
                                                            ctgGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAgg                                                                      >  2:1018330/1‑151 (MQ=255)
                                                                acacGATTGCCTCACTCTTCGACCACCCTTTTTTCCTTTTTTGCGCCCGGCGCGCTACGGCGGGAAAGTTTGCGCCC‑TTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCa                                                                        <  2:1438599/145‑1 (MQ=255)
                                                                acacGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGcc                                                                   >  2:35603/1‑150 (MQ=255)
                                                                   ataTTGTCTTACTTTCCGATAAGCCTTTTTATCTCTTTTTCTCCCGGCCCGTTACGGCGGGAAAGATTGCGCTC‑TTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGa                                                                <  1:1243810/148‑1 (MQ=255)
                                                                      ttGACTGACTCTTCGATCCCCCTTTTTATCTTTTTTCCTGCCTGCGCCCTCCGGCGGGAAAGTTTGCCCCC‑TCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAat                                                              <  2:57547/149‑1 (MQ=25)
                                                                       tgactgacTCTTCGATCAGCCCTTTGTTCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATAtt                                                           >  2:576276/1‑151 (MQ=255)
                                                                            gcctctTCGATCAGCCCTTTGATCTTTTTTGCTGCCTGCGCGCTACGGCTGGAAAGTTTGCGCAC‑TTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGt                                                        <  1:1269211/147‑1 (MQ=255)
                                                                            gactCTCCGACCACCCCTTTGATCTTTTTTGCTGCCGGCGCGCTCCGGCGGGCAGGATTGCCCAC‑TCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATAtt                                                           <  1:1131836/146‑1 (MQ=25)
                                                                               tctTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTgg                                                    >  1:135876/1‑150 (MQ=255)
                                                                                 ttCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGaa                                                  >  1:68741/1‑150 (MQ=255)
                                                                                   cgtccAGCCTTTTTTTCTCTTTTGCTCCGTCCCCTCTGCGGCGGGAAAGATTGCCCCC‑TCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTg                                                     <  2:123844/141‑1 (MQ=18)
                                                                                       cgcccttttGTTCTTTTTTCCTGCCGGCGCGCTACGGCGGGAAAGATTGCCCAC‑TCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCa                                              <  1:746998/142‑1 (MQ=25)
                                                                                         ccccTTTGATCTTTTTTGCTCCCCGCGCGCTACGGCTGGAAAGATGGCGCAC‑TCCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTcc                                         <  2:512357/150‑1 (MQ=37)
                                                                                          ccttttattcttTTTTCCTGCCTGCGCGCTCCGGCGGAAAAGTTTGCCCCC‑TCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTcc                                         <  1:564896/142‑1 (MQ=21)
                                                                                          cccTTTGCTCTCTTTTGCTCCCCGCGCGCTACGGCGGCCAAGAGTGCGCCC‑TTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCg                                        <  1:335743/151‑1 (MQ=25)
                                                                                             tttgtttttttttGCTCCCCCCGCCCTCGGGCGGGAAAGTTTGCCCCC‑TCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTc                                      <  2:214263/141‑1 (MQ=21)
                                                                                               tgttcttTTTTGCTGCCTGCGCCCTCCGGGTGGCAAGTTTGCCCCC‑TCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGa                                    <  2:904588/147‑1 (MQ=25)
                                                                                               tGATCTTTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCCCCC‑TCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGa                                    <  2:58266/150‑1 (MQ=255)
                                                                                                  tcttTTTTCCTCCCTGCGCGCTCCGGCTGGAAAGTTTGCCCCC‑TCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGAt                                   <  1:783830/148‑1 (MQ=25)
                                                                                                     ctTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGc                              >  1:192492/1‑150 (MQ=255)
                                                                                                      ttttGCGGCCTGCGCGCTCCGGCGGGCAAGATCGCGCCC‑TCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAAc                                <  2:1053507/147‑1 (MQ=25)
                                                                                                      tttcGCCGCCGGCGCGATACGGCTGGCACG‑TTGCCCCCTTTCCCCTGCCACTACCGCAAACCCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATa                                  <  2:568060/145‑1 (MQ=21)
                                                                                                       tttGCTCCCTGCGCCCTCGGGCGGGAAAGTTTGCGCCC‑TCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATaa                                 <  2:635319/145‑1 (MQ=25)
                                                                                                          gctgcCTGCGCGCTCCGGCGGGAAAGTTTGCGCAC‑TCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATa                         <  1:484421/150‑1 (MQ=37)
                                                                                                           ctgcctgcGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGCTAACGCTGAt                          >  2:516886/1‑148 (MQ=255)
                                                                                                            tgcctgcGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATg                      >  2:177304/1‑151 (MQ=255)
                                                                                                              cctgcGCGCTACGGCTGACAAGATCGCCCAC‑TCCACCTGCCACTCCCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTc                    <  1:140021/151‑1 (MQ=255)
                                                                                                              cctccgcgcccgggCGGGCAAGTTCCCGCCC‑TCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGt                     <  2:219288/138‑1 (MQ=25)
                                                                                                               ctgcGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCg                   >  1:354646/1‑151 (MQ=255)
                                                                                                               ctgcGCGCTACGGCTGGCAAGATTGCGCAC‑GTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCg                   >  2:364363/1‑151 (MQ=255)
                                                                                                                 gcgcgcTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGcc                 >  2:1006150/1‑151 (MQ=255)
                                                                                                                       tcgggCGGGCAAGATTGCGCAC‑TCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtt             <  1:122968/146‑1 (MQ=255)
                                                                                                                       accggcggcaaaGTTTGCCCCC‑TCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCa                <  2:804253/136‑1 (MQ=25)
                                                                                                                        aCGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAttttt          >  2:426356/1‑151 (MQ=255)
                                                                                                                          ggCGGGCAAGTTTCCCCCC‑TCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtttttt         <  2:1158502/150‑1 (MQ=25)
                                                                                                                           gCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGc       >  2:553108/1‑151 (MQ=255)
                                                                                                                                cAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCt  >  1:327750/1‑151 (MQ=255)
                                                                                                                                cAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCt  >  2:93067/1‑151 (MQ=255)
                                                                                                                                 aaGATTGCGCCC‑TTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCt  <  1:35608/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                      
CACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCT  >  NC_000913/1969650‑1969928

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: ‑