Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,538,1090AC42.6% 104.8 / 28.5 69G219G (GGT→GGG) nanMN‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (29/10);  new base C (0/29);  total (29/40)
Fisher's exact test for biased strand distribution p-value = 4.87e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537960‑4538258
                                                                                                                                                          |                                                                                                                                                     
aTCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CAc                                                                                                                                                      <  2:75452/151‑1 (MQ=255)
  ccccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCt                                                                                                                                                    <  2:559484/151‑1 (MQ=255)
    ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCttt                                                                                                                                                  <  1:852557/151‑1 (MQ=255)
     cGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCttt                                                                                                                                                  >  1:901854/1‑150 (MQ=255)
     cGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCttt                                                                                                                                                  >  2:104925/1‑150 (MQ=255)
     cGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTa                                                                                                                                                 <  2:365504/151‑1 (MQ=255)
     cGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTa                                                                                                                                                 <  2:639560/151‑1 (MQ=255)
         aaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTca                                                                                                                                             <  1:1574/151‑1 (MQ=255)
            ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACaa                                                                                                                                          <  2:870082/151‑1 (MQ=255)
              gCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACAAcc                                                                                                                                        <  1:743425/151‑1 (MQ=255)
                cAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACAAccgc                                                                                                                                      <  1:410121/151‑1 (MQ=255)
                      cGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACAACCGCCGCAcc                                                                                                                                >  2:600979/1‑151 (MQ=255)
                       gCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGCTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGGTTTAT‑CACCTTTATTCACAACCGCCGCAc                                                                                                                                 >  2:192764/1‑149 (MQ=255)
                                 gatgaGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGt                                                                                                                     >  2:731685/1‑151 (MQ=255)
                                   tgaGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTa                                                                                                                    >  2:655418/1‑150 (MQ=255)
                                     agaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTAcc                                                                                                                  >  2:200019/1‑150 (MQ=255)
                                         gggAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGc                                                                                                             >  2:261582/1‑151 (MQ=255)
                                         gggAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGc                                                                                                             >  2:857281/1‑151 (MQ=255)
                                          ggAGCAAGTTTTTTCCATTTAAATTTTTTCCCGGGAAAATAAAGTTAAAAACCGCCACCCTTCCGC‑ACCCCGGTTTGGCC‑TCCGCCTTA‑AAAAACCC‑GGTTTTTCC‑CCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                            <  2:395856/151‑1 (MQ=255)
                                            aGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGAt                                                                                                          >  1:24734/1‑151 (MQ=255)
                                             gCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGAtt                                                                                                         >  1:894608/1‑151 (MQ=255)
                                              cAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                                        >  1:10422/1‑151 (MQ=255)
                                               aaGTTTATTCCATTTTAAATTATTACCGGTGAAATAAATTTAAAATCCGGCATCCGTTCGCAACCCTGGTTGGG‑C‑TTCGCCATA‑AATAAGCA‑AGGTTTTAT‑CCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg                                                                                                       <  1:478547/151‑1 (MQ=255)
                                                   tttttCCTTTTTAAATTTTGCCGGGTGAAACAATGTCAAAATCCGCCCTCCTTCCGCACCCCTGTTTTGG‑C‑‑TCCCCAAT‑AAAAAGCCAAGGTTTTAC‑CCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                                        <  2:149664/143‑1 (MQ=255)
                                                   tttttCCATTTAAAATTATTCCGGGTAAAATAATTTTAAAATCCGCCTCCCTTCCACAATCCTGGTTTGG‑T‑TTCCCCATA‑AAAACGCC‑GGTTTTTAT‑CCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                     <  2:197969/146‑1 (MQ=255)
                                                     aTTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTCGG‑C‑TTCGCCATT‑AATAAGCC‑AGGCTTTAT‑CACCTTTATTCACAAACGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAg                                                                                                   >  1:810948/1‑149 (MQ=255)
                                                       tCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt                                                                                                >  1:710644/1‑150 (MQ=255)
                                                                aaTTATTCCGGGTGAAATCAAGTACAAATCCGGCACCGTTTCCCACTCCT‑GTTTGG‑G‑TTCCCCCTTAAATAACCC‑AGGTTTTAT‑CCCTTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTg                                                                                      <  2:885327/151‑1 (MQ=255)
                                                                   tattaCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg                                                                                   >  1:947531/1‑151 (MQ=255)
                                                                          ggTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                              >  1:91353/1‑149 (MQ=255)
                                                                           gaaaaacaaatttCAATTCCGGCCTCCTTTCCCAACCCTGGTTGGG‑T‑TCCCCCTTT‑AAAAAGCA‑TGTTTTTAC‑CCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAgg                                                                           <  2:371234/140‑1 (MQ=255)
                                                                                  cAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                     >  1:551258/1‑150 (MQ=255)
                                                                                   aaGTTAAAATACGGCATCCGTTCGCAATCCTGTTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAC‑CCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                     <  1:77184/149‑1 (MQ=255)
                                                                                    aGTTAAAATACGGCATCCGTCCGCAATCCTGGTTTGG‑C‑TCCGCCATT‑AATAAGCC‑AGGTTTTAT‑CCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAg                                                                  <  1:716994/151‑1 (MQ=255)
                                                                                        aaaaTGCGGCACCCTTTCCCACCCCTGTTTGGC‑C‑TTCCCAATA‑AATACGCC‑GGGTTTTAT‑CCCCTTTTTCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGaca                                                               <  2:642220/149‑1 (MQ=255)
                                                                                         aaaTAGGGAGTCCGTTCAAAATCCAGTTTTGC‑C‑TACGCAAGT‑AAGAAGCC‑TGTTTTTAT‑CCCTTTTTTACAAAGCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                             <  1:358709/151‑1 (MQ=255)
                                                                                          aaTACGGCACCGGTTCGCATCTCTGTTTTGG‑C‑TCCGCATTT‑AAAACGCC‑A‑GTTTTATCCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaaca                                                            <  2:923251/151‑1 (MQ=255)
                                                                                           atccGGCATCCGTTCGCATTCCTGTTTTGG‑C‑TCCGCAATT‑AAAAAGCC‑GGTTTTTAC‑CCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAg                                                           <  2:490258/148‑1 (MQ=255)
                                                                                           aTACGGCATCCGTTCGCAATCCTGTTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGTTTTTAC‑CCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAg                                                           <  1:828728/151‑1 (MQ=255)
                                                                                               gccacccGTTCCAATTCCTGTTTGGG‑C‑TTCGCCTTA‑AATAAGCC‑AGTTTTTAT‑CCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAg                                                           <  2:323562/142‑1 (MQ=255)
                                                                                                  aTCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                     >  1:327161/1‑150 (MQ=255)
                                                                                                    ccGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                                  >  1:694230/1‑151 (MQ=255)
                                                                                                    ccGTCCGCAATCCTGGTTTGG‑C‑TCCGCCATT‑AATAAGCC‑AGGTTTTAT‑CCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAt                                                   <  1:649247/150‑1 (MQ=255)
                                                                                                    ccGTCCGAAACCCTGTTTTGG‑C‑TCCGCCATA‑AATAACCC‑AGTTTTTAC‑CCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTa                                                    <  2:171972/149‑1 (MQ=255)
                                                                                                     cGTTCGCAATCCGGGTTGGG‑C‑TTCGCCATT‑AATAAGCC‑AGTTTTTAT‑CCCTTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                                 <  1:931723/151‑1 (MQ=255)
                                                                                                          gCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                             >  1:745688/1‑150 (MQ=255)
                                                                                                               ccTGTTTGGG‑T‑TCCGCCATT‑AAAAACCA‑AGTTTTTAT‑CCCCTTTATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGaaa                                           <  2:686501/147‑1 (MQ=255)
                                                                                                                  ggTTTGG‑C‑TCCGCCATT‑AATAAGCC‑AGTTTTTAT‑CCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaata                                         <  1:960824/146‑1 (MQ=255)
                                                                                                                     ttGG‑CTTTCGCCATA‑AATAGGCC‑AGTTTTTAT‑CCCTTTTTTTCACAGCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtc                                      <  1:397041/147‑1 (MQ=255)
                                                                                                                     ttGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                                 >  2:23909/1‑151 (MQ=255)
                                                                                                                        g‑C‑TCCGCTATT‑AAAAAGCA‑AGGTTTTAT‑CTCTTTTCTTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                                <  1:828214/149‑1 (MQ=255)
                                                                                                                            ttCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt                            >  1:505325/1‑151 (MQ=255)
                                                                                                                                 cATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                       >  1:741202/1‑151 (MQ=255)
                                                                                                                                   tt‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                      >  1:878806/1‑150 (MQ=255)
                                                                                                                                   tt‑AATAACCC‑AGTTTTTCC‑CCCCTTTATACACACCCGCCGCCCCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                       <  1:561352/149‑1 (MQ=255)
                                                                                                                                    ttaaaagccc‑AGTTTTTAT‑CCCTTTTATTCCAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                       <  1:857294/141‑1 (MQ=255)
                                                                                                                                      aataaGCC‑AGTTTTTAT‑CCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa                    <  1:446236/150‑1 (MQ=255)
                                                                                                                                      aataaCCC‑AGTTTTTAT‑CCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa                    <  1:201567/150‑1 (MQ=255)
                                                                                                                                       aaaatcc‑AGTTTATCT‑CCCCTTTATACCAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt                  <  2:871655/146‑1 (MQ=255)
                                                                                                                                       aaaaggcgatttttaAT‑CACTTTTATTCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa                    <  2:227812/139‑1 (MQ=255)
                                                                                                                                       aaaaggccggtttttAT‑CCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaaa                   <  1:602751/140‑1 (MQ=255)
                                                                                                                                            cc‑AGGTTTTAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg              >  2:226043/1‑150 (MQ=255)
                                                                                                                                             ccgggTTTTAC‑CCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg              <  2:396265/147‑1 (MQ=255)
                                                                                                                                             c‑AGGTTTTAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg              >  1:228154/1‑149 (MQ=255)
                                                                                                                                               gggTTTTAT‑CACCCTTATTCACAACCGCCGCCCCAGCCGTTCTGCACGAGGGCGATTCTCCAGCGTAACTCCATTGCTGTGCTGAGGGAACAAAAGACATCCGAGACTTATTGAAGACTTATTCTTCAGCTTTTTTGACTAACTAGTGAGc           >  2:122060/2‑151 (MQ=255)
                                                                                                                                                 tttttAC‑CCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt               <  1:2873/144‑1 (MQ=255)
                                                                                                                                                  ttttAT‑CCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAt         <  1:893109/150‑1 (MQ=255)
                                                                                                                                                    ttAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg       >  1:501665/1‑150 (MQ=255)
                                                                                                                                                    ttAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg       >  1:923817/1‑150 (MQ=255)
                                                                                                                                                          aCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtttt  >  1:438775/1‑150 (MQ=255)
                                                                                                                                                          |                                                                                                                                                     
ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATT‑AATAAGCC‑AGGTTTTAT‑CACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537960‑4538258

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