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

AAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCT  >  NC_000913/366378‑366676
                                                                                                                                                       |                                                                                                                                                      
aaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGATTTCGGATTTGGCGCCAgggggg                                                                                                                                                      >  1:922358/1‑151 (MQ=255)
  gCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGATTTGGGGCCCgggggg                                                                                                                                                      >  2:38792/1‑149 (MQ=255)
       gggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttct                                                                                                                                               <  2:1015882/151‑1 (MQ=255)
         gTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttcttt                                                                                                                                             <  2:541004/151‑1 (MQ=255)
          tGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTTTTGGGCGCCCGGGGGGtttttttttt                                                                                                                                            >  2:379529/1‑151 (MQ=255)
                    gtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTgag                                                                                                                                  <  1:54001/151‑1 (MQ=255)
                     tgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGACAGGCGGTTTTTGCATTTGGCGGCCCGGGGGGtttttttt                                                                                                                                              >  2:697109/1‑139 (MQ=255)
                        gCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACgg                                                                                                                              <  1:409078/151‑1 (MQ=255)
                        gCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACg                                                                                                                               <  2:655972/150‑1 (MQ=255)
                         cTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCTGGTTGTGAGACCGGGGGGGtttttttttttccaagtggaaggg                                                                                                                              >  2:203772/1‑137 (MQ=255)
                          tAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACggg                                                                                                                             >  1:61465/1‑150 (MQ=255)
                           aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGGATTTGGGGGCAGGGGGGTTTTTTTTTTTAACAGTggggaggg                                                                                                                             >  2:716004/1‑142 (MQ=255)
                           aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACggg                                                                                                                             <  1:423530/149‑1 (MQ=255)
                            aCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCaa                                                                                                                          <  2:774735/151‑1 (MQ=255)
                             cTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGGGTTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTgagaacgggaa                                                                                                                          >  2:494618/1‑143 (MQ=255)
                               cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACa                                                                                                                        <  1:1051494/150‑1 (MQ=255)
                                acaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGc                                                                                                                      <  1:860252/151‑1 (MQ=255)
                                  aTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTg                                                                                                                    <  1:685944/151‑1 (MQ=255)
                                     aaTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGAtt                                                                                                                 <  1:610693/151‑1 (MQ=255)
                                      aTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTg                                                                                                                <  2:133651/151‑1 (MQ=255)
                                       ttgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGATTTGGGGCCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGAtt                                                                                                                 >  1:424400/1‑149 (MQ=255)
                                        tgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGc                                                                                                               <  1:229529/150‑1 (MQ=255)
                                          cgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCt                                                                                                            <  2:793010/151‑1 (MQ=255)
                                               cgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGCCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAAAAGCTGATTGCCCttttcc                                                                                                       >  2:913622/1‑147 (MQ=255)
                                                    aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTGTTTGGGGCCCGGGGGGGTTTTCTTTTTTCCAGTGAGACGGGGAACAGCCGGTTGTCCTTCACcccc                                                                                                    >  2:432463/1‑149 (MQ=255)
                                                    aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTg                                                                                                  <  2:971660/151‑1 (MQ=255)
                                                             tttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTTCGGGTTTGGCGCCGGGGGGGGTTTTTTTTTTACCCGGGAGACGGGCGACACCTGGATGCCCCTCACCGCCCGGgcc                                                                                              >  2:373708/1‑146 (MQ=255)
                                                               tCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCCGGGGGGGTTTTTTTTTCACCCGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgag                                                                                          >  1:288967/1‑148 (MQ=255)
                                                                ccAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTTGTGTTTGGGGCCCGGGGGGGTTTTTTTTTTACCCGGGAGACGGGGCACAGGCGATTTCCCCTTACCCCCCGGCCCCtagag                                                                                        >  1:955970/1‑149 (MQ=255)
                                                                ccAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgagaga                                                                                       <  2:891824/150‑1 (MQ=255)
                                                                ccAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTGCGGATTTGGCGGCCGGGGGGGTTTTTTTTTTACCCGGGAGACGGGCAACAGCTGATTGCCCTTCACCCCCCGGCCCTgagag                                                                                        >  1:461178/1‑149 (MQ=255)
                                                                  aGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTTGCATTTGGGGGCGGGGGGGGTTTTTTTTTTAACCTTGAGACGAGTAACAGGCGATTGCCCCTCACCGCCTGGCGCTGcgtgaggt                                                                                    >  1:679192/1‑149 (MQ=255)
                                                                  aGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTTCGGATTGGGGGGCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGAGTGTCCTTCACCCCCCGGCCCCGagagaggt                                                                                    >  2:408045/1‑149 (MQ=255)
                                                                   gTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTg                                                                                   <  2:1107003/151‑1 (MQ=255)
                                                                     cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTGGGGGCCCGGGGGGGTTTTCTTTTTCCCCGTGAGACGGGCAACACCTGATTGCCCCTCACCGCCTGGCCCCGagggaggt                                                                                    >  1:727327/1‑146 (MQ=255)
                                                                     cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTGCGTATTGGGGGCCGGGGGGGGTTTTTTTTTTCCCAGGGCGACGGGGAACAGCTGATTGGCCTTCACCGCCTGGTCCCGTGAGAGtt                                                                                    >  1:292328/1‑148 (MQ=255)
                                                                      gggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGGAACAGCTGATTGTCCCTTACCGCCTGGCCCCGggagag                                                                                      >  1:507328/1‑145 (MQ=255)
                                                                        gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTGGGGGGCAGGGGGGGTTTTTTTTTTACCAGGGAGAGGGGGAACAGCTGATTGCCCTTCACCCCCTGGCCCCGTGAGAGTTtcag                                                                                >  2:438078/1‑149 (MQ=255)
                                                                         aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCGCCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCaa                                                                             <  2:711976/151‑1 (MQ=255)
                                                                             cTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTGGGCGCCCGGGGGGGTTTTCTTTTTCCCAGGGAGACGGGCAACAGCTGATTGCCCTTTCCCCCCCCGCCCCGCGGGAGGTGCAGCAAACgg                                                                         >  2:513359/1‑151 (MQ=255)
                                                                                  gTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCa                                                                     <  1:367686/150‑1 (MQ=255)
                                                                                   tGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCAGGGGGGGTTTTCTTTTTACCAGTGAGACGGGCAAAAGCTGATTGCCCTTCACCCCCTGGCCCTGAGAGAGGTGCAGCAAGCGGTCCAc                                                                    >  2:328351/1‑150 (MQ=255)
                                                                                      cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGTGTAATTGGCGCCAGGGGGGGTTTTCTTTTTACCAGGGAGACGGGCCACAGCTGATTGCCCTTCCCCGCCCGGCCCTCAGGGAGGTTCAGCAAACGGGCCACGc                                                                  >  2:491407/1‑149 (MQ=255)
                                                                                       aGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCTTTTTTGGGGCCCAGGGGGGTTTTCTTTTTTCCAGTGAGAGGGGCCAAAGCCGATTTCCCCTCACCCGCTAGCCCTTTGAGAGGTGGAGCAAACGGTccccaccgg                                                               >  1:428308/1‑144 (MQ=255)
                                                                                          tGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCCTCACCCCCTGGCCCTGAGGGAGGTGCAGCAAACGGGCCCCCCCGGttt                                                            >  1:283552/1‑151 (MQ=255)
                                                                                          tGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGATTGCGGATTGGGGGCCCGGGGGGGTTTTCTTTTCCCCCGGGGGACGGGGAACAGCTGATTTCCCTTCACCGCCCGGCCCTTGGAGAGTTGCGGCAAGCGGGGCCCGCTGGgtt                                                            >  1:883949/1‑151 (MQ=255)
                                                                                           gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGCCAGGGGGGGTTTTTTTTTTACCAGGGAGAAGGGCAACAGCTGATTGCCCTTCACCCCCTCGCCCTGGGAGAGTTGGAGCAAGCGGGCCCCGCTGGgtt                                                            >  1:175546/1‑150 (MQ=255)
                                                                                                    aaTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagcag                                                  <  1:662946/151‑1 (MQ=255)
                                                                                                      tCGGCCAACGCGCGGGGAGAGGCGGGTTTCGTATTTGGCGGCCGGGGGGGTTTTTTTTTCACCCGTGAGACGGGGAACAACTGATTGCCCCTTACCCCCTGGCCCTGAGAGAGGTGCAGCAAACGGTCCCCGCCGGGTTGCCCCCGCAgg                                                 >  1:757119/1‑150 (MQ=255)
                                                                                                       cGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGc                                                <  2:424397/150‑1 (MQ=255)
                                                                                                        ggCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCg                                               <  1:379525/150‑1 (MQ=255)
                                                                                                          ccAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaaa                                            <  2:1023107/151‑1 (MQ=255)
                                                                                                           cAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaaaa                                           <  2:956005/151‑1 (MQ=255)
                                                                                                            aaCGCGCGGGGAGAGGCGGTTTGTGTATTTGGCGCCCGGGGGGGTTTTCTTTTTACCAGTGAGACGGGGAAAAAGTGGTTGCCCCTCACCCCCTGGGCCCGAGAGAGGTGGAGCAAGCGGGCCACACTGGCTTGCCCCCGCCGGCGaaaa                                           >  2:243372/1‑150 (MQ=255)
                                                                                                                cgcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCt                                       <  2:114/150‑1 (MQ=255)
                                                                                                                      gagaGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAt                                <  1:38793/151‑1 (MQ=255)
                                                                                                                      gagaGGCGGTTTGCGTATTGGGCGCCAGGGGGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTg                                  >  1:372817/1‑149 (MQ=255)
                                                                                                                          ggCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtg                            <  2:468248/151‑1 (MQ=255)
                                                                                                                          ggCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtg                            <  1:267971/151‑1 (MQ=255)
                                                                                                                                ttGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTaa                       <  1:715986/150‑1 (MQ=255)
                                                                                                                                  gCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg                     <  1:187349/150‑1 (MQ=255)
                                                                                                                                      aTTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg                >  1:933898/1‑151 (MQ=255)
                                                                                                                                        ttGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg                >  1:960210/1‑150 (MQ=255)
                                                                                                                                           ggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat            <  1:735984/151‑1 (MQ=255)
                                                                                                                                             cgcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa          <  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)
                                                                                                                                                       |                                                                                                                                                      
AAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCT  >  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: