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→GGGnanMN‑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

ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537960‑4538258
                                                                                                                                                          |                                                                                                                                                     
aTCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCAc                                                                                                                                                      <  2:75452/151‑1 (MQ=255)
  ccccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCt                                                                                                                                                    <  2:559484/151‑1 (MQ=255)
    ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCttt                                                                                                                                                  <  1:852557/151‑1 (MQ=255)
     cGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCttt                                                                                                                                                  >  1:901854/1‑150 (MQ=255)
     cGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCttt                                                                                                                                                  >  2:104925/1‑150 (MQ=255)
     cGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTa                                                                                                                                                 <  2:365504/151‑1 (MQ=255)
     cGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTa                                                                                                                                                 <  2:639560/151‑1 (MQ=255)
         aaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTca                                                                                                                                             <  1:1574/151‑1 (MQ=255)
            ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACaa                                                                                                                                          <  2:870082/151‑1 (MQ=255)
              gCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAAcc                                                                                                                                        <  1:743425/151‑1 (MQ=255)
                cAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAAccgc                                                                                                                                      <  1:410121/151‑1 (MQ=255)
                      cGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCAcc                                                                                                                                >  2:600979/1‑151 (MQ=255)
                       gCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGCTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGGTTTATCACCTTTATTCACAACCGCCGCAc                                                                                                                                 >  2:192764/1‑149 (MQ=255)
                                 gatgaGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGt                                                                                                                     >  2:731685/1‑151 (MQ=255)
                                   tgaGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTa                                                                                                                    >  2:655418/1‑150 (MQ=255)
                                     agaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTAcc                                                                                                                  >  2:200019/1‑150 (MQ=255)
                                         gggAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGc                                                                                                             >  2:261582/1‑151 (MQ=255)
                                         gggAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGc                                                                                                             >  2:857281/1‑151 (MQ=255)
                                          ggAGCAAGTTTTTTCCATTTAAATTTTTTCCCGGGAAAATAAAGTTAAAAACCGCCACCCTTCCGCACCCCGGTTTGGCCTCCGCCTTAAAAAACCCGGTTTTTCCCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                            <  2:395856/151‑1 (MQ=255)
                                            aGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGAt                                                                                                          >  1:24734/1‑151 (MQ=255)
                                             gCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGAtt                                                                                                         >  1:894608/1‑151 (MQ=255)
                                              cAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                                        >  1:10422/1‑151 (MQ=255)
                                               aaGTTTATTCCATTTTAAATTATTACCGGTGAAATAAATTTAAAATCCGGCATCCGTTCGCAACCCTGGTTGGGCTTCGCCATAAATAAGCAAGGTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg                                                                                                       <  1:478547/151‑1 (MQ=255)
                                                   tttttCCTTTTTAAATTTTGCCGGGTGAAACAATGTCAAAATCCGCCCTCCTTCCGCACCCCTGTTTTGGCTCCCCAATAAAAAGCCAAGGTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                                        <  2:149664/143‑1 (MQ=255)
                                                   tttttCCATTTAAAATTATTCCGGGTAAAATAATTTTAAAATCCGCCTCCCTTCCACAATCCTGGTTTGGTTTCCCCATAAAAACGCCGGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                     <  2:197969/146‑1 (MQ=255)
                                                     aTTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTCGGCTTCGCCATTAATAAGCCAGGCTTTATCACCTTTATTCACAAACGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAg                                                                                                   >  1:810948/1‑149 (MQ=255)
                                                       tCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt                                                                                                >  1:710644/1‑150 (MQ=255)
                                                                aaTTATTCCGGGTGAAATCAAGTACAAATCCGGCACCGTTTCCCACTCCTGTTTGGGTTCCCCCTTAAATAACCCAGGTTTTATCCCTTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTg                                                                                      <  2:885327/151‑1 (MQ=255)
                                                                   tattaCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg                                                                                   >  1:947531/1‑151 (MQ=255)
                                                                          ggTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                              >  1:91353/1‑149 (MQ=255)
                                                                           gaaaaacaaatttCAATTCCGGCCTCCTTTCCCAACCCTGGTTGGGTTCCCCCTTTAAAAAGCATGTTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAgg                                                                           <  2:371234/140‑1 (MQ=255)
                                                                                  cAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                     >  1:551258/1‑150 (MQ=255)
                                                                                   aaGTTAAAATACGGCATCCGTTCGCAATCCTGTTTTGGCTTCGCCATTAATAAGCCAGGTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                     <  1:77184/149‑1 (MQ=255)
                                                                                    aGTTAAAATACGGCATCCGTCCGCAATCCTGGTTTGGCTCCGCCATTAATAAGCCAGGTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAg                                                                  <  1:716994/151‑1 (MQ=255)
                                                                                        aaaaTGCGGCACCCTTTCCCACCCCTGTTTGGCCTTCCCAATAAATACGCCGGGTTTTATCCCCTTTTTCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGaca                                                               <  2:642220/149‑1 (MQ=255)
                                                                                         aaaTAGGGAGTCCGTTCAAAATCCAGTTTTGCCTACGCAAGTAAGAAGCCTGTTTTTATCCCTTTTTTACAAAGCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                             <  1:358709/151‑1 (MQ=255)
                                                                                          aaTACGGCACCGGTTCGCATCTCTGTTTTGGCTCCGCATTTAAAACGCCAGTTTTATCCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaaca                                                            <  2:923251/151‑1 (MQ=255)
                                                                                           atccGGCATCCGTTCGCATTCCTGTTTTGGCTCCGCAATTAAAAAGCCGGTTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAg                                                           <  2:490258/148‑1 (MQ=255)
                                                                                           aTACGGCATCCGTTCGCAATCCTGTTTTGGCTTCGCCATTAATAAGCCAGTTTTTACCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAg                                                           <  1:828728/151‑1 (MQ=255)
                                                                                               gccacccGTTCCAATTCCTGTTTGGGCTTCGCCTTAAATAAGCCAGTTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAg                                                           <  2:323562/142‑1 (MQ=255)
                                                                                                  aTCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                     >  1:327161/1‑150 (MQ=255)
                                                                                                    ccGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                                  >  1:694230/1‑151 (MQ=255)
                                                                                                    ccGTCCGCAATCCTGGTTTGGCTCCGCCATTAATAAGCCAGGTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAt                                                   <  1:649247/150‑1 (MQ=255)
                                                                                                    ccGTCCGAAACCCTGTTTTGGCTCCGCCATAAATAACCCAGTTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTa                                                    <  2:171972/149‑1 (MQ=255)
                                                                                                     cGTTCGCAATCCGGGTTGGGCTTCGCCATTAATAAGCCAGTTTTTATCCCTTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                                 <  1:931723/151‑1 (MQ=255)
                                                                                                          gCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                             >  1:745688/1‑150 (MQ=255)
                                                                                                               ccTGTTTGGGTTCCGCCATTAAAAACCAAGTTTTTATCCCCTTTATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGaaa                                           <  2:686501/147‑1 (MQ=255)
                                                                                                                  ggTTTGGCTCCGCCATTAATAAGCCAGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaata                                         <  1:960824/146‑1 (MQ=255)
                                                                                                                     ttGGCTTTCGCCATAAATAGGCCAGTTTTTATCCCTTTTTTTCACAGCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtc                                      <  1:397041/147‑1 (MQ=255)
                                                                                                                     ttGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                                 >  2:23909/1‑151 (MQ=255)
                                                                                                                        gCTCCGCTATTAAAAAGCAAGGTTTTATCTCTTTTCTTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                                <  1:828214/149‑1 (MQ=255)
                                                                                                                            ttCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt                            >  1:505325/1‑151 (MQ=255)
                                                                                                                                 cATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                       >  1:741202/1‑151 (MQ=255)
                                                                                                                                   ttAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                      >  1:878806/1‑150 (MQ=255)
                                                                                                                                   ttAATAACCCAGTTTTTCCCCCCTTTATACACACCCGCCGCCCCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                       <  1:561352/149‑1 (MQ=255)
                                                                                                                                    ttaaaagcccAGTTTTTATCCCTTTTATTCCAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                       <  1:857294/141‑1 (MQ=255)
                                                                                                                                      aataaGCCAGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa                    <  1:446236/150‑1 (MQ=255)
                                                                                                                                      aataaCCCAGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa                    <  1:201567/150‑1 (MQ=255)
                                                                                                                                       aaaatccAGTTTATCTCCCCTTTATACCAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt                  <  2:871655/146‑1 (MQ=255)
                                                                                                                                       aaaaggcgatttttaATCACTTTTATTCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa                    <  2:227812/139‑1 (MQ=255)
                                                                                                                                       aaaaggccggtttttATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaaa                   <  1:602751/140‑1 (MQ=255)
                                                                                                                                            ccAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg              >  2:226043/1‑150 (MQ=255)
                                                                                                                                             ccgggTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg              <  2:396265/147‑1 (MQ=255)
                                                                                                                                             cAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg              >  1:228154/1‑149 (MQ=255)
                                                                                                                                               gggTTTTATCACCCTTATTCACAACCGCCGCCCCAGCCGTTCTGCACGAGGGCGATTCTCCAGCGTAACTCCATTGCTGTGCTGAGGGAACAAAAGACATCCGAGACTTATTGAAGACTTATTCTTCAGCTTTTTTGACTAACTAGTGAGc           >  2:122060/2‑151 (MQ=255)
                                                                                                                                                 tttttACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt               <  1:2873/144‑1 (MQ=255)
                                                                                                                                                  ttttATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAt         <  1:893109/150‑1 (MQ=255)
                                                                                                                                                    ttATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg       >  1:501665/1‑150 (MQ=255)
                                                                                                                                                    ttATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg       >  1:923817/1‑150 (MQ=255)
                                                                                                                                                          aCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtttt  >  1:438775/1‑150 (MQ=255)
                                                                                                                                                          |                                                                                                                                                     
ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  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: