Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,538,1090AC45.2% 136.1 / 29.9 84G219G (GGT→GGG) nanMN‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (36/10);  new base C (0/38);  total (36/48)
Fisher's exact test for biased strand distribution p-value = 5.79e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATT‑ACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537962‑4538258
                                                                                                                                                        |                                                                                                                                                       
ccccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATT‑ACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCAcc                                                                                                                                                       <  2:292716/150‑1 (MQ=255)
  ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATT‑ACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑ttn                                                                                                                                                   <  2:1311437/151‑2 (MQ=255)
    gCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATT‑ACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑At                                                                                                                                                <  2:1259557/151‑1 (MQ=255)
     cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATT‑ACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑Att                                                                                                                                               <  1:110191/151‑1 (MQ=255)
      aaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATT‑ACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTc                                                                                                                                              <  2:1059761/151‑1 (MQ=255)
       aaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATT‑ACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTca                                                                                                                                             <  1:308718/151‑1 (MQ=255)
        aaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATT‑ACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTcac                                                                                                                                            <  2:1279758/151‑1 (MQ=255)
          ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATT‑ACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACaa                                                                                                                                          <  2:845815/151‑1 (MQ=255)
           cGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATT‑ACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAAc                                                                                                                                         <  2:1225237/151‑1 (MQ=255)
                     gCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCTTTTTAAA‑TTATT‑ACCGGTGAAATCAATTTCAAATACGGCATCCGTTCGCAAT‑CCTGTTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCCCC‑TTT‑ATTCACAACCGCCGCACCa                                                                                                                               <  2:4760/151‑1 (MQ=255)
                      ccATCTGGTGATGAGACGGGACCAAGTTTATTCCATTTTAAA‑TTATT‑ACCGGTGAAATCAAGTTAAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGTTTTT‑ATCCCT‑TTT‑ATTCACAACCGCCGCACCAg                                                                                                                              <  2:1468652/151‑1 (MQ=255)
                                  gagcCGGGAGCAAGCTTTTTCCATTTAAAA‑TTATT‑CCGGGTAAAATAAAGTTCAAATCCGGCTTCCGTTCGCAAT‑CCTGTTTTGCCTTCGCC‑ATTAATAAGCCAGTTTTT‑ATCCCC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTAc                                                                                                                   <  2:852631/150‑1 (MQ=255)
                                  gagaCGGGAGCAAGCTTATTCCATTTTAAA‑TTATT‑ACCGGTGAAATCAAGTCCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTAcc                                                                                                                  >  2:1290698/1‑151 (MQ=255)
                                   agaCGGGAGCAAGCTTATTCCATTTTAAA‑TTATT‑ACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTAcc                                                                                                                  >  2:355050/1‑150 (MQ=255)
                                         gAGCAAGCTTATTCCATTTTAAA‑TTATT‑ACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                            >  2:731734/1‑150 (MQ=255)
                                          aGCAAGCTTATTCCATTTTAAA‑TTATT‑ACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGAt                                                                                                          >  2:911462/1‑151 (MQ=255)
                                          aGCAAGCTTATTCCATTTTAAA‑TTATT‑ACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGAt                                                                                                          >  2:842710/1‑151 (MQ=255)
                                           gCAAGCTTTTTCCATTTTAAA‑TTATT‑CCCGGGAAAATCAAGTTAAAATACGGCATCCGTTCGCAAT‑CCTGTTTTGGCTTCGCC‑ATTAAAAAGCCAGGTTTT‑ATCCCC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGAtt                                                                                                         <  1:1063045/151‑1 (MQ=255)
                                           acaacttttttCCTTTTTAAA‑TTTTT‑CCCGGTAAAATCAATTTAAAATCCGCCCCCCGTTCCAAAC‑CCGGGTTTGGCTTCCCC‑ATAAAAAAGCCAGTTTTT‑ATCCCT‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGAtt                                                                                                         <  2:698728/142‑1 (MQ=255)
                                            cAAGCTTATTCCATTTTAAA‑TTATT‑ACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                                        >  1:1030066/1‑151 (MQ=255)
                                             aaccTTATTTCCTTTTAAA‑ATATTCTCCGGTCAATTAAATTTCAAATACGTGCCCCATCCGCACT‑CCTGTTTTGCCTTCGCC‑TAAAATACGCCGGTTTTT‑CTACCT‑TTT‑CTTCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                                        <  2:536175/148‑1 (MQ=255)
                                               gttttttCCTTTTTAAT‑TTATT‑ACCGGGAAAACAAATTTAAAATCCGGCTCCCGTCCGCATC‑CGTGGTTTGCCTCCCCT‑ATAAAAACGCCGGTTTTT‑ACCCCC‑TTT‑TCCCCCACCCCCCGCCCCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                     <  2:1432378/146‑1 (MQ=255)
                                               gCTTATTCCATTTAAAA‑TTATT‑CCCGGTAAAATCAAGTTCAAATCCGGCATCCTTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGTTTTT‑ATCCCC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                     <  1:1322353/151‑1 (MQ=255)
                                                ttttttCCTTTTTAAT‑TTATT‑CCCGGTAAAATCAATTTAAAATACGCCCTCCTTTCGCAAT‑CTTGGTTTGCCTTCCCC‑ATAAATAACCCGGTTTTT‑ATCCCT‑TTT‑ATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                     <  2:262376/146‑1 (MQ=255)
                                                 ttATTCCATTTTAAA‑TTATT‑ACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                     >  1:1155142/1‑149 (MQ=255)
                                                    tttcATTTTAAATTTATT‑CCCGGTAAAATAAATTTAAAATCCGGC‑TTCGTTCGCCATCCCTGTTTTGTTTCCGCC‑TTAAATAAGCC‑TGTTTT‑TTCCCC‑TTTAATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGc                                                                                                  <  2:242535/147‑1 (MQ=255)
                                                     tCCTTTTAAAT‑TTATT‑CCGGGTAAAATAAATTTAAAATCGGGCATCCTTCCGCATT‑CTTGGTTTGTTTTCCCA‑ATTAAAAAGCC‑GGTTTT‑TTCCCCTTTT‑ATTCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt                                                                                                <  2:1201607/150‑1 (MQ=255)
                                                     tCCATTTTAAA‑TTATT‑ACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt                                                                                                >  2:1400069/1‑150 (MQ=255)
                                                     cccATTTTAAA‑TTATG‑CCCGGTAAAACCAATTTAAAAAACGGCATCCTTTCGCAAT‑CCTGTTTGGGCTCCGCC‑A‑TAAAGAGCCCGGTTTTAACCACC‑TTT‑ATTCCAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa                                                                                               <  2:41206/150‑1 (MQ=255)
                                                         tttaAAA‑TTATT‑CCGGTTGAAATCAATTTAAAATCCGCCATCCTTTCGCATC‑CTTGTTTTGACTTCGCC‑ATAAATAACCCAGGTTTT‑ACCCCT‑TTT‑TTCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAAc                                                                                             <  1:83313/149‑1 (MQ=255)
                                                           ttAAA‑TTATT‑CCCGGTAAAATCAATTTAAAATCCGCCATCCTTTCGAAAT‑CCTGGTTTGGCTTCCCC‑ATTAATAAGCCAGGTTTT‑ACCCCC‑TTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                          <  2:1417739/150‑1 (MQ=255)
                                                             aaa‑TTATT‑ACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                       >  1:764039/1‑151 (MQ=255)
                                                                  tttttcccGGTAAAATAAAGTTAAAATACGCCTTCCGTCCCCAAT‑CCGGTTTGGGTTTCCCC‑TTTAATAGGCCAGTTTTT‑ATCCCC‑TTT‑ATCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGc                                                                                     <  2:101830/144‑1 (MQ=255)
                                                                   att‑cCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGTTTTGGCTTCCCC‑ATTAATAAGCCAGGTTTT‑ATCCCC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgt                                                                                  <  1:601004/151‑1 (MQ=255)
                                                                    tt‑aCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgtg                                                                                 >  2:183225/1‑151 (MQ=255)
                                                                       aCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt                                                                               >  1:738424/1‑151 (MQ=255)
                                                                        ccGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                              >  1:469163/1‑151 (MQ=255)
                                                                             gAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                          >  2:235763/1‑150 (MQ=255)
                                                                             aaaaaCAATTTAAAATACGCCACCCTTTCCCACC‑CCTGGTTGGGCTTCGCC‑ATAAATAACCCAGGTTTT‑ATCCCC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                          <  2:45318/149‑1 (MQ=255)
                                                                                 aaaaGTTCAAATCCGGCTCCCTTTCGCAAC‑CCTGTTTTGGCTGCGCT‑ATTAAAAACCCAGTTTTT‑TCCCCC‑TTT‑ATACGCAACCGCCGCTCCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGTATCaa                                                                     <  2:1099623/149‑1 (MQ=255)
                                                                                  cAATTTAAAATCCGGCACCCGTTCGCAAT‑CCTGTTTTGGCTTCGCC‑ATTAATAAGCCAGTTTTT‑ACCCCC‑TTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                     <  2:46028/150‑1 (MQ=255)
                                                                                  cAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                     >  2:219701/1‑150 (MQ=255)
                                                                                  aaaGTTAAAATCCGGCATCCGTTCGCAAT‑CCTGTTTTGGCTTCGCC‑ATTAATAACCCAGTTTTT‑ATCCCC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                                    <  2:187670/150‑1 (MQ=255)
                                                                                      ttaaaagccGGCACCCGTTCGCATT‑CCTGTTTTGGCTTCGCC‑ATTAATAAGCCAGTTTTT‑ATCCCT‑TTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGa                                                                 <  2:1391390/142‑1 (MQ=255)
                                                                                       taaaaTCCGGCATCCGTCCCCAAT‑CCTGGTTGGGCTTCGCC‑ATAAATAAGCCAGTTTTT‑ATCCCT‑TTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGan                                                                <  2:1/148‑2 (MQ=255)
                                                                                         aaaTACGGCATCCGTTCTCAAT‑CCTGGTTTTGCTTCGCC‑ATTAATAAACAAGGGTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCATCCGTTCCGGACCAGGGCGTTTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                             >  2:673960/1‑151 (MQ=255)
                                                                                         aaaTACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                             >  1:1194728/1‑151 (MQ=255)
                                                                                            agcGGCATCCGTTCACAA‑‑CCAGGGTTGGTTTGGCCAATTAATACGCCAGTTTTC‑AACCCC‑TTT‑TTCCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGa                                                          <  2:716607/149‑1 (MQ=255)
                                                                                             ccGGCTTCCTTTCGCAAT‑CCTGTTTTGGCTTCCCC‑ATTAATAACCCAGTTTTT‑ATCCCC‑TTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa                                                         <  2:912931/150‑1 (MQ=255)
                                                                                              cGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa                                                         >  1:687693/1‑150 (MQ=255)
                                                                                                 cATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGTTTTT‑ATCCCC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                     <  1:184364/151‑1 (MQ=255)
                                                                                                  accggTTCGCAAT‑CCTGGTTTGGCTTCCCC‑ATTAATAAGCCGGGTTTA‑ATCCCC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTa                                                    <  1:108465/147‑1 (MQ=255)
                                                                                                  aTCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                     >  1:467838/1‑150 (MQ=255)
                                                                                                      gTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                                 >  1:358968/1‑150 (MQ=255)
                                                                                                        ccGCAAC‑CCTGTTTTGGTTCCGCC‑TTTAATAAGCCAGTTTTT‑ACCCCT‑TTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaag                                              <  1:517023/150‑1 (MQ=255)
                                                                                                         cccaaacccTGTTTGGGCTTCGCA‑TAAAAAAAGCCAGTTTTT‑ATCCCC‑TTA‑ATCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaag                                              <  1:189337/144‑1 (MQ=255)
                                                                                                           gcatcccggttttGGCTGCGCC‑ATTAATAAGCCAGTTTTT‑ATCCCC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaagaa                                            <  1:132267/141‑1 (MQ=255)
                                                                                                            aaT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCCCC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                          <  1:56049/151‑1 (MQ=255)
                                                                                                            aaT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                          >  2:135138/1‑151 (MQ=255)
                                                                                                              acccTGTTTTGGCTTCCCC‑ATAAATAAGCCAGTTTTT‑ATCCCC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaata                                         <  1:919392/149‑1 (MQ=255)
                                                                                                                  tGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtct                                     >  1:1303128/1‑151 (MQ=255)
                                                                                                                   ggTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtctt                                    >  1:1330494/1‑151 (MQ=255)
                                                                                                                      ttGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                                 >  2:603261/1‑151 (MQ=255)
                                                                                                                       gggCTTCCCC‑ATTAATAAGCCGGCTTTT‑ATCCCT‑TTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                                <  1:623629/150‑1 (MQ=255)
                                                                                                                        ggCTTCGCC‑ATTAAAAAGCCAGGTTTT‑ATCCCT‑TTT‑ATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                                <  1:747694/150‑1 (MQ=255)
                                                                                                                          cTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttt                             >  1:1297536/1‑151 (MQ=255)
                                                                                                                            tcccc‑ATTAATAAGCCAGTTTTT‑ACCCCC‑TTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt                            <  2:993935/147‑1 (MQ=255)
                                                                                                                            tCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttttt                           >  2:1197805/1‑151 (MQ=255)
                                                                                                                             cGCC‑ATAAATAAGCCAGGTTTT‑ATCCCC‑TTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttttt                           <  1:1468515/150‑1 (MQ=255)
                                                                                                                              gCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg                         >  2:250548/1‑151 (MQ=255)
                                                                                                                               cc‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt                        >  1:1228685/1‑151 (MQ=255)
                                                                                                                                c‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                       >  1:1347296/1‑151 (MQ=255)
                                                                                                                                c‑ATTAAAAAGCGAGTTTTT‑ATCCCC‑TTT‑ATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                       <  1:338457/151‑1 (MQ=255)
                                                                                                                                  aTTAACAAGCAAGTTTTT‑GTCCCT‑TTT‑TTTCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                       <  2:351705/150‑1 (MQ=255)
                                                                                                                                   ttAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                      >  2:532919/1‑150 (MQ=255)
                                                                                                                                   ttAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                      >  1:1286549/1‑150 (MQ=255)
                                                                                                                                    taataaGCCGGTTTTT‑ATCCCC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa                     <  1:705472/150‑1 (MQ=255)
                                                                                                                                     aataaGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa                    >  1:440890/1‑150 (MQ=255)
                                                                                                                                     aataaGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa                    >  1:1197016/1‑150 (MQ=255)
                                                                                                                                     aataaGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa                    >  2:1016938/1‑150 (MQ=255)
                                                                                                                                          gCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg              >  2:758182/1‑151 (MQ=255)
                                                                                                                                           ccAGTTTTT‑ATCCCC‑TTT‑ATTCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGa             <  1:911469/151‑1 (MQ=255)
                                                                                                                                            cAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg              >  2:1209571/1‑149 (MQ=255)
                                                                                                                                               gTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAt         >  1:1335803/1‑151 (MQ=255)
                                                                                                                                                        aCC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtttt  >  1:35671/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                       
CCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATT‑ACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537962‑4538258

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