Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,065,5800TA48.1% 60.1 / 53.4 83T63S (ACC→TCC) yihOputative sulphoquinovose importer
Reads supporting (aligned to +/- strand):  ref base T (26/16);  new base A (1/38);  total (27/56)
Fisher's exact test for biased strand distribution p-value = 3.01e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGA  >  NC_000913/4065433‑4065717
                                                                                                                                                         |                                                                                                                                            
aacaaaAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTaaa                                                                                                                                           >  2:1008304/1‑151 (MQ=255)
    aaaGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAAc                                                                                                                                      <  2:1197520/151‑1 (MQ=255)
    aaaGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAAc                                                                                                                                      <  2:544000/151‑1 (MQ=255)
      aGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACtt                                                                                                                                    >  2:495268/1‑151 (MQ=255)
              gCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGccacca                                                                                                                            >  1:1379742/1‑151 (MQ=255)
               cGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAg                                                                                                                           >  2:1044637/1‑151 (MQ=255)
                 aTAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGaa                                                                                                                         <  1:1165547/151‑1 (MQ=255)
                   aaaCTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAg                                                                                                                       >  1:223300/1‑151 (MQ=255)
                    aaCTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGa                                                                                                                     <  2:110947/151‑1 (MQ=255)
                    aaCTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGGGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGa                                                                                                                     >  1:1422773/1‑151 (MQ=255)
                      cTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGGCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCCGTAAACGCGGAAAAAAACTTCTCCACCAGAAAGATa                                                                                                                   >  2:262020/1‑151 (MQ=255)
                      cTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATa                                                                                                                   >  1:60218/1‑151 (MQ=255)
                       tGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATaa                                                                                                                  <  1:778904/151‑1 (MQ=255)
                        ggAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAAt                                                                                                                 >  2:1071162/1‑151 (MQ=255)
                             gTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCg                                                                                                            >  2:818387/1‑150 (MQ=255)
                                gCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCAt                                                                                                        >  2:706364/1‑151 (MQ=255)
                                gCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGGGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCAt                                                                                                        >  1:1237035/1‑151 (MQ=255)
                                           cccgggAACCGCTAAATATAAAAAAAAGGGAAAAATTTGGCTTTTAGCCCCAATTTTTTCCCGCGAGGCGAAAAAAACTCCGGGGAGCATTTCGGTAAACGCGAAAAAAACCTCCCCCACAAAAAAAATAATCCCGCCATAGTAAGCAGGc                                                                                             <  2:931398/147‑1 (MQ=255)
                                                gacccGCTGCATATAAAAAAAAGGAACGAAATTGGCCTTTAGACCAAATATTTTTCCGAGAGTCGAGTAAAAACCCGGTGAGCATATCGGAAAACGGGAAAAAAAACTCCGCCCCCAAAAAAAAAACCCCGCCATAGTAAGCAGGCATcc                                                                                         <  2:1564753/147‑1 (MQ=255)
                                                gAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATcc                                                                                         >  2:63448/1‑150 (MQ=255)
                                                 accctcTGAATAAAAAAAAAAGGAAAAAAATTGGCTTTTAGCCCAAATTTTTTCCCGGGAGGCGAAAAAAAACCCGGTGACCATATCGGAAACCGCGAAAAAAACCTCCCCCACCAAAAAAATAACCCCGCCATAGTAAGCAGGCATccc                                                                                        <  2:92539/145‑1 (MQ=255)
                                                     gCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAAtt                                                                                   >  1:50376/1‑151 (MQ=255)
                                                     cttgaataaAAAAAAAAGGGACAAATTTTGCTTTTAGCCCAATTTTTTTCCCGCGAGTCGAGTAAAAACCCGGGGAGCATATGGGTAAACGGGGAAAAAACCTCCGCCACAAGAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAAtt                                                                                   <  1:435045/146‑1 (MQ=255)
                                                      ttGCATAAAAAAAAAAGGGAAAAAATTTGCCTTTAGGCCAAATATTTTCCCGCGAGTCGAGAAAAAATCCGGTGAGCATATCGGAAAACGCGGCAAAAACCTCCGCCACCAGAAAGAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTc                                                                                  <  2:605177/150‑1 (MQ=255)
                                                      ctgaaaaaaaaataaaGGGACGAATTTTTCCTTGGGCCCAAATTTTTTTCCGCGAGTCGAGAAAAAATCCGGTGAGCTTTTCGGTAAACGCGGAAAAAAACTTCGCCCCCAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTc                                                                                  <  1:93564/143‑1 (MQ=255)
                                                      ctacaaaaaaaataaaGGGAAGAAATTTGCTTTTAGGCCAAATATTTTCCCGCGAGTGGAGAAAAAACCCGGGGAGAATATCGGAAAACGCGGAAAAAAACTCCGCCACAAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTc                                                                                  <  2:76233/143‑1 (MQ=255)
                                                      cTGAATAAAAAAAAAAGGGACAAATTTTCCTTTTAGGCCAAATTTTTTCCCGGGAGTCGAGAAAAAACCCGGTGAGCATATGGAAAAACGCGGAAAAAAACTTCGCCACCAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAAtt                                                                                   <  1:129831/150‑1 (MQ=255)
                                                       tGCATAAAAAAAAAAGGGACAAATTTTGCCTTTAGGCCAAATTTTTTCCCGCGAGTCGAGAAAAAACCGGGTAAAAATACCGAAAAACGCGGAAAAAACCTCCGCCCCAAAAAAGAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCa                                                                                 <  2:386948/151‑1 (MQ=255)
                                                       ctaaaaaaaaaaaaaAGGACAGAAATTAGCTTTAAGACCAAATATTTTCCCGCGGATAGAGTAAAAACCCGGGGCGCATTTGAATAACCGGGGAAAAAACCTTCACAACAAAAAAAAAACTCCCGCCATAGTAAGCAGGCATCCCTAATTc                                                                                  <  2:168590/144‑1 (MQ=255)
                                                          atatAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCg                                                                              >  2:262882/1‑151 (MQ=255)
                                                                aataaaGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATaa                                                                        >  2:871944/1‑151 (MQ=255)
                                                                   aaaGGGCCAAATTTGGCCTTTAGCCCCAATATTTTTCCGCGAGTGGAGAAAAAACCCGGTGAGCATATGGGTAAACGGGGAAAAAAACTTCGCCCCCAGAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAt                                                                     <  1:555997/151‑1 (MQ=255)
                                                                   aaaGGGACGAAATTTGCTTTAGGCCCAAATTTTTTCCCGGGAGTCGAGAAAAAACCCGGTGACCATATCGGAAACCGCGGAAAAAACCTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAt                                                                     <  2:1473726/151‑1 (MQ=255)
                                                                   aaaGGGACAAAATTTCCCTTTAGGCCAATTTTTTTCCCGGGAGGCGAGAAAAAATCCGGTGAGCATATGGGTAAACGCGGTAAAAACCTTCCCCACAAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAt                                                                     <  1:1212558/151‑1 (MQ=255)
                                                                      gggACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTg                                                                  >  2:1239053/1‑151 (MQ=255)
                                                                       ggACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTga                                                                 >  2:106233/1‑151 (MQ=255)
                                                                                tttccttttAGGCCAATTTTTTTCCGGGGTGTGGAAAAAAAACCCGGGGAGAATATGGAAAACCGGGAAAAAAAATTTCGCCACAAGAAAGAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATAc                                                         <  1:1580694/145‑1 (MQ=255)
                                                                                tttGCTTTTAGGCCAAATATTTTTCCGCGAGTTGATTAAAACTCCGGTGAGCATTTCGATAAACGCGGTAAAAAACTTCCCCCCCAAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATAc                                                         <  1:819293/151‑1 (MQ=255)
                                                                                  tGCTTTGAGGCCCATTTTTTTCGAGGGAGGCAAGAAAAAGCCCGGGGAGAATTTCGAAAAACGCGGTAAAAACCTCCACCACAAGAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACag                                                       <  2:212016/151‑1 (MQ=255)
                                                                                  tGCGTTTAGGCCAAATTTTTTCCCGCGAGTCGAGTAAAAACCGGGGGAGAATGCCGGAAAACGCGGAAAAAACCTCCCCCACCAGAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACag                                                       <  2:1486295/151‑1 (MQ=255)
                                                                                  tGCCTTTAGGCCAATTTTTTTTCCGCGAGTGGAATAAAAATCGGGTGGGCATACCGGTAAACGCGGAAAAAAACTCCGCCCCAAAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACag                                                       <  2:400412/151‑1 (MQ=255)
                                                                                    cttttGGGCCAAATTTTTTTCCGCGAGTCGAGAAAAACCCGGGTGAGAATTTCGGAAAACGCGAAAAAAACCTCCCCCCCAAAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACagag                                                     <  2:171939/149‑1 (MQ=255)
                                                                                    cttttAGCCCAATTATTTTTCCGAGAGTCGAGAAAAAACCCGGTGAGATTATCGGAAAACGCGGAAAAAAACTCCGCCACAAAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACagag                                                     <  2:224815/149‑1 (MQ=255)
                                                                                     ccttggggcaaaaTATTTTGCCGCGGGTTAAGAAAAAACCGGGGGACAATTTGGCAAAGCGGGATAAAAAATTTCGCCACGAAAAAAAAAATCCCGCCATAGTGAGCAGGCATACCTAATTCATCGGTATAAAATTTGAGGAGATACagag                                                     <  1:1250069/140‑1 (MQ=255)
                                                                                     cTTTTGGCCCAAATTTTTTCCCGGGATTCGAGAAAAACCCCGGTGAGCATATCGGAAACCGCGGAAAAAAACTCCCCCCCAAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACagag                                                     <  1:344447/151‑1 (MQ=255)
                                                                                     cTTTAGGCCAAATTTTTTTCCGGGAGGGGAGAAAAAATCGGGGGAGCATATGGGAAACCGGGGAAAAAAACTCCCCCCCAAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGt                                                    <  1:924456/151‑1 (MQ=255)
                                                                                       ttAGGCCAAATATTTTTCCGCGAGTCGAGAAAAAACCCGGGGAGCATTTCGGTAAACGCGGTAAAAAACTTCCCCCCCAAAAAGAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTAc                                                  <  1:127178/151‑1 (MQ=255)
                                                                                        tAGCCCAAATATTTTTCCGCGAGTCAAGAAAAAACCGGGTGAGCATATCGGAAAGCGCGAAAAAAACCTCCCCCACCAAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTAcc                                                 <  2:905326/151‑1 (MQ=255)
                                                                                           gCCAAATTTTTTCCCGGGAGGCGAGAAAAAACCGGGTGAGAATTGCGGAAAACGGGGAAAAAAACCCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGa                                               <  1:262949/150‑1 (MQ=255)
                                                                                           gCCAAATTTTTTCCCGGGAGCCGAATAAAAACCCGGTGAGCATATGGGAAAACGGGAAAAAAAACTCCCCCCCCAAAAAGAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGa                                               <  1:706735/150‑1 (MQ=255)
                                                                                            cccAATTATTTTCCCGCGAGTAGAGAAAAAACCCGGTGAGCATATCGGTAAACGCGAAAAAAAACTCCCCCACCAAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGa                                               <  2:76378/150‑1 (MQ=255)
                                                                                            cccAATTTTTTCCCGGGAGTGGAGAAAAAACCCGGTGAGAATATCGGAAACCGCGGAAAAAAACTTCGCCCCCAAAAAGAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATg                                             <  1:955718/151‑1 (MQ=255)
                                                                                             ccaaaTTTTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATACCGGTAAACGCGGTAAAAAACTCCGCCACAAGAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATg                                             <  2:144485/148‑1 (MQ=255)
                                                                                              caaaTATTTTCCCGCGAGTCGAGTAAAAAGCCGGTGAGCATATCGGTAAACGCGTAAAAAAACTTCGCCACCAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGc                                            <  2:892487/149‑1 (MQ=255)
                                                                                              cAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCa                                           >  2:1307964/1‑151 (MQ=255)
                                                                                                 atatTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga                                         >  2:1622433/1‑151 (MQ=255)
                                                                                                 atatTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga                                         >  1:131102/1‑151 (MQ=255)
                                                                                                    tttttcCCGGGGGTCGAGAAAAAACCCGGTGGACATATCGGAAAACGCGGAAAAAACCTCCCCCACAAAAAAGAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGc                                       <  2:1286100/150‑1 (MQ=255)
                                                                                                             cgggattcgagaAAAAATCGGGTGATCATATCGGTAAACGGGGAAAAAACCTCCGCCACAAAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAAt                                <  2:248825/145‑1 (MQ=255)
                                                                                                             cgcgatccaagaaaaaaCCCGGTGACCTTATCGTAAACCGCGGTAAAAACCTCCCCCCCCAAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAAt                                <  2:282236/139‑1 (MQ=255)
                                                                                                             cgagaggcggaaaaaaaCCCGGAAAACATATCGGTAAACGGGAAAAAAAACTCCGCCCCAAAAAGAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAAt                                <  1:812168/139‑1 (MQ=255)
                                                                                                             cacaagtcgagtAAAACACCGGTGAAAATATCGGAAACGTCGGTAAAAACCTCCCCCCCCAAAAATAAATCTCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAAt                                <  1:642668/147‑1 (MQ=255)
                                                                                                               cgtgtaaagaaaaaaCCCGGGGAGATTCCCGGTAAACGTGGAAAAAAACCTCCCCACCAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAAtt                               <  1:427091/140‑1 (MQ=255)
                                                                                                               cgagtcgagaAAAAACCGGGTGAGCATATCGGAAAACGCGGAAAAAACCTCCGCCACCAGAAGGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAAtt                               <  2:187778/150‑1 (MQ=255)
                                                                                                                     gagtAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCg                        >  1:101627/1‑151 (MQ=255)
                                                                                                                             acccGGTGAGCATATCGGTAAACGCGATAAAAAACTCCGCCCCCAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGcc                <  1:1054775/149‑1 (MQ=255)
                                                                                                                              ccggggnngCANACCNNNNNNNGCNAAAANNNNNNNCCCCTCAAGAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGccc               <  2:1151088/143‑1 (MQ=255)
                                                                                                                              cccGGTGATCATAGCGTAAAACGCGAAAAAAAACTACCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATCCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGccc               <  2:1552445/150‑1 (MQ=255)
                                                                                                                                   tagacATACCGGAAAACGGGGTAAAAAACTTCCCCACAAGAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCccatacc          <  2:1543798/147‑1 (MQ=255)
                                                                                                                                   tGAGCATATCGGTAAACGCGGAAAAAACCTTCCCCCCCAAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCccatacc          <  1:201800/151‑1 (MQ=255)
                                                                                                                                   tGAGCATATCGGTAAACGCGGAAAAAAACTCCGCCACAAAAAAGAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCccatac           <  2:1127708/150‑1 (MQ=255)
                                                                                                                                   tGAGCATATCGGAAAACGGGGAAAAAACCTTCGCCACAAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCccatac           <  1:117664/150‑1 (MQ=255)
                                                                                                                                   tGAGCATATCGGAAAACGGGGAAAAAACCTCCGCCCCAAAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCccatacc          <  1:825974/151‑1 (MQ=255)
                                                                                                                                   tGAGAATAGCGTAAAACACGAAAAAACACTTGGCCACAAAAAAAATAGTCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCccatac           <  2:687336/150‑1 (MQ=255)
                                                                                                                                    gagaATATCGGTAACCGCGGAAAAAAACTCCGCCACAAGAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCccatacc          <  2:1456018/150‑1 (MQ=255)
                                                                                                                                    gAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCccatacca         >  2:1291900/1‑151 (MQ=255)
                                                                                                                                    gAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCccatacca         >  2:760657/1‑151 (MQ=255)
                                                                                                                                     agaatatCGGAAAACGCGGAAAAAAACTCCGCCACAAGAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCccatacca         <  1:241117/147‑1 (MQ=255)
                                                                                                                                     aGCATATCGGAAAACGCGATAAAAACCTTCCCCCCAAGAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCccatacca         <  2:223244/150‑1 (MQ=255)
                                                                                                                                      gCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCccataccata       >  1:210029/1‑151 (MQ=255)
                                                                                                                                      gCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCccataccata       >  1:113544/1‑151 (MQ=255)
                                                                                                                                      acaatatgggaaaCGGGGGAAAAAACCTTCCCCACAAAAAAGAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCccatacca         <  1:1239631/140‑1 (MQ=255)
                                                                                                                                           tCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGa  >  1:158269/1‑151 (MQ=255)
                                                                                                                                                         |                                                                                                                                            
AACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGA  >  NC_000913/4065433‑4065717

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