Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913366,5260TG43.3% 113.9 / 18.8 67T329P (ACC→CCC) lacIlactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (5/33);  new base G (29/0);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 4.05e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCT  >  NC_000913/366378‑366676
                                                                                                                                                       |                                                                                                                                                      
aaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGATTTCGGA‑TTTGGC‑G‑CCAgggggg                                                                                                                                                      >  1:922358/1‑151 (MQ=255)
  gCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGA‑TTTGGG‑G‑CCCgggggg                                                                                                                                                      >  2:38792/1‑149 (MQ=255)
       gggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTttttct                                                                                                                                               <  2:1015882/151‑1 (MQ=255)
         gTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTttttcttt                                                                                                                                             <  2:541004/151‑1 (MQ=255)
          tGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTT‑TTGGGC‑G‑CCCGGGGGGtttttttttt                                                                                                                                            >  2:379529/1‑151 (MQ=255)
                    gtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTgag                                                                                                                                  <  1:54001/151‑1 (MQ=255)
                     tgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGACAGGCGGTTTTTGCA‑TTTGGCGG‑CCCGGGGGGtttttttt                                                                                                                                              >  2:697109/1‑139 (MQ=255)
                        gCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACgg                                                                                                                              <  1:409078/151‑1 (MQ=255)
                        gCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACg                                                                                                                               <  2:655972/150‑1 (MQ=255)
                         cTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCTGG‑TTGTGA‑GACCGGGGGGGtttttttttttccaagtggaaggg                                                                                                                              >  2:203772/1‑137 (MQ=255)
                          tAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGG‑G‑CCCGGGGGGGTTTTTTTTTTACCAGTGAGACggg                                                                                                                             >  1:61465/1‑150 (MQ=255)
                           aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGGA‑TTTGGG‑G‑GCAGGGGGGTTTTTTTTTTTAACAGTggggaggg                                                                                                                             >  2:716004/1‑142 (MQ=255)
                           aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACggg                                                                                                                             <  1:423530/149‑1 (MQ=255)
                            aCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCaa                                                                                                                          <  2:774735/151‑1 (MQ=255)
                             cTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGGG‑TTTGGC‑G‑CCCGGGGGGGTTTTTTTTTTACCAGTgagaacgggaa                                                                                                                          >  2:494618/1‑143 (MQ=255)
                               cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACa                                                                                                                        <  1:1051494/150‑1 (MQ=255)
                                acaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGc                                                                                                                      <  1:860252/151‑1 (MQ=255)
                                  aTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTg                                                                                                                    <  1:685944/151‑1 (MQ=255)
                                     aaTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGAtt                                                                                                                 <  1:610693/151‑1 (MQ=255)
                                      aTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTg                                                                                                                <  2:133651/151‑1 (MQ=255)
                                       ttgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGA‑TTTGGG‑G‑CCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGAtt                                                                                                                 >  1:424400/1‑149 (MQ=255)
                                        tgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGc                                                                                                               <  1:229529/150‑1 (MQ=255)
                                          cgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCt                                                                                                            <  2:793010/151‑1 (MQ=255)
                                               cgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGG‑G‑CCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAAAAGCTGATTGCCCttttcc                                                                                                       >  2:913622/1‑147 (MQ=255)
                                                    aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTG‑TTTGGG‑G‑CCCGGGGGGGTTTTCTTTTTTCCAGTGAGACGGGGAACAGCCGGTTGTCCTTCACcccc                                                                                                    >  2:432463/1‑149 (MQ=255)
                                                    aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTg                                                                                                  <  2:971660/151‑1 (MQ=255)
                                                             tttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTTCGGG‑TTTGGC‑G‑CCGGGGGGGGTTTTTTTTTTACCCGGGAGACGGGCGACACCTGGATGCCCCTCACCGCCCGGgcc                                                                                              >  2:373708/1‑146 (MQ=255)
                                                               tCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGC‑G‑CCCGGGGGGGTTTTTTTTTCACCCGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgag                                                                                          >  1:288967/1‑148 (MQ=255)
                                                                ccAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTTGTG‑TTTGGG‑G‑CCCGGGGGGGTTTTTTTTTTACCCGGGAGACGGGGCACAGGCGATTTCCCCTTACCCCCCGGCCCCtagag                                                                                        >  1:955970/1‑149 (MQ=255)
                                                                ccAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgagaga                                                                                       <  2:891824/150‑1 (MQ=255)
                                                                ccAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTGCGGA‑TTTGGC‑G‑GCCGGGGGGGTTTTTTTTTTACCCGGGAGACGGGCAACAGCTGATTGCCCTTCACCCCCCGGCCCTgagag                                                                                        >  1:461178/1‑149 (MQ=255)
                                                                  aGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTTGCATTTGGGG‑G‑CGGGGGGGGTTTTTTTTTTAACC‑TTGAGACGAGTAACAGGCGATTGCCCCTCACCGCCTGGCGCTGcgtgaggt                                                                                    >  1:679192/1‑149 (MQ=255)
                                                                  aGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTTCGGA‑TTGGGG‑G‑GCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGAGTGTCCTTCACCCCCCGGCCCCGagagaggt                                                                                    >  2:408045/1‑149 (MQ=255)
                                                                   gTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTg                                                                                   <  2:1107003/151‑1 (MQ=255)
                                                                     cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTGGGG‑G‑CCCGGGGGGGTTTTCTTTTTCCCCGTGAGACGGGCAACACCTGATTGCCCCTCACCGCCTGGCCCCGagggaggt                                                                                    >  1:727327/1‑146 (MQ=255)
                                                                     cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTGCGTA‑TTGGGG‑G‑CCGGGGGGGGTTTTTTTTTTCCCAGGGCGACGGGGAACAGCTGATTGGCCTTCACCGCCTGGTCCCGTGAGAGtt                                                                                    >  1:292328/1‑148 (MQ=255)
                                                                      gggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGG‑G‑CCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGGAACAGCTGATTGTCCCTTACCGCCTGGCCCCGggagag                                                                                      >  1:507328/1‑145 (MQ=255)
                                                                        gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTGGGG‑G‑GCAGGGGGGGTTTTTTTTTTACCAGGGAGAGGGGGAACAGCTGATTGCCCTTCACCCCCTGGCCCCGTGAGAGTTtcag                                                                                >  2:438078/1‑149 (MQ=255)
                                                                         aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCGCCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCaa                                                                             <  2:711976/151‑1 (MQ=255)
                                                                             cTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTGGGC‑G‑CCCGGGGGGGTTTTCTTTTTCCCAGGGAGACGGGCAACAGCTGATTGCCCTTTCCCCCCCCGCCCCGCGGGAGGTGCAGCAAACgg                                                                         >  2:513359/1‑151 (MQ=255)
                                                                                  gTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCa                                                                     <  1:367686/150‑1 (MQ=255)
                                                                                   tGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGC‑G‑CCAGGGGGGGTTTTCTTTTTACCAGTGAGACGGGCAAAAGCTGATTGCCCTTCACCCCCTGGCCCTGAGAGAGGTGCAGCAAGCGGTCCAc                                                                    >  2:328351/1‑150 (MQ=255)
                                                                                      cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGTGTA‑ATTGGC‑G‑CCAGGGGGGGTTTTCTTTTTACCAGGGAGACGGGCCACAGCTGATTGCCCTTCCCCGCCCGGCCCTCAGGGAGGTTCAGCAAACGGGCCACGc                                                                  >  2:491407/1‑149 (MQ=255)
                                                                                       aGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCTTT‑TTTGGG‑G‑CCCAGGGGGGTTTTCTTTTTTCCAGTGAGAGGGGCCAAAGCCGATTTCCCCTCACCCGCTAGCCCTTTGAGAGGTGGAGCAAACGGTccccaccgg                                                               >  1:428308/1‑144 (MQ=255)
                                                                                          tGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGC‑G‑CCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCCTCACCCCCTGGCCCTGAGGGAGGTGCAGCAAACGGGCCCCCCCGGttt                                                            >  1:283552/1‑151 (MQ=255)
                                                                                          tGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGATTGCGGA‑TTGGGG‑G‑CCCGGGGGGGTTTTCTTTTCCCCCGGGGGACGGGGAACAGCTGATTTCCCTTCACCGCCCGGCCCTTGGAGAGTTGCGGCAAGCGGGGCCCGCTGGgtt                                                            >  1:883949/1‑151 (MQ=255)
                                                                                           gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGG‑G‑CCAGGGGGGGTTTTTTTTTTACCAGGGAGAAGGGCAACAGCTGATTGCCCTTCACCCCCTCGCCCTGGGAGAGTTGGAGCAAGCGGGCCCCGCTGGgtt                                                            >  1:175546/1‑150 (MQ=255)
                                                                                                    aaTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagcag                                                  <  1:662946/151‑1 (MQ=255)
                                                                                                      tCGGCCAACGCGCGGGGAGAGGCGGGTTTCGTA‑TTTGGC‑G‑GCCGGGGGGGTTTTTTTTTCACCCGTGAGACGGGGAACAACTGATTGCCCCTTACCCCCTGGCCCTGAGAGAGGTGCAGCAAACGGTCCCCGCCGGGTTGCCCCCGCAgg                                                 >  1:757119/1‑150 (MQ=255)
                                                                                                       cGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGc                                                <  2:424397/150‑1 (MQ=255)
                                                                                                        ggCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCg                                               <  1:379525/150‑1 (MQ=255)
                                                                                                          ccAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaaa                                            <  2:1023107/151‑1 (MQ=255)
                                                                                                           cAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaaaa                                           <  2:956005/151‑1 (MQ=255)
                                                                                                            aaCGCGCGGGGAGAGGCGGTTTGTGTA‑TTTGGC‑G‑CCCGGGGGGGTTTTCTTTTTACCAGTGAGACGGGGAAAAAGTGGTTGCCCCTCACCCCCTGGGCCCGAGAGAGGTGGAGCAAGCGGGCCACACTGGCTTGCCCCCGCCGGCGaaaa                                           >  2:243372/1‑150 (MQ=255)
                                                                                                                cgcgGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCt                                       <  2:114/150‑1 (MQ=255)
                                                                                                                      gagaGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAt                                <  1:38793/151‑1 (MQ=255)
                                                                                                                      gagaGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGGGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTg                                  >  1:372817/1‑149 (MQ=255)
                                                                                                                          ggCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtg                            <  2:468248/151‑1 (MQ=255)
                                                                                                                          ggCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtg                            <  1:267971/151‑1 (MQ=255)
                                                                                                                                ttGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTaa                       <  1:715986/150‑1 (MQ=255)
                                                                                                                                  gCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg                     <  1:187349/150‑1 (MQ=255)
                                                                                                                                      a‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg                >  1:933898/1‑151 (MQ=255)
                                                                                                                                        ttGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg                >  1:960210/1‑150 (MQ=255)
                                                                                                                                           ggC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat            <  1:735984/151‑1 (MQ=255)
                                                                                                                                             c‑g‑cCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa          <  2:1041185/151‑1 (MQ=255)
                                                                                                                                                  cAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaaaa        >  2:848776/1‑148 (MQ=255)
                                                                                                                                                  cAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACAt       >  1:343676/1‑151 (MQ=255)
                                                                                                                                                    gggTGTTTTTTCTTTTCCCCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGTTCCAGCCGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGa     <  2:440636/151‑1 (MQ=255)
                                                                                                                                                    gggTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAAAATg      >  1:356637/1‑150 (MQ=255)
                                                                                                                                                       tggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCt  >  1:360626/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
AAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCT  >  NC_000913/366378‑366676

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