Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,538,1090AC41.8% 214.3 / 42.1 123G219G (GGT→GGGnanMN‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (57/14);  new base C (0/51);  total (57/66)
Fisher's exact test for biased strand distribution p-value = 8.61e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTG  >  NC_000913/4537960‑4538253
                                                                                                                                                              |                                                                                                                                                   
aTCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCa                                                                                                                                                     <  2:1294725/150‑1 (MQ=255)
  ccccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCAcc                                                                                                                                                   <  2:856140/150‑1 (MQ=255)
    ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCttt                                                                                                                                               >  2:190202/1‑151 (MQ=255)
    ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCttt                                                                                                                                               <  2:932944/151‑1 (MQ=255)
       cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTAtt                                                                                                                                            >  2:1376272/1‑151 (MQ=255)
        aaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTAtt                                                                                                                                            <  2:1264927/150‑1 (MQ=255)
        aaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTc                                                                                                                                          >  1:480073/1‑151 (MQ=255)
        aaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTc                                                                                                                                          >  1:1550361/1‑151 (MQ=255)
        aaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTc                                                                                                                                          >  2:1638720/1‑151 (MQ=255)
            ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTcaca                                                                                                                                      >  2:87509/1‑150 (MQ=255)
            ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTcaca                                                                                                                                      >  1:127417/1‑150 (MQ=255)
            ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACaa                                                                                                                                     <  2:388858/151‑1 (MQ=255)
             cGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAAc                                                                                                                                    >  1:548085/1‑151 (MQ=255)
              gCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAAcc                                                                                                                                   >  2:431914/1‑151 (MQ=255)
               ccAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAAccg                                                                                                                                  >  2:1596103/1‑151 (MQ=255)
                 aGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAAccgcc                                                                                                                                <  1:213443/151‑1 (MQ=255)
                 aGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAAccgcc                                                                                                                                <  2:1105251/151‑1 (MQ=255)
                   cTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACcgccgc                                                                                                                              >  2:204493/1‑151 (MQ=255)
                   cTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACcgccgc                                                                                                                              >  2:934201/1‑151 (MQ=255)
                   cTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTAAAATCAAGTTCAAATACGGCATCCGTTCGCAATCTGGTTTTGTCTTCGCCATTAATAAGCCAGTTTTTATCACCTTTATTCACAAccgccg                                                                                                                               <  1:334255/150‑1 (MQ=255)
                    tACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACcgccgc                                                                                                                              >  2:1482394/1‑150 (MQ=255)
                           tCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCg                                                                                                                      >  2:141813/1‑151 (MQ=255)
                            cTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGt                                                                                                                     >  2:227722/1‑151 (MQ=255)
                               gTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTcc                                                                                                                  >  2:314019/1‑151 (MQ=255)
                                tgatgaGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCg                                                                                                                 >  1:1204588/1‑151 (MQ=255)
                                    gaggCGGGAGCAAGTTTTTTCCTTTTAAAATTATTCCCGGTAAAATCAAGTTAAAATCGGGCATCCTTCCGCAATCCTGTTTTGGCTTCGCAATAAAAAACCCAGTTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTAcc                                                                                                             <  2:665440/151‑1 (MQ=255)
                                    gagaCGGGAGCAATTTTTTTCCATTTAAATTTATTCCCGGTAAAAAAAAGTTCAAATCGGGCATCCTTTCGCATCCCTGTTTTGGTCTCGCCATTAAAAAGCCAGTTTTTACCACTTTTATCCCAAACCGCCGCACCAGCCGTTCCGTAcc                                                                                                             <  2:550029/151‑1 (MQ=255)
                                                aGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc                                                                                                 >  2:1572350/1‑151 (MQ=255)
                                                 agttttttcATTTTAAATTTATTCCCGGAAAAATAAATTAAAAATACGCCACCCTTCCCCACCCCTGTTTTGGCTTCCCAATAAAAAAGCCGGTTTTTACCCCTTTTATCAACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc                                                                                                 <  1:255663/143‑1 (MQ=255)
                                                   ttATTCCTTTTTAATTTTTTCCCGGTGAAATAAATTTCAAATTCGCCCTCCTTCCCAAACCCTGTTTGGGCTCCGCCTTAAATACGCCGGTTTTTATCTCCTTTTTTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCa                                                                                               <  2:188145/150‑1 (MQ=255)
                                                      ttCCTTTTAAATTTATTACCGGAGAAAAAAATTAAAAATACGGCACCCGTGCGAAACCCGGGTTGGGCTTCGCCATAAAAAAGCGAGTTTTTTTCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGc                                                                                             <  2:1418661/149‑1 (MQ=255)
                                                      ttCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt                                                                                           >  1:1480372/1‑151 (MQ=255)
                                                             ttAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                     >  1:1352681/1‑150 (MQ=255)
                                                               aaaTTTTTACCGGAGAACCAATGTCAAAAACCTGCCTCCGTACGCAACCGGGGTTTGGGTCACCAATAAAAAACCTAGTTTTTATCCCTTTTTTTCCCATCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                  <  2:900399/151‑1 (MQ=255)
                                                                 attattACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGc                                                                                >  2:655552/1‑151 (MQ=255)
                                                                   ttattcCGGGTGAAATAAAGTTAAAATCCGCCTCCCGTTCGCAACCTTGTTTGGTCTCCCCAATAAAAAAGCCGGTTTTTATCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGc                                                                                <  2:1204681/150‑1 (MQ=255)
                                                                        aCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt                                                                          >  1:29374/1‑151 (MQ=255)
                                                                           ggtaaaataaattaaaaaTCCGCCACCCTTCCGCAATCCTGTTTTGGCTTCGCAATTAATAAGCCAGTTTTTATCCCTTTTTTTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAg                                                                       <  2:220894/136‑1 (MQ=255)
                                                                           ggagAATTCAAGTTAAAATCGGGCACCCTTTCGCAATCGTTTTTGGGCTCCACAATTAATAACAAAGTTTTTATCCCTTTTATACAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAg                                                                       <  1:679930/148‑1 (MQ=255)
                                                                               aatccaattcaaaaTCCGCCTTCCGTTCCCACTCTTGTTTTGGCTCCCCCATAAAAAGCCCAGGTTTTATCCCTTTTATCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGa                                                                    <  1:1469315/139‑1 (MQ=255)
                                                                                    aagtaaaaaTACGCCTTCCTTTCCTATCCCTTGTTGGCCTTCTCTATAAATAGCCAAGTTTTTACCCCTTTTATCCCCATCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCa                                                                 <  1:403933/142‑1 (MQ=255)
                                                                                    aaGTTAAAATCCGGCCCCCTTTCCCATTCCTGTTTTGGCTTCCCCATTAAAAAGCCAGTTTTTACCCCCTTTTTCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaaa                                                              <  1:1312496/151‑1 (MQ=255)
                                                                                     aattaaaactccGGCGTCCTTACACAATCTTGTTTTGCCTCCGCCATTAAAAAGCAATTTTTTACCCCTTTTATCAACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCa                                                                 <  1:1578274/137‑1 (MQ=255)
                                                                                     aattCAAAATCGGCCATCCGTTCGAATTCCTGTTTGGGCTTCCCAATTAATAGCCCAGTTTTTATCCCCTTTATCCTAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                <  1:1394447/147‑1 (MQ=255)
                                                                                       ttCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                           >  2:1210571/1‑151 (MQ=255)
                                                                                         cAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGaca                                                          >  1:1598421/1‑150 (MQ=255)
                                                                                          aaaTGCGGCCTCCGTCCGCAATCGTGGTTTGGTTCCGCCTTAAATAACCCGGTTTTTATCCCTTTTATCCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                        <  1:1552894/151‑1 (MQ=255)
                                                                                          aaaTCCGCCACCCGTTCGCAATCCTGTTTTGGCTCCCCCATTAATAAGCCGGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                           <  2:716348/148‑1 (MQ=255)
                                                                                          aaaTACGGCCCCCTTCCGCAACCCTGTTTTGGCTTCCCCTTTAAAAACCCAGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                        <  1:570599/151‑1 (MQ=255)
                                                                                            aTACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                        >  2:685502/1‑149 (MQ=255)
                                                                                             aaaaacgactCCGTTCCGCAATCTGGTTTTGCCTTCGCCATAAATAACCCAGTTTTTACCCCTTTTTTACAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                           <  1:655493/139‑1 (MQ=255)
                                                                                                 ggCATCCTTCCGCAATCTTGTTTTGGCTCCCCAATTAAAAAGCCAGTTTTTATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa                                                    <  1:92910/149‑1 (MQ=255)
                                                                                                 ggCATCCGTTCGAAATCCGGTTTTGCTTTCGCCATTAAAAAGCCGGTTTTTATCCCTTTTATCCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAAc                                                  <  2:693942/151‑1 (MQ=255)
                                                                                                     aTCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                >  1:81144/1‑150 (MQ=255)
                                                                                                      cctctgtccaaaTCCTGGTTTTGCCTCCCCCATAAAAAAGCAGGTTTTTCCCCCTTTTACCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa                                                    <  1:1258364/137‑1 (MQ=255)
                                                                                                        cGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                            >  2:890257/1‑151 (MQ=255)
                                                                                                        cGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                            >  1:627356/1‑151 (MQ=255)
                                                                                                          gtccGCAATCTTGGTTTGGCTTCGCCATTAATAACCCAGTTTTTACCCCTTTTATCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                <  1:201901/143‑1 (MQ=255)
                                                                                                          ctcgcacattcTGGTTTTGGTTCCCCATATAAAAAGCATGTTTTTTCCCCCTTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCACTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                <  2:1668794/137‑1 (MQ=255)
                                                                                                           tCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCTGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                           >  1:827730/1‑149 (MQ=255)
                                                                                                           tCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                           >  1:197758/1‑149 (MQ=255)
                                                                                                             caaaacccttttttGTCTTCCCCATTAAAAAGCCATTTTTTATCCCTTTTATCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                           <  1:90751/137‑1 (MQ=255)
                                                                                                                aaTCCTGGTTTGTCTCCCCAATAAAAACGCCAGGTTTTACCACCTTTTCCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaag                                         <  2:77870/147‑1 (MQ=255)
                                                                                                                   ccTGGTTTGGTTCCCCAATAAATAACCCAGTTTTTCCCCCCTTTTTCCCAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaata                                    <  1:575148/149‑1 (MQ=255)
                                                                                                                   ccTGGTTGGGTTGCGCAATTAATAACCCGGTTTTTATCCCTTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaa                                          <  1:890873/143‑1 (MQ=255)
                                                                                                                    cTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtc                                 >  2:1457882/1‑151 (MQ=255)
                                                                                                                     tGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtct                                >  2:640960/1‑151 (MQ=255)
                                                                                                                      gttttGGCTTCCCCATAAATAAGCCAGTTTTTATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataa                                   <  2:341419/145‑1 (MQ=255)
                                                                                                                      gttttGGCTCCCCAATAAATAACCCAGTTTTTTCCCCTTTAATCAACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACt                                                 <  1:1594701/131‑1 (MQ=255)
                                                                                                                      ggTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtctt                               >  2:544834/1‑151 (MQ=255)
                                                                                                                        tttGGTTCCGCTATTAAAAACCCAGTTTTTATCACCTTTATCAACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtct                                <  2:524144/148‑1 (MQ=255)
                                                                                                                        tttGCCTTCCCAATAAAAAGCCCAGTTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataa                                   <  1:1261745/145‑1 (MQ=255)
                                                                                                                        tttGCCTCCGCCATTAAAAGGCCAGTTTTTACCCCTTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtctt                               <  1:42866/149‑1 (MQ=255)
                                                                                                                         ttGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                            >  1:877716/1‑151 (MQ=255)
                                                                                                                              cTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttt                        >  2:891527/1‑151 (MQ=255)
                                                                                                                                cccgcATTAAATAGCCCGGGTTTTTCCCCTTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttt                        <  1:1376178/146‑1 (MQ=255)
                                                                                                                                  cccaataaaaaaGCCAGTTTTTACCCCTTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                            <  1:87505/134‑1 (MQ=255)
                                                                                                                                     aTTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                  >  1:1651172/1‑150 (MQ=255)
                                                                                                                                        taataaGCCAGGATTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa                >  2:515417/1‑150 (MQ=255)
                                                                                                                                         aataaGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa               >  1:990482/1‑150 (MQ=255)
                                                                                                                                         aataaGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa               >  2:1150563/1‑150 (MQ=255)
                                                                                                                                         aaaacgCCAGTTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                 <  1:1482295/143‑1 (MQ=255)
                                                                                                                                            aataaccctttttttATCCCTTTTATTCACATCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa                <  1:987688/138‑1 (MQ=255)
                                                                                                                                            aaaaggcgtgtttttCTCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttttt                      <  1:1638615/132‑1 (MQ=255)
                                                                                                                                            aaaacgccggtttttATCCCTTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaaa              <  1:163663/140‑1 (MQ=255)
                                                                                                                                            aaaaacccggtttttATCCCTTTTTTTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa                <  1:658859/138‑1 (MQ=255)
                                                                                                                                            aaGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt             >  1:93304/1‑149 (MQ=255)
                                                                                                                                              aaccAGGTTTTTATCACATTTATTAAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGGTTTTTTGTCAAAATa            <  1:1430116/148‑1 (MQ=255)
                                                                                                                                               cccAGTTTTTACCCCTTTTTTCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg         <  1:366829/150‑1 (MQ=255)
                                                                                                                                               accccggTTTTTCCCCCTTTATTCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt                       <  2:47926/134‑1 (MQ=255)
                                                                                                                                                 ccgggTTTAATCCCCTTTCTCCACATCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa                <  2:1647465/140‑1 (MQ=255)
                                                                                                                                                  ccggtttttACCCCCTTTATTCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGa        <  2:920048/146‑1 (MQ=255)
                                                                                                                                                     gTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAt    >  2:648472/1‑151 (MQ=255)
                                                                                                                                                     gTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAt    >  2:473559/1‑151 (MQ=255)
                                                                                                                                                       ttttATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAtt   >  2:409029/1‑151 (MQ=255)
                                                                                                                                                        tttATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg  >  2:5650/1‑151 (MQ=255)
                                                                                                                                                        tttATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg  >  2:254743/1‑151 (MQ=255)
                                                                                                                                                         ttATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg  <  2:501009/150‑1 (MQ=255)
                                                                                                                                                         ttATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg  >  1:1323041/1‑150 (MQ=255)
                                                                                                                                                         ttACCACTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg  <  2:554729/150‑1 (MQ=255)
                                                                                                                                                              |                                                                                                                                                   
ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTG  >  NC_000913/4537960‑4538253

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