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

GGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCT  >  NC_000913/1969647‑1969942
                                                                                                                                                      |                                                                                                                                                    
ggTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCt                                                                                                                                                   >  1:1088159/1‑151 (MQ=255)
     ccGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGcc                                                                                                                                                <  2:481308/149‑1 (MQ=255)
       gAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACt                                                                                                                                             >  2:894168/1‑150 (MQ=255)
        aaCTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTAc                                                                                                                                           >  1:585344/1‑151 (MQ=255)
                        cATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGcccc                                                                                                                           >  1:713586/1‑151 (MQ=255)
                          tGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTg                                                                                                                         >  2:360208/1‑151 (MQ=255)
                           ggcggcggcggcgTTATTCCCCATTTTCGACCCCGGCTGGACACGATTGACTGCCTCTTCGACCAGCCTTTTGATCTCTTTTGCTGCCTGCGCGCTCGGGTTGGCAAGATTGCGCACTCCCCCTGCCACTACCGCAAAACCACGCCCCt                                                                                                                          <  2:826951/145‑1 (MQ=255)
                           ggTCGGGGCGGGGTTTTTCACCATTTCCGACCCCTGTGGGACACAATTGACTAACTCTTCGATCCCCCCTTTGATCTTTTTTCCTCCCTGCCCGCTCCGGCTGGCAAGATTGCGCACTCCCCCTGCCACTACCGCAAAACCACGCCCCt                                                                                                                          <  2:605678/149‑1 (MQ=255)
                               gcggcggcgTTATTCCCCAGTTCCGACCCCTGTTGGACCCGATTGCCTGACTCTTCGACCACCCCTTTGATCTTTTTGGCTGCCCGCGCGCTACGCCGGGCAATTTTGCGCCCTTCGCCTGCCACTACCGCAAAACCACGCCCCTGctct                                                                                                                     <  1:485061/150‑1 (MQ=255)
                               gcggcggcgTTATTCACCAGTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTTTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGctct                                                                                                                     <  2:549180/150‑1 (MQ=255)
                                    ggggTTATTCCCCATTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTTTTTTGCTGCCTGCGCGCTCCGGCTGGCAAGATTGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGgc                                                                                                                <  2:788742/147‑1 (MQ=255)
                                       gTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGcgcg                                                                                                             >  2:459370/1‑150 (MQ=255)
                                        ttATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGcgcg                                                                                                             >  2:562140/1‑149 (MQ=255)
                                           ttCACCAGTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCTTTTGATCTCTTTGGCTGCCTGCGCGCTACGGCTGGAAAGATTGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGcc                                                                                                         <  2:830005/150‑1 (MQ=255)
                                            tCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGc                                                                                                       >  1:559210/1‑151 (MQ=255)
                                             cACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCt                                                                                                      >  1:1182016/1‑151 (MQ=255)
                                                cAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTcc                                                                                                   >  2:618384/1‑151 (MQ=255)
                                                   tttCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACt                                                                                                >  1:1124992/1‑151 (MQ=255)
                                                    tccaaccccTGCTGGCCACTTTTGCCTGCCTCTTCGACCCCCCTTTTGATCTTTTTTGCTGCCTGCGCGCTCCGGCGGGAAAGTTTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACt                                                                                                <  1:28476/144‑1 (MQ=255)
                                                        accccTGCTGGACACGATTGACTGCCTCTTCGATCACCCTTTTGTTCTTTTTTGCTGCCTGCGCGCTCGGGCGGGCAAGATGGCGCATTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTg                                                                                               <  2:874652/145‑1 (MQ=255)
                                                           ccTGCTGGCCACGATTGACTGACTCTTCGATCCGCCCTTTGCCCTTTTTGGCTCCCCGCCCGCTGGGGCGGGCAAGATCGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCg                                                                                            <  1:66502/147‑1 (MQ=255)
                                                              gctgGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTc                                                                                       >  2:1198011/1‑149 (MQ=255)
                                                                   acacGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGcc                                                                                 >  1:424456/1‑150 (MQ=255)
                                                                    cacGTTTGACTGCCTCTCCGATCCCCCCTTTGATCTTTTTGGCTGCCGGCGCGCTACGGCGGGAAAGTTGGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAg                                                                                     <  1:929262/145‑1 (MQ=255)
                                                                      cGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGa                                                                              >  2:389298/1‑150 (MQ=255)
                                                                        aTTGACTGACTCTGCGTCCAGCCTTTTAACCTCTTTTGCGGCCGCCGCGCTCGGGCGGGAAAGATCGCCCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAg                                                                               <  2:756706/147‑1 (MQ=255)
                                                                                  ttctCGATCGCCCTTTGAATCTTTTTGCCTGCCGGGGCGCTGCGGGTGGAAAGTTTGCCCTCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTg                                                                        <  2:528097/141‑1 (MQ=18)
                                                                                    tccaacccgCCTTTTGATCTCTTTGCCGCCCGCGCCACTACGGCGGGCAGGTTGCCGCCCGCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCt                                                                    <  1:52477/138‑1 (MQ=16)
                                                                                      cgctCACCCTTTTCATCTTTTTTGCTCCCGGCGCCCGCGGGCTGGAAAGATTCCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTg                                                                   <  1:274234/142‑1 (MQ=18)
                                                                                        aTCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGCTTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAg                                                              >  1:991211/1‑148 (MQ=255)
                                                                                               ctttggtcctctTTTGCCCCCTCCGCCCTACGGCCGCAAATTTGCGCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAg                                                              <  1:869155/134‑1 (MQ=18)
                                                                                                tttGATCTTTTTTGCTCCCCGCGCGCTCCGGCTGCCAAGATTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGt                                                     <  2:218579/149‑1 (MQ=25)
                                                                                                tttGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCg                                                   >  2:896199/1‑151 (MQ=255)
                                                                                                 tttgcttttttttGTTGCCGCCGCCCTGCGGCGGGAAGGTTTGCCCTCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTc                                                    <  2:429743/141‑1 (MQ=21)
                                                                                                 ttGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTc                                                    >  2:1059137/1‑149 (MQ=255)
                                                                                                  tGATCTTTTTTGCTCCCGGCGCGCTCGGGCTGGCAAGATTGCGCATTCCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGa                                                  <  2:1088195/150‑1 (MQ=37)
                                                                                                   gATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGAt                                                 >  2:222477/1‑150 (MQ=255)
                                                                                                        tctTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGc                                            >  2:947722/1‑151 (MQ=255)
                                                                                                          ttttGCTCCCGGCCCGCGGCGGCGGCAAAGATTCCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAAc                                              <  1:188734/147‑1 (MQ=21)
                                                                                                          ttttGCGCCCTGCGCGCAACGCCTGGAAAGTTTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGa                                                  <  2:1190724/143‑1 (MQ=25)
                                                                                                           tttGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTg                                          >  2:196792/1‑150 (MQ=255)
                                                                                                               ctgcctgcGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAAt                                     >  2:337466/1‑151 (MQ=255)
                                                                                                               ctgccggcGCGCTACGGCGGGCAAGTTTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATa                                       <  2:1159353/149‑1 (MQ=37)
                                                                                                                  ccggcgcgcTACGCCGGCCAAGTTGCCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTc                                  <  1:833259/148‑1 (MQ=25)
                                                                                                                   cggcgcgcTCCGGCGGGAAATATTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCg                                 <  2:1114870/149‑1 (MQ=25)
                                                                                                                     gcgcgcTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGcc                               >  1:990677/1‑151 (MQ=255)
                                                                                                                            tACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtttt                         >  2:121231/1‑151 (MQ=255)
                                                                                                                              cGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTTTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtttttt                       >  1:1093874/1‑151 (MQ=255)
                                                                                                                               ggCTGGAAAGTTTGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTc                                  <  2:650078/139‑1 (MQ=25)
                                                                                                                               ggCGGGCAAGTTTCCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAttt                          <  2:1118305/147‑1 (MQ=25)
                                                                                                                               ggCGGGAAATATTGCGCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAttttt                        <  1:603216/149‑1 (MQ=255)
                                                                                                                               ggCGGGAAAGTTTGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGc                                <  2:1073048/141‑1 (MQ=255)
                                                                                                                                gCTGGAAAGATTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGc                     <  2:871569/151‑1 (MQ=255)
                                                                                                                                   ggaaaGTTTGCCCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTg                      <  2:1128885/144‑1 (MQ=255)
                                                                                                                                   ggCAAGATTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCg                    <  2:989266/149‑1 (MQ=255)
                                                                                                                                       aagtttgcccctttCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGc                 <  2:430970/137‑1 (MQ=255)
                                                                                                                                        gATTGCCCTCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCg                    <  2:855829/144‑1 (MQ=255)
                                                                                                                                             cgcACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCg         >  1:801942/1‑150 (MQ=255)
                                                                                                                                             cgcACTCCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCtggtgg           <  2:598095/148‑1 (MQ=255)
                                                                                                                                                cACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGAtc      >  2:793381/1‑151 (MQ=255)
                                                                                                                                                    tCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGCCGCGAGCCGCTTCCACTGCAGCATTCAGGGCCAGAATCTTCGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCt  >  2:825666/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
GGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCT  >  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: