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→GGGtapmethyl‑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

CACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCT  >  NC_000913/1969650‑1969928
                                                                                                                                                 |                                                                                                                                      
cACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGcc                                                                                                                                  >  1:1061944/1‑151 (MQ=255)
    gAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACt                                                                                                                               >  1:638991/1‑150 (MQ=255)
     aaCTGACAATATCGATCATGGTCGCGGCGGCGTTATTCCCCAGTTTCGACCCCTGCTGGACACGATTGACTGCCTCTTCGATCAGCCCTTTGACCTCTTTGGCTGCCTGCGCGTTACGGCTGGAAAGATTGCGCACTTCCCCTGCCACTa                                                                                                                              <  2:89379/150‑1 (MQ=255)
     aaCTGACAATATCGATCATGGTCGCGGCGGCGTTATTCCCCAGTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCCCCTGCCACTa                                                                                                                              <  2:369863/150‑1 (MQ=255)
     aaCTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTa                                                                                                                              <  1:3675/150‑1 (MQ=255)
           cAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCaa                                                                                                                        >  1:22428/1‑150 (MQ=255)
            aaTATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCaaa                                                                                                                       >  2:687470/1‑150 (MQ=255)
                     cATGGTCGCGGCGGCGTTATTCCCCAGTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTGCGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCAc                                                                                                                  <  1:1008976/146‑1 (MQ=255)
                             cggcggcgTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTcc                                                                                                     >  1:791729/1‑151 (MQ=255)
                                cggcgTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGACTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGgc                                                                                                  >  2:1258226/1‑151 (MQ=255)
                                 ggctTTATTCCCCAGTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTTTTTTGCTGCCTGCGCGCTCCGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGgcg                                                                                                 <  1:1370729/151‑1 (MQ=255)
                                         tCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGc                                                                                         >  1:294663/1‑151 (MQ=255)
                                              atttccGACCCCTCCGGGACCCGATTGCCGTCCTCTTCCCCCAGCCCTTTGACCTCTTTCGCTCCCCGCCCGTGTCGGCTGAAAAGTTCGCTCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTcc                                                                                     <  1:1101455/145‑1 (MQ=255)
                                                ttccGACCCCTGCTGGACCCGATTGCCTGACTCTTCGATCAGCCTTTTGATCTTTTTTGCTGCCTGCGCGCTACGGCTGGCAAGTTTGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCa                                                                                    <  1:1049540/146‑1 (MQ=255)
                                                   cGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGcc                                                                               >  2:922201/1‑151 (MQ=255)
                                                     aaCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGc                                                                             >  1:253456/1‑151 (MQ=255)
                                                     aaCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGc                                                                             >  2:537826/1‑151 (MQ=255)
                                                     aaCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTACACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGc                                                                             >  1:1145015/1‑151 (MQ=255)
                                                       cccTGCGGGACACTTTTGACTCACTCTCCGATCACCCTTTTGATCTTTTTTGCTGCCGGCGCGCTACGGGGGGAAAGTTGGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAt                                                                           <  2:781886/151‑1 (MQ=255)
                                                        ccTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAtt                                                                          >  1:1337963/1‑151 (MQ=255)
                                                          tgctgGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAtt                                                                          >  1:595286/1‑149 (MQ=255)
                                                            ctgGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAgg                                                                      >  2:1018330/1‑151 (MQ=255)
                                                                acacGATTGCCTCACTCTTCGACCACCCTTTTTTCCTTTTTTGCGCCCGGCGCGCTACGGCGGGAAAGTTTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCa                                                                        <  2:1438599/145‑1 (MQ=255)
                                                                acacGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGcc                                                                   >  2:35603/1‑150 (MQ=255)
                                                                   ataTTGTCTTACTTTCCGATAAGCCTTTTTATCTCTTTTTCTCCCGGCCCGTTACGGCGGGAAAGATTGCGCTCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGa                                                                <  1:1243810/148‑1 (MQ=255)
                                                                      ttGACTGACTCTTCGATCCCCCTTTTTATCTTTTTTCCTGCCTGCGCCCTCCGGCGGGAAAGTTTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAat                                                              <  2:57547/149‑1 (MQ=25)
                                                                       tgactgacTCTTCGATCAGCCCTTTGTTCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATAtt                                                           >  2:576276/1‑151 (MQ=255)
                                                                            gcctctTCGATCAGCCCTTTGATCTTTTTTGCTGCCTGCGCGCTACGGCTGGAAAGTTTGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGt                                                        <  1:1269211/147‑1 (MQ=255)
                                                                            gactCTCCGACCACCCCTTTGATCTTTTTTGCTGCCGGCGCGCTCCGGCGGGCAGGATTGCCCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATAtt                                                           <  1:1131836/146‑1 (MQ=25)
                                                                               tctTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTgg                                                    >  1:135876/1‑150 (MQ=255)
                                                                                 ttCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGaa                                                  >  1:68741/1‑150 (MQ=255)
                                                                                   cgtccAGCCTTTTTTTCTCTTTTGCTCCGTCCCCTCTGCGGCGGGAAAGATTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTg                                                     <  2:123844/141‑1 (MQ=18)
                                                                                       cgcccttttGTTCTTTTTTCCTGCCGGCGCGCTACGGCGGGAAAGATTGCCCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCa                                              <  1:746998/142‑1 (MQ=25)
                                                                                         ccccTTTGATCTTTTTTGCTCCCCGCGCGCTACGGCTGGAAAGATGGCGCACTCCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTcc                                         <  2:512357/150‑1 (MQ=37)
                                                                                          ccttttattcttTTTTCCTGCCTGCGCGCTCCGGCGGAAAAGTTTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTcc                                         <  1:564896/142‑1 (MQ=21)
                                                                                          cccTTTGCTCTCTTTTGCTCCCCGCGCGCTACGGCGGCCAAGAGTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCg                                        <  1:335743/151‑1 (MQ=25)
                                                                                             tttgtttttttttGCTCCCCCCGCCCTCGGGCGGGAAAGTTTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTc                                      <  2:214263/141‑1 (MQ=21)
                                                                                               tgttcttTTTTGCTGCCTGCGCCCTCCGGGTGGCAAGTTTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGa                                    <  2:904588/147‑1 (MQ=25)
                                                                                               tGATCTTTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGa                                    <  2:58266/150‑1 (MQ=255)
                                                                                                  tcttTTTTCCTCCCTGCGCGCTCCGGCTGGAAAGTTTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGAt                                   <  1:783830/148‑1 (MQ=25)
                                                                                                     ctTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGc                              >  1:192492/1‑150 (MQ=255)
                                                                                                      ttttGCGGCCTGCGCGCTCCGGCGGGCAAGATCGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAAc                                <  2:1053507/147‑1 (MQ=25)
                                                                                                      tttcGCCGCCGGCGCGATACGGCTGGCACGTTGCCCCCTTTCCCCTGCCACTACCGCAAACCCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATa                                  <  2:568060/145‑1 (MQ=21)
                                                                                                       tttGCTCCCTGCGCCCTCGGGCGGGAAAGTTTGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATaa                                 <  2:635319/145‑1 (MQ=25)
                                                                                                          gctgcCTGCGCGCTCCGGCGGGAAAGTTTGCGCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATa                         <  1:484421/150‑1 (MQ=37)
                                                                                                           ctgcctgcGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGCTAACGCTGAt                          >  2:516886/1‑148 (MQ=255)
                                                                                                            tgcctgcGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATg                      >  2:177304/1‑151 (MQ=255)
                                                                                                              cctgcGCGCTACGGCTGACAAGATCGCCCACTCCACCTGCCACTCCCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTc                    <  1:140021/151‑1 (MQ=255)
                                                                                                              cctccgcgcccgggCGGGCAAGTTCCCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGt                     <  2:219288/138‑1 (MQ=25)
                                                                                                               ctgcGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCg                   >  1:354646/1‑151 (MQ=255)
                                                                                                               ctgcGCGCTACGGCTGGCAAGATTGCGCACGTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCg                   >  2:364363/1‑151 (MQ=255)
                                                                                                                 gcgcgcTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGcc                 >  2:1006150/1‑151 (MQ=255)
                                                                                                                       tcgggCGGGCAAGATTGCGCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtt             <  1:122968/146‑1 (MQ=255)
                                                                                                                       accggcggcaaaGTTTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCa                <  2:804253/136‑1 (MQ=25)
                                                                                                                        aCGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAttttt          >  2:426356/1‑151 (MQ=255)
                                                                                                                          ggCGGGCAAGTTTCCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtttttt         <  2:1158502/150‑1 (MQ=25)
                                                                                                                           gCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGc       >  2:553108/1‑151 (MQ=255)
                                                                                                                                cAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCt  >  1:327750/1‑151 (MQ=255)
                                                                                                                                cAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCt  >  2:93067/1‑151 (MQ=255)
                                                                                                                                 aaGATTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCt  <  1:35608/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                      
CACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCT  >  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: