Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913366,5260TG44.8% 80.4 / 26.3 58T329P (ACC→CCC) lacIlactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (12/20);  new base G (26/0);  total (38/20)
Fisher's exact test for biased strand distribution p-value = 1.34e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACA  >  NC_000913/366382‑366670
                                                                                                                                                |                                                                                                                                                
cTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTATTGGGGGCCAGGGGGGgttt                                                                                                                                            >  1:183372/1‑151 (MQ=255)
 tGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCCGGGGGGgtttt                                                                                                                                           >  1:372605/1‑151 (MQ=255)
 tGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGAGGAGGGGGGTTGTTTATTGGGGGCCCGGGGGGgtttt                                                                                                                                           >  1:242329/1‑151 (MQ=255)
  ggggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttc                                                                                                                                          <  2:809755/151‑1 (MQ=255)
   gggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttct                                                                                                                                         <  2:416555/151‑1 (MQ=255)
     gTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGTATTTGGGGGCAGGGGGGGttttcttt                                                                                                                                       >  1:173655/1‑151 (MQ=255)
         cTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCAGGGGGGGttttttttt                                                                                                                                      >  2:119834/1‑148 (MQ=255)
           aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTAcc                                                                                                                                  >  2:117002/1‑150 (MQ=255)
           aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCAcc                                                                                                                                  <  1:301821/150‑1 (MQ=255)
           aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCa                                                                                                                                 <  2:315394/151‑1 (MQ=255)
            aTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCCGGGGGGGTTTTTTTTTTAcc                                                                                                                                  >  2:721942/1‑149 (MQ=255)
               agtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGGCAGGGGGGGTTTTCTTTTTACCa                                                                                                                                 >  2:490516/1‑147 (MQ=255)
                 tgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCAGGGGGGGTTTTCTTTTCACCAGTg                                                                                                                              >  1:766251/1‑148 (MQ=255)
                  gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTGCGGGATTGGCGGCCGGGGGGtttttttttttccccggt                                                                                                                              >  1:304033/1‑139 (MQ=255)
                     cTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACgg                                                                                                                        <  2:739379/150‑1 (MQ=255)
                      tAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTTTTTTGCGCCCGGGGGGGTTTTTTTTTTTCCCCTGAGACagg                                                                                                                       >  1:706962/1‑150 (MQ=255)
                       aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACggg                                                                                                                       >  1:304065/1‑149 (MQ=255)
                        aCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTGGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACggg                                                                                                                       >  2:774333/1‑148 (MQ=255)
                           cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACa                                                                                                                  <  1:122128/150‑1 (MQ=255)
                            acaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCAGGGGGGGTTTTCTTTTTACCAGTGAGACGGGCAACa                                                                                                                  >  2:820381/1‑149 (MQ=255)
                                  aTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGGGAAGGGGaa                                                                                                                    >  2:743151/1‑141 (MQ=255)
                                  aTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTg                                                                                                          <  2:758506/151‑1 (MQ=255)
                                           cgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCAcc                                                                                                 <  1:826159/151‑1 (MQ=255)
                                            gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTTGTATTGGGCGGCCGGGGGGGTTTTTTTTTTTCCAGTGAGAAGGGGAACAGCCGGTTTccccttcccc                                                                                                >  1:921428/1‑144 (MQ=255)
                                                 cTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGGGAGACGGGGAACAACTGATTGCCCCTCACcccc                                                                                              >  1:772296/1‑148 (MQ=255)
                                                  tGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCAGGGGGGGTTTTTTTTTTACCCGTGAGAAGGGCAACAGCTGATTGCCCTTCACcccc                                                                                              >  1:818920/1‑147 (MQ=255)
                                                       gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGGATTTGGCGGCCGGGGGGGTTTTTTTTTTACCAGTGGGGAGGGCAACAGCTGATTGGCCCTCACCCCCCGGccc                                                                                        >  1:846690/1‑148 (MQ=255)
                                                           tCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgagag                                                                                  <  1:735758/150‑1 (MQ=255)
                                                                        ccTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCg                                                                    <  2:796374/151‑1 (MQ=255)
                                                                           gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGt                                                                  <  2:346991/150‑1 (MQ=255)
                                                                                        cATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGGGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCCGGCCCTGAGAGAGTTGCAGCAAACGGTCCACGCTGGgtt                                                      >  2:171857/1‑149 (MQ=255)
                                                                                         aTTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCCGGGGGGGTTTTTTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTTACCCCCCGGCCCTGAGAGAGGTGCAGCAAACGGTCCCCCCTGGGTTGcc                                                   >  2:717601/1‑151 (MQ=255)
                                                                                                aaTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGCCCGGGGGGGTTTTCTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCCTCACCCCCTGGCCCCGCGAGAGGTGCAGCAAGCGGGGCCCGCCGGGTTTCCCccgc                                              >  1:260250/1‑149 (MQ=255)
                                                                                                   cGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCAGGGGGGGTTTTTTTTTCACCAGGGAGACGGGGAACAGCTGATTGCCCCTCACCCCCCGGCCCTGGGAGAGGTGCAGCAAGCGGTCCACGCCGGGTTTCCCCCGCAgg                                           >  1:813305/1‑149 (MQ=255)
                                                                                                      ccAACGCGCGGGGAGAGGCGGTTTTTGTATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTTCCCCTCACCCCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCCCGCCGGGTTGCCCCCGCAGGCg                                         >  1:707090/1‑148 (MQ=255)
                                                                                                         aCGCGCGGGGAGAGGCGGTTTTCGTATTTGGCGCCAGGGGGGTTTTTCTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCCGGCCCTGAGAGAGGTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaaaa                                     >  2:840543/1‑149 (MQ=255)
                                                                                                            cgcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTg                                <  2:962815/151‑1 (MQ=255)
                                                                                                            cgcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTg                                <  1:378587/151‑1 (MQ=255)
                                                                                                                  gagaGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAt                          >  1:665206/1‑151 (MQ=255)
                                                                                                                    gagGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtgg                        >  2:730664/1‑151 (MQ=255)
                                                                                                                     agGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtg                         >  1:491264/1‑149 (MQ=255)
                                                                                                                       gCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtgg                     <  2:312822/151‑1 (MQ=255)
                                                                                                                          gTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                  <  2:751147/151‑1 (MQ=255)
                                                                                                                          gTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                  <  2:128912/151‑1 (MQ=255)
                                                                                                                            ttGCGTATTGGGCGCCAGGGGGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAc                >  1:261152/1‑151 (MQ=255)
                                                                                                                              gCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg               <  1:513051/150‑1 (MQ=255)
                                                                                                                                gTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcg            <  2:14332/151‑1 (MQ=255)
                                                                                                                                  aTTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg          >  1:41694/1‑151 (MQ=255)
                                                                                                                                   ttGCGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGATGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGCTTGCCCCAGCAGGCGACAATCCTGTTTGATGGTGGTTAAcggcgg           >  1:393419/1‑149 (MQ=255)
                                                                                                                                    tGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg          >  1:208342/1‑149 (MQ=255)
                                                                                                                                      ggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat      <  2:432223/151‑1 (MQ=255)
                                                                                                                                       gggcCAGGGTTGTTTTTCTTTTCCCCAGTGAGACGGGCAACAGCTGATTGCCCTTCCCCGGCTGGCCCTGAGGGAGTTGCCGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGAtata     <  2:801724/149‑1 (MQ=255)
                                                                                                                                        cgcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGAtata     >  2:825038/1‑150 (MQ=255)
                                                                                                                                         gcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa    >  2:771645/1‑150 (MQ=255)
                                                                                                                                         gcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa    >  1:626750/1‑150 (MQ=255)
                                                                                                                                           cAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa    >  2:429768/1‑148 (MQ=255)
                                                                                                                                            aGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa    >  1:314178/1‑147 (MQ=255)
                                                                                                                                            aGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa  >  2:33548/1‑149 (MQ=255)
                                                                                                                                              ggtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGGTATaaaa  >  1:415802/1‑145 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACA  >  NC_000913/366382‑366670

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