Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,538,1090AC43.1% 115.9 / 24.6 72G219G (GGT→GGGnanMN‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (34/7);  new base C (0/31);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 6.46e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537964‑4538258
                                                                                                                                                      |                                                                                                                                                      
ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCttt                                                                                                                                                   >  2:418509/1‑151 (MQ=255)
ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCttt                                                                                                                                                   <  1:392161/151‑1 (MQ=255)
  gCAAAACCGCCAGCTACGCCATCTGGTGATGAGCCGGGAGCAAGCTTATTCCATTTTAAATTATTCCCGGTGAAATCAAGTTCAAATACCGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTa                                                                                                                                                  <  1:283478/150‑1 (MQ=255)
      aaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTcac                                                                                                                                            >  1:450093/1‑151 (MQ=255)
            cAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAAccgc                                                                                                                                      <  2:144450/151‑1 (MQ=255)
                tACGCCATCTGGCGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACcgccgc                                                                                                                                   >  1:690481/1‑150 (MQ=255)
                 aCGCCATCGGGTGATGAGACGGGAGCAAGCTTATTCCTTTTAAAATTATTCCCGGTAAAATCAAGTTAAAATCCGGCATCCTTTCGCAACCCTGGTTGGGCTTCGCCATTAATAAGCCGGGTTTTACCCCCTTAATTCACAACcgccgc                                                                                                                                   <  1:795061/149‑1 (MQ=255)
                     cATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGc                                                                                                                             >  2:450414/1‑151 (MQ=255)
                     cATCTGGTGAGGAGAGGGGAGCAAGTTTTTTCCTTTTTAAATTATTACCGGTAAAATCAAGTCAAAATCGGCCATCCGTCCGCAATCCTGTTTTGGCTTCGCCATAAATAAGCCGGGTTTTATCACCTTTATTCACAACCGCCGCACCAGc                                                                                                                             <  1:850797/151‑1 (MQ=255)
                        cTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGt                                                                                                                          >  1:77672/1‑151 (MQ=255)
                                  gaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTAcc                                                                                                                  >  1:672866/1‑149 (MQ=255)
                                     gggAGCAATCTTTTTCCATTTAAAATTATTCCGGGTAAAATCAATTTAAAATACGGCCTCCGTCCCCATCCCGGGTTGGCCTCCGCCATAAAAAACCCAGTTTTTACCCCTTTTATACACAACCGCCGCACCAGCCGTTCCGTACCAGGGc                                                                                                             <  1:184581/151‑1 (MQ=255)
                                     gggAGCAAGTTTTTTCCTTTTTAAATTATTCCCGGTAAAATCAAGTTCAAATACGCCATCCTTCCGCAATCCTGTTTTGGTTTCCCCATAAATAAGCCAGTTTTTATCCCTTTTATTCAAAACCGCCGCACCAGCCGTTCCGTACCAg                                                                                                                <  2:153225/148‑1 (MQ=255)
                                     gggAGCAACCTTTTCCTTTTTAAAATTATTCCGGGAAAAACAAATTCAAAATCCGGCACCCTTTCACATCCTGGTTTGGGCTTCCCCATTAAAACGCCAGTTTTTACCCCTTTTATTCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGc                                                                                                             <  1:460567/151‑1 (MQ=255)
                                      ggAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                            >  1:1034495/1‑151 (MQ=255)
                                           aaGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg                                                                                                       >  1:442717/1‑151 (MQ=255)
                                           aaGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg                                                                                                       >  2:473945/1‑151 (MQ=255)
                                           aaGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg                                                                                                       >  1:825965/1‑151 (MQ=255)
                                            aGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc                                                                                                      >  2:340503/1‑151 (MQ=255)
                                             gCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                     >  2:1055763/1‑151 (MQ=255)
                                             agtttttCCAATTTAAATTTTTACCGGTGAAAATAAATTAAAAAGCCGCCCCCCTTCCGAAACCTGGTTTGGGCTTCCCCTTAAATAACCAGGTTTTAACCCCTTTTTTCCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc                                                                                                      <  1:824862/146‑1 (MQ=255)
                                                        tttAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                          >  1:285024/1‑151 (MQ=255)
                                                         ttAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                          >  2:242383/1‑150 (MQ=255)
                                                          taaaaTTATTCCCGGTAAAAACAAGTAAAAATAGGGCCCCCGTCCCAATTCCTGTTTTGTCTCCGCAATAAAAACCCCGGTTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                          <  2:606314/148‑1 (MQ=255)
                                                            aaaTTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                       >  2:211806/1‑151 (MQ=255)
                                                               ttattaCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCt                                                                                    >  1:923652/1‑151 (MQ=255)
                                                                 attaCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgt                                                                                  >  1:1029614/1‑151 (MQ=255)
                                                                  ttcccGGTAAAATCAATTTCAAATACGGCATCCGTTCGCAATCCTGTTTTGCCTTCGCCATTAATAAGCCAGGTTTTACCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgt                                                                                  <  2:813790/147‑1 (MQ=255)
                                                                     cgggTAAAATCAAGTACAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCCCCATTAATAAGCCGGTTTTTATCCCTTTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgtg                                                                                 <  2:676068/146‑1 (MQ=255)
                                                                     ccGGTAAAATCAATTTAAAATCCGGCACCCGTTCGCAACCCTGTTTTGGTTTCGCCATTAAAAAGCCAGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                              <  2:843130/151‑1 (MQ=255)
                                                                        ggtaaaacaaatttaaaaTCCGGCCCCCGTTCGCAATCCTGGTTGGGTTTCGCCTTTAATAAGCCAGTTTTTATCCCCTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGa                                                                             <  1:279757/135‑1 (MQ=255)
                                                                        ggtaaaaCAAATTCAAATTCCGCCATCCTTTCGCAATCCTGTTTTGTTTTCCCCATAAAAACGCGGGTTTTTATCCCCTTTTTTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGa                                                                             <  2:349200/147‑1 (MQ=255)
                                                                         tGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                          >  1:573692/1‑151 (MQ=255)
                                                                            attCAATTTAAAATCCGGCATCCTTTCGCAATCCTGGTTTGGCTCCCCCATTAATAACCCAGTTTTTACCCCCTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc                                                                       <  1:750738/149‑1 (MQ=255)
                                                                               cAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                     >  1:387423/1‑150 (MQ=255)
                                                                               cAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                     >  1:445081/1‑150 (MQ=255)
                                                                               cAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                     >  1:39890/1‑150 (MQ=255)
                                                                                aaGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                                    >  2:831900/1‑150 (MQ=255)
                                                                                aaGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                                    >  2:384999/1‑150 (MQ=255)
                                                                                 aatttaaaaTCCGGCTTCCGTTCGCACTCCTGTTTTGGCTCCGCCATTAAAAAGCCAGTTTTTATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCa                                                                      <  1:450326/142‑1 (MQ=255)
                                                                                 aatttaaaaTCCGCCATCCTTTCGAAATCCTGTTTTGGTTTCCCTATTAATAGGCCAGTTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCa                                                                      <  2:688442/142‑1 (MQ=255)
                                                                                   ttCAAATACGGCATCGTTTCGAAATCCGGTTTGGGCTTCCCCATTAATAGGCCAGTTTTTACCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGa                                                                 <  2:61861/150‑1 (MQ=255)
                                                                                       aaaTCCGCCTTCCTTTCCCAACCCTGTTTTGGCTTCCCCTTTAATAACCCAGTTTTTATCCCCTTTTTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                                <  2:149117/148‑1 (MQ=255)
                                                                                       aaaTACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                             >  1:1763/1‑151 (MQ=255)
                                                                                         aTACGCCATCCTTCCGCAACCCTGTTTTGGTTCCGCTTTTAAAAACCAAGTTTTTACCACTTTTATTCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaaca                                                            <  2:144180/150‑1 (MQ=255)
                                                                                              gcctccttttCCCAATCCTGGTTGGGCTCCCCCTTTAAAAAACCCGGTTTAATCCCTTTTATCCCAAGCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGa                                                          <  1:297411/140‑1 (MQ=255)
                                                                                              gccATCCTTTCCCATTCCTGTTTTGGCTTCCCAATAAATAGGCCGGTTTTTATCCCTTTTATCAAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaaa                                                        <  1:142687/148‑1 (MQ=255)
                                                                                                 acccGTTCCCATTCCTGTTTTGCCTTCGCATTAAATAAGCCGGGTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACt                                                      <  2:118953/147‑1 (MQ=255)
                                                                                                 aTCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                     >  2:811126/1‑150 (MQ=255)
                                                                                                   cctttCGCAATCCTGGTTTGGCTTCCCCATAAAAAAGCCGGTTTTTACCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                     <  2:672981/145‑1 (MQ=255)
                                                                                                   ccGTCCGCAATCCTGGTTTGGTTCCCCCATTAAAAACCCGGGTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTa                                                    <  2:366740/149‑1 (MQ=255)
                                                                                                      ttCGCAATCCTGGTTGGGCTTCCCCATTAAAAAGCCAGTTTTTACCCCTTTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                                  <  2:20761/148‑1 (MQ=255)
                                                                                                      ttCGCAACCTGGTTTTGGTTTCGCCATAAATAAGCCAGGTTTTATCCCCTTTATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                                <  2:982355/150‑1 (MQ=255)
                                                                                                      tcccaaaTCCTGTTTTGGCTCCCCCATAAATAAGCCAGGTTTTATCACTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAg                                                           <  1:328433/134‑1 (MQ=255)
                                                                                                       tCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                                >  1:287453/1‑149 (MQ=255)
                                                                                                          cAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaagaa                                            >  2:710806/1‑150 (MQ=255)
                                                                                                                tGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCAATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtct                                     >  1:853119/1‑151 (MQ=255)
                                                                                                                 ggTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtctt                                    >  1:3136/1‑151 (MQ=255)
                                                                                                                    tttGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                                  >  2:712490/1‑151 (MQ=255)
                                                                                                                     ttGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                                 >  2:66753/1‑151 (MQ=255)
                                                                                                                          cTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttt                             >  2:1056409/1‑151 (MQ=255)
                                                                                                                           ttCGCCATAAAAAACCCGGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                                <  2:1024751/147‑1 (MQ=255)
                                                                                                                             cGCCATTAATAAGCCGGGTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttttt                           <  2:1012575/150‑1 (MQ=255)
                                                                                                                                  ttAAAAACCCAGTTTTTATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                       <  2:880906/149‑1 (MQ=255)
                                                                                                                                    aaaaacccggtttttACCCCTTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt                        <  1:897822/135‑1 (MQ=255)
                                                                                                                                     ataaGCCGGTTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt                  <  2:754722/151‑1 (MQ=255)
                                                                                                                                          ccAGTTTTTACCCCCTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt               <  1:1055588/149‑1 (MQ=255)
                                                                                                                                            aGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc           >  1:801703/1‑151 (MQ=255)
                                                                                                                                             gttttaATCCCTTTAATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCa          <  1:537917/149‑1 (MQ=255)
                                                                                                                                               tttaATCCCCTTTTTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGa             <  1:1052729/146‑1 (MQ=255)
                                                                                                                                                     cACTTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAttt   <  2:578134/150‑1 (MQ=255)
                                                                                                                                                     cACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtttt  >  2:909162/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
CCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537964‑4538258

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: