Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8750AT42.1% 81.4 / 23.5 58N272K (AAT→AAAcynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (20/13);  new base T (0/24);  total (20/38)
Fisher's exact test for biased strand distribution p-value = 5.33e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTC  >  NC_000913/357725‑358006
                                                                                                                                                         |                                                                                                                                      
aCTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTa                                                                                                                                        <  1:739448/151‑1 (MQ=255)
     aCTCGTATGTCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCGGCCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACTCCATGTTCAAAAATGCCTTCGCGACGGCTGTCTGCCAGCTAttgtt                                                                                                                                   <  2:793562/151‑1 (MQ=14)
     aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                                                                                   >  1:957946/1‑151 (MQ=255)
        cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTcc                                                                                                                                <  1:1116491/150‑1 (MQ=255)
        cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTcc                                                                                                                                <  2:177564/150‑1 (MQ=255)
        cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTc                                                                                                                                 <  2:808279/149‑1 (MQ=255)
         gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgc                                                                                                                              >  1:1070631/1‑151 (MQ=255)
            ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc                                                                                                                           <  2:172441/151‑1 (MQ=255)
             gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc                                                                                                                           >  1:1172195/1‑150 (MQ=255)
                 gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                        <  1:1232236/149‑1 (MQ=255)
                  aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                    <  1:1301364/151‑1 (MQ=255)
                    aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                   <  1:1205811/150‑1 (MQ=255)
                    aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                  >  1:1303713/1‑151 (MQ=255)
                       cGTAGCGCATCAGGCAATTCCACCCGCTGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                 <  2:1219159/149‑1 (MQ=255)
                             gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc                                                                                                          >  2:554059/1‑150 (MQ=255)
                                tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCTCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctctc                                                                                                      >  1:1057746/1‑151 (MQ=255)
                                  aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                     >  1:1185853/1‑150 (MQ=255)
                                        ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCCCTCCAGGAGTg                                                                                              >  2:674692/1‑151 (MQ=255)
                                          ccACCCGCAGCCTTGTGTCGGCGGCTCCCGTGTTTCTTTCCCGCCAACACCCGCGCTTTTACCAACCGCCTGGCGAAAAAATCCCTCCGGGGGGGTTGTCTGCCAGTTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                            <  1:1324225/151‑1 (MQ=255)
                                           cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                          >  2:715539/1‑151 (MQ=255)
                                              ccGCAGCCTTGTCTCGGCGCCTCCTGTTATTTATTCCCGCAAACACCCGCGCATTAATCCAACGCCTTGTGAAAAATGTCTTCCGCGGCGCCTGTCTGCAATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagtag                                                                                                <  2:171586/143‑1 (MQ=255)
                                               cccagccatgtgCCGGCGGCTGCCGTGATCCATCCCCCCCAACACCCGCGCTTTTACCCACCCCGATGTAAAAAAGTTCGTCGGCGGCGCCTGCCGCCCCGCTATTTTCCCCGCCAAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                        <  2:1456042/141‑1 (MQ=255)
                                               cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                       >  1:547115/1‑150 (MQ=255)
                                                   gcccTGTGTCAGCGCCTACCGTTTTTCTTTCCCGCCAACACCCGCGCTTTTACCCACCGCCATGTAAAAAAATGCCTTCGCGGCGGCTTTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                   <  1:657477/148‑1 (MQ=255)
                                                      cTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                       >  1:585558/1‑143 (MQ=255)
                                                         tttgggCGTCTACCGTGATTCTTCCCCCCCAACACCCGCGCTTTTATCCAGCCCCATGGGAAAAAGTCCCTCCGCGGCGGCTGTGTCCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaa                                                                              <  1:664317/144‑1 (MQ=255)
                                                         tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggagggggaa                                                                             >  2:692297/1‑140 (MQ=255)
                                                              gcgcCTCCCGTGTTTCATTCCCCCCACCACCCGCGCTTTTCCCCACCCCTTTGTAAAAAAATCCTTCGGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                           <  2:98300/147‑1 (MQ=255)
                                                                ggCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaaaaaaattt                                                                       >  1:858129/1‑142 (MQ=255)
                                                                ggCTACCGTAATTCTTTCCCGCCAACACCCGGGCTTTTACCCACCCCCCTGTAAAAAAACGCCTGCGCGGGGGCTGCCTCCCAGTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                                                             <  1:15446/143‑1 (MQ=255)
                                                                gcctcccGTGTTTCTTTCCCGCCACCACCCCCGCTTTTACCAACCGCCAGGTCAAAAAAGCCTTCCGCGGCGGTTGCCTGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaa                                                                               <  2:230456/136‑1 (MQ=255)
                                                                gcctcccGTGATTAATTCCCGCAACCACCCGCGCTTTTACCCACCCCCATGTGAAAAAAGCCCTTCGCGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCt                                                                         <  1:949156/142‑1 (MQ=255)
                                                                 gCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGgaaaaaattt                                                                       >  1:181921/1‑141 (MQ=255)
                                                                  cccccGTGTTTCTTTCCCGCCAACACCCCCGATTTTACCCACCGCTATGTCAAAAAATCCTTCCGCGGCGGCTTTCTGCCAGTTTTTTTCCCGCCACAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAg                                                                    <  2:1217260/147‑1 (MQ=255)
                                                                  cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAGAAGaaaaaaactt                                                                       >  2:1323704/1‑140 (MQ=25)
                                                                    ccggttaTTAATTCCCGCCAACACCCCCCCATTTACCCACCCCATTGTGAAAAAATCCCTCGGCGGCGGCTTTCGGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGccc                                                                 <  2:1346563/145‑1 (MQ=255)
                                                                    aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAACTTTAAAccc                                                                 >  1:598663/1‑151 (MQ=255)
                                                                      cggtttttaTTTCCCCCAACAACCCGCGCATTTACCCACCCCCTTGTAAAAAATTCCCTCCGGGGCGGTTGTCGGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCtt                                                                        <  2:1254825/135‑1 (MQ=255)
                                                                               ttccccccAAAACCCGCCCTTTTACCCACCCCAATGTAAAAAAATCCCTGCGCGGGGGCTGTCTGCCAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                                                          <  1:592878/141‑1 (MQ=255)
                                                                               tccccGCCACCACCCGCTTTTTTACCCCACGCCATGTAAAAAAATCCTTCCGCGGGGGCTGTCTGCCAGTTTTTTTCCCGCCACACAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGACCGGCAtgttg                                                        <  2:608206/146‑1 (MQ=255)
                                                                                  cccccAACACCCGCGCTTTTACCCACCGCCTTGTGAAAAATTGCTTTGGGGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgt                                                     <  1:252997/146‑1 (MQ=255)
                                                                                   ccccaaccacccccgcTTTTACCCACCCCCATGGGAAAAATTCCCTGCGGGGCGGTTGTCGCCAATTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                                                          <  2:954952/131‑1 (MQ=255)
                                                                                       aacaccCGCGCATTTACCCACCCCCTTGTGAAAAAATTCCTTCGGGGCGGCTTTCTGCCATTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAt                                               <  1:1368108/150‑1 (MQ=255)
                                                                                             ccgcacttttccccACCGCCATGTAAAAAAGTGCATGGGCGGGGGCGGTCGGCGGTTTTTTTGCCCGCCACAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                                              <  2:688549/134‑1 (MQ=255)
                                                                                               gcgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAAAAAAAAAACTTTAACCCCTCCTGTTTTGTTGGAATccc                                            >  2:1309730/1‑144 (MQ=25)
                                                                                                       aTCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAATTTTAACCCCTCCATTTGTGTTGCGAACGCTGCTggggaaaa                                >  2:774083/1‑144 (MQ=25)
                                                                                                         ccaccgccagggaaaaaaaTGCCTGCGCGGCGCCTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                  <  2:1276166/132‑1 (MQ=255)
                                                                                                             ccccttgtaaaaaaaTGCTTCCGCGGCGGTTGCCGGCCAGTTTTTTTCCCCCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                              <  2:852803/136‑1 (MQ=255)
                                                                                                               ccagggaaaaaaaTGCCTCCGGGGCGGTTTTCGGCCAGTTTTTTTCCCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                          <  2:751548/140‑1 (MQ=255)
                                                                                                                   gtcaaaaaaTCCCTGCGGGGCGGTTGTCTCCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg                      <  1:62683/144‑1 (MQ=255)
                                                                                                                   gtaaaaaaaGGCCTCCGCGGGGCTTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg                      <  2:75230/144‑1 (MQ=255)
                                                                                                                    taaaaaaaTCCTTCCGCGGCGCTTGTCTGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                     <  2:476208/145‑1 (MQ=255)
                                                                                                                    taaaaaaaTCCCTCCGCGGCGGCTTTCTCCCATTTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATAGCTGCTGGTAACAATGTGGAAAGGGa                    <  1:97284/146‑1 (MQ=255)
                                                                                                                    aacAAAAATGCCTCCGCGGCGGTTGTCTGCCGGTTATTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg                  <  2:1143154/149‑1 (MQ=255)
                                                                                                                         aaaTGCCTGCGGGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg             <  2:588244/151‑1 (MQ=255)
                                                                                                                          aaTGCCTTCGGGGGCGCCTTTCTCCCGTCTTTTTTCCCCCCAAACAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg             <  2:1513976/151‑1 (MQ=255)
                                                                                                                               cTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTTAAACCCCCCTGTTTTGTGGGGATTGCTGCTGGGAACAAAGGGGGAAGGGGAgggggggg           >  2:1343207/1‑140 (MQ=255)
                                                                                                                                       ggcggcTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTc  <  1:437978/150‑1 (MQ=255)
                                                                                                                                       ggcggcTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGACGAGGAAAAGAATATTACCCGACATATTGTGTGGGGCTCGGCGCTGGGGACAATATTGAAAAGGGAGGGCGGGGAGTCAGc    >  2:166503/1‑148 (MQ=255)
                                                                                                                                                         |                                                                                                                                      
ACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTC  >  NC_000913/357725‑358006

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