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

TCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTG‑CTGGACACGA‑TTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGT  >  NC_000913/1969649‑1969931
                                                                                                                                                    |                                                                                                                                         
tCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTG‑CTGGACACGA‑TTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGc                                                                                                                                      >  2:991184/1‑151 (MQ=255)
  aCCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTG‑CTGGACACGA‑TTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGcc                                                                                                                                     >  1:791278/1‑150 (MQ=255)
    cGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTG‑CTGGACACGA‑TTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACt                                                                                                                                  >  2:819826/1‑151 (MQ=255)
      aaCTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTG‑CTGGACACGA‑TTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCCCTTCACCTGCCACTAc                                                                                                                                >  2:648382/1‑151 (MQ=255)
            cAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTG‑CTGGACACGA‑TTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCaa                                                                                                                           >  1:1041648/1‑150 (MQ=255)
             aaTATCGATCAGGGTCGCGGCGGCGTTATTCCCCAGTTCCGACCCCTG‑CTGGCCACGT‑TTGACTGACCCTTCGACCCGCCCTTTGACCTTTTTTGCTCCCGGCCCGCTACGGCTGGCAAGA‑TGGCCCCCTTCACCTGCCACTACCGCaaaa                                                                                                                         <  2:400235/151‑1 (MQ=255)
                atCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTG‑CTGGACACGA‑TTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCa                                                                                                                      <  1:1022009/151‑1 (MQ=255)
                  cGATCATGGTCGCGGGGGCTTTATTCACCAGTTTCGACCCCTG‑CTGGACACGA‑TTGACGGACTCTTCGATCAGCCCTTTGATCTTTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TGGCGCCCTCCCCCTGCCACTACCGCAAAACCACg                                                                                                                    <  1:628259/151‑1 (MQ=255)
                   gATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTG‑CTGGACACGA‑TTGACGGACTATTCGATCAGCCCTTTGATCTCTTGTGCTGCCTGCGCGCTACGGCTGGCAAGG‑TTGTGCACTTCATCTGCCACTACCGCAAAACCAc                                                                                                                     >  2:1177452/1‑149 (MQ=255)
                             gcggcggcgTTATTCACCAGTTTCGAACCCTG‑CTGGACACGA‑TTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCtctc                                                                                                         >  1:420809/1‑151 (MQ=255)
                              cggcggcgTTATTCACCAGTTTCGAACCCTG‑CTGGACACGA‑TTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTcc                                                                                                        >  2:561939/1‑151 (MQ=255)
                               ggcggcgTTATTCACCAGTTTCGAACCCTG‑CTGGACACGA‑TTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTcc                                                                                                        >  1:492647/1‑150 (MQ=255)
                                gcggcgTTATCCACCA‑TTTCCAACCCCGCCGGGACACGT‑TTGACTGTCTCTTCGTTCACCCCTTTGCTCTCTTTCCCTGCCTGCCCGCTGGGGCTGGCAAGA‑GCCCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGctct                                                                                                          <  2:362916/147‑1 (MQ=255)
                                      ttATTCCCCAGTTTCGACCCCGG‑CTGGCCACGT‑TTGACTGCCTCTTCGACCAGCCCTTTGATCTCTTTTGCTGCCTGCGCCCTCCGGCGGGCAAGA‑TTGCCCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGcgcg                                                                                                  <  1:705723/149‑1 (MQ=255)
                                        aTTCACCAGTTTCGAACCCTG‑CTGGACACGA‑TTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGcc                                                                                              >  2:709522/1‑151 (MQ=255)
                                           cACCAGTTTCGAACCCTG‑CTGGACACGA‑TTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCt                                                                                           >  2:61961/1‑151 (MQ=255)
                                             ccAGTTTCGACCCCTG‑CTGGACACTA‑TTGACTGCCTCTTCGATCACCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGTTGCCAAGT‑TTGCCCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGcc                                                                                              <  1:806239/146‑1 (MQ=255)
                                                 tttCGACCCCTG‑CGGGACACGT‑TTCACTGCCTCTCCAACCACCCTTTTGATCTTTTTTGCTCCCGGCGCCCGGCGCCGGGAAAGA‑TGGCGCTCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCa                                                                                       <  1:747363/149‑1 (MQ=255)
                                                 tttCGAACCCTG‑CTGGACACGA‑TTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACt                                                                                     >  2:634911/1‑151 (MQ=255)
                                                 ttccaaccccTG‑CTGGAC‑CGATTTGACTGCCTCTCCGCCCACCCCTTGGATCTTTTTTCCGCCCGGCGCCCCCGGGTTGCCAAGT‑TGGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCt                                                                                           <  1:15789/138‑1 (MQ=255)
                                                        cccTG‑CTGGACACGA‑TTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAt                                                                              >  1:112913/1‑151 (MQ=255)
                                                        cccTG‑CGGGACAGGT‑TTGAGTCACTCTCCGTTCGCCCTTTTGACCTCTTTTGCTGCCGGCGCCCTACGGCTGGCAAGA‑TGGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTg                                                                                    <  1:1105773/145‑1 (MQ=255)
                                                               tgGACGCGT‑TTACCTGTCTCTCCGACCAGCCTTTTGACCTTTTTTGCTGCCGGCCCGGTAGGGCGGGAAAGT‑TGCCGCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTc                                                                            <  1:444928/147‑1 (MQ=255)
                                                                ggACACGT‑TTGACTGCCTCTTCGTGCAGCCTTTTGACCTTTTTGGCTGCCGCCGCGCTAGGGCTGGCAAGA‑TTGCCCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAgg                                                                         <  1:214084/149‑1 (MQ=255)
                                                                   cacGA‑TTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAg                                                                          >  2:602278/1‑145 (MQ=255)
                                                                      gA‑TTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAat                                                                 >  1:850839/1‑151 (MQ=255)
                                                                         ttgccgaactCTTCGACCAGCCTTTTGTTCTTTTTTGCGGCCTGCGCGCGCCGGCGGGCAAAT‑TTGCCCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCa                                                                     <  1:368050/138‑1 (MQ=255)
                                                                         ttaactgcctCTTCGATCACCCTTTTGATTTTTTTTGCGGCCTGCGCGCTACGGCTGGCAAGT‑TTGCCCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAata                                                                <  1:920556/147‑1 (MQ=37)
                                                                            actgcctCTTCGATCACCCTTTTGATCTTTTTGGCGGCCGGCGCGCTCCGGCGGGAAAGA‑TTGCGCCCTCCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAatat                                                               <  1:542141/148‑1 (MQ=37)
                                                                                actCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTg                                                        >  2:890426/1‑151 (MQ=255)
                                                                                 ctctTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTgg                                                       >  2:120143/1‑151 (MQ=255)
                                                                                    ttcgaccccccttttGATCTTTTTTGCGGCCGGCGCGCTCCGGCGGGAAAGA‑TTGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTc                                                          <  2:797672/133‑1 (MQ=25)
                                                                                    ttcaaccacccttttGATCTTTTTTGCCGCCTGCGCGCTCCGGCGGGCAAGA‑TTGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTgg                                                       <  1:467218/136‑1 (MQ=25)
                                                                                    ttCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGaa                                                     >  2:635336/1‑150 (MQ=255)
                                                                                      cGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGa                                                      >  1:96923/1‑147 (MQ=255)
                                                                                      cGATCACCCTTTTGATCTTTTTTCCCGCCGGCGCGCTCCGGCTGGCAATT‑TTGCGCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGc                                                  <  1:78174/151‑1 (MQ=25)
                                                                                       gATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGc                                                  >  1:813037/1‑150 (MQ=255)
                                                                                          cacccttttGTTCTTTTTTGCTCCCGGCGCGCTACGGC‑GGAAAGATTTCCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAAt                                               <  2:985448/144‑1 (MQ=21)
                                                                                          cAGCCCTTTGTTCTTTTTTGCGGCCTGCGCGCTGCGGCTGGAAAGA‑TTGCGCCCTCCCCCTGCCGCTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAAtt                                              <  1:526115/151‑1 (MQ=25)
                                                                                          cAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAAtt                                              >  1:948254/1‑151 (MQ=255)
                                                                                          cAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAAtt                                              >  2:261342/1‑151 (MQ=255)
                                                                                              ccTTTGATCTTTTTTGCTGCCTGCGCGCTCCGGCTGGCAAGA‑TCGCCCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTc                                             <  1:1082563/148‑1 (MQ=37)
                                                                                              ccTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCCTTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCg                                           >  1:598924/1‑150 (MQ=255)
                                                                                                tttGATCTTTTTTGCTGCCGGCGCGCTCCGGCGGGCAAGA‑TTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGt                                          <  1:305674/149‑1 (MQ=37)
                                                                                                     tctcTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAAc                                   >  1:520076/1‑151 (MQ=255)
                                                                                                         ttttGCTCCCTGCGCGCTACGGCTGGAAAGT‑TTCCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAAc                                   <  1:63577/147‑1 (MQ=37)
                                                                                                             gctgcCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGCAAGCAATTCCGTCGATAACGCTGATa                            >  2:1080972/1‑150 (MQ=255)
                                                                                                                gcctgcGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGt                        >  2:464754/1‑151 (MQ=255)
                                                                                                                 ccggcgcgcTCGGGCTGGAAGGT‑TTCCCCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGt                        <  2:332397/147‑1 (MQ=25)
                                                                                                                   tgcGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGc                     >  2:1023865/1‑151 (MQ=255)
                                                                                                                      gcgcTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCa                   >  1:1086764/1‑150 (MQ=255)
                                                                                                                      gcgcTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCa                   >  1:1020785/1‑150 (MQ=255)
                                                                                                                       cgcTACGGTGGGCAATT‑TTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCa                   <  1:888968/149‑1 (MQ=25)
                                                                                                                       cgcTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAt                 >  2:1049958/1‑151 (MQ=255)
                                                                                                                             ggCTGGAAAGA‑TGCCCCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAttttt             <  1:819861/149‑1 (MQ=25)
                                                                                                                                 ggaaaaa‑TTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGc          <  1:407919/142‑1 (MQ=25)
                                                                                                                                    aaGA‑TTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGa        <  1:76164/147‑1 (MQ=255)
                                                                                                                                     aatatcgcgcACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGc      <  1:847300/143‑1 (MQ=255)
                                                                                                                                      gA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCtggt  >  1:1091326/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                         
TCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTG‑CTGGACACGA‑TTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGT  >  NC_000913/1969649‑1969931

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