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

CCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGC  >  NC_000913/4537964‑4538249
                                                                                                                                                   |                                                                                                                                               
ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCttt                                                                                                                                            >  2:134522/1‑151 (MQ=255)
  gCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑At                                                                                                                                         >  2:916868/1‑151 (MQ=255)
   cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑Att                                                                                                                                        >  2:622402/1‑151 (MQ=255)
   cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑At                                                                                                                                         <  2:139413/150‑1 (MQ=255)
    aaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑Att                                                                                                                                        <  2:28160/150‑1 (MQ=255)
    aaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTc                                                                                                                                       >  1:12801/1‑151 (MQ=255)
        ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTc‑aca                                                                                                                                   >  1:862930/1‑150 (MQ=255)
         cGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑c                                                                                                                                <  2:905627/151‑1 (MQ=255)
              gCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑ccgccg                                                                                                                           >  1:1144075/1‑151 (MQ=255)
                  ctccATCTGGTGATGAGCCGGGAGCAAGCTTATTCCATTTTAATTTATTACCGGTGAAATCAAGTTAAAATACGGCTTCCGTTCGCAATCCTGTTTTGGCTTCGCC‑TTTAATAAG‑CCGGGTTTTATCCCCTTT‑ATTC‑ACAA‑Ccgccgc                                                                                                                          <  2:1133127/146‑1 (MQ=255)
                      aTCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGcc                                                                                                                   >  1:5939/1‑151 (MQ=255)
                              atgaGACGGGAGCAAGCTTATTCCATTTTAAATTATTCCCGGTAAAATCAAGTTAAAATCCGGCATCCTTTCGCAATCCGGTTTTGGTTTCGCA‑TTTAATAAG‑CCAGGTTTTACCCCTTTT‑ATCC‑ACAA‑CCGCCGCACCAGCCGTTCCGTa                                                                                                           <  1:1144530/151‑1 (MQ=255)
                                gagaCGGGAGCAAGTTTTTTCCATTTTAAATTATTCCCGGTAAAATAAATTTAAAATACGGCACCCGTTCGCAACCCTGTTTTGGCTCCCCC‑ATAAATAAG‑CCAGTTTTTATCCCTTTT‑ATCC‑ACAA‑CCGCCGCACCAGCCGTTCCGTAc                                                                                                          <  2:834885/150‑1 (MQ=255)
                                         gCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGAtt                                                                                                >  2:409035/1‑151 (MQ=255)
                                          cAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                               >  2:862070/1‑151 (MQ=255)
                                           aaGCTTATTCTATTTAAAATTATTCCCGGTAAAATAAAGTTAAAATACGGCACCCTTCCGCAATCCTGTTTTGGCTTCGCC‑ATAAAAAGG‑CCAGGTTTTATCCCCTTT‑ATCC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                               <  2:1164386/150‑1 (MQ=255)
                                           aaGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg                                                                                              >  2:1029462/1‑151 (MQ=255)
                                             gtttATTCCATTTTAAATTATTCCGGGTGAAATCAAGTTCAAACACGGCATCCGTCCGCAATCCTGTTTGGGCTTCGCC‑ATTAATAAG‑CCAGTTTTTACCCCCTTT‑ATCC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc                                                                                             <  1:11047/148‑1 (MQ=255)
                                              cTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                            >  1:333356/1‑150 (MQ=255)
                                               ttATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                            >  2:595635/1‑149 (MQ=255)
                                                tATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGc                                                                                         >  1:780473/1‑151 (MQ=255)
                                                    ccATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt                                                                                       >  1:745628/1‑149 (MQ=255)
                                                     cATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAAc                                                                                    >  2:921329/1‑151 (MQ=255)
                                                       ttttAAATTATTCCCGGTAAATTCAAGTTAAAACACGGCATCCGTTCGCATTTCTGGTTTGGTTTCCCA‑ATTAATAGC‑CCAGGTTTTACCCCTTTT‑ATTC‑AAAA‑CCGCCGCGCCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt                                                                                   <  2:89139/150‑1 (MQ=255)
                                                       tttaAAATTATTACCGGTAAAATAAAGTTAAAATACGGCATCCGTTCGCAACCCGGTTTGGGTTTCGCC‑ATAAATAAG‑CCAGGTTTTATCCCTTTT‑ATCC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa                                                                                      <  2:40384/147‑1 (MQ=255)
                                                         ttAAATTATTACGGGTGAAACAAATTTAAAATCCGCCATCCGTCCGCAATCCTGTTTTGGCTCCGCC‑ATAAAAAAG‑CCAGTTTTTACCCCCTTT‑ATCC‑CCAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                 <  2:436830/150‑1 (MQ=255)
                                                         ttAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTCCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                 >  1:755577/1‑150 (MQ=255)
                                                           aaaTTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                              >  2:641321/1‑151 (MQ=255)
                                                                    cgggtaaaacaaaGTTAAAATAGGGCACCCTTTCGCAACCCGGTTTGGGCTTCGCC‑ATAAAAAAG‑CCAGTTTTTACCCCCTTT‑ATCC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt                                                                      <  1:231934/139‑1 (MQ=255)
                                                                    ccGGTGAAATCAATTTAAAATACGGCATCCTTTCGCAACCCTGTTTTGGCTTCGCC‑ATAAAAAAG‑CCGGTTTTTATCCCCTTT‑ATCC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt                                                                       <  1:827551/149‑1 (MQ=255)
                                                                    ccGGTAAAATCAATTTCAATTCCGGCACCCTTTCCCAATCCTGGTTTGGTTCCCCC‑ATTAATAAC‑CCGGTTTTTATCCCCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt                                                                       <  1:231333/149‑1 (MQ=255)
                                                                        taaaaTCAAGTCAAATTACGGCATCCGTCCTCAATCCTTTTTTTGTTTCGCC‑TTTAAAAAG‑CCAGGTTTTATCCCTTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAgg                                                                  <  1:464644/148‑1 (MQ=255)
                                                                          aaaaCAAGTTAAATTACGGCATCCGTTCGCAATCCGGTTTTGGCTTCGCC‑ATTAAAAAG‑CCAGGTTTTATCCCCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGa                                                                <  1:947794/150‑1 (MQ=255)
                                                                          aaaTCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt                                                               >  1:870448/1‑151 (MQ=255)
                                                                           aaTCAAGTTCAAATACGGCATCCGTTCACAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc                                                              >  1:377848/1‑151 (MQ=255)
                                                                              cAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                            >  2:1096687/1‑150 (MQ=255)
                                                                                 gTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                           >  1:155079/1‑148 (MQ=255)
                                                                                  ttCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                       >  1:187387/1‑151 (MQ=255)
                                                                                    cAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGaca                                                      >  2:897122/1‑150 (MQ=255)
                                                                                     aaatcgggCACCGGTTCGCAATCCTGTTTGGGCTTCGCA‑ATTAAAAGG‑CCAGTTTTTACCCCCTTT‑ATCC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                       <  1:1064024/142‑1 (MQ=255)
                                                                                        tACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                    >  1:98620/1‑148 (MQ=255)
                                                                                          cGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa                                                >  2:1018061/1‑150 (MQ=255)
                                                                                           ggCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAAc                                              >  1:795050/1‑151 (MQ=255)
                                                                                           ggCACCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCCCCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa                                                <  1:910769/149‑1 (MQ=255)
                                                                                                ccGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                         >  1:422868/1‑151 (MQ=255)
                                                                                                ccGTTCGAAACCCTGTTTTGGTTTCGCC‑TTAAAAAAG‑CCAGTTTTTACCCCCTTT‑ATCC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaaa                                               <  1:223890/145‑1 (MQ=255)
                                                                                                    tCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                       >  2:243179/1‑149 (MQ=255)
                                                                                                      gCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                    >  1:522473/1‑150 (MQ=255)
                                                                                                      cgaaccccggTTTGGGTTCCCC‑AAAAAAAAACCCAGTTTTTACCCCGTTT‑ATCC‑AAAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa                                                <  2:860425/131‑1 (MQ=255)
                                                                                                      cccaTCCCGGTTTGGCCTCCCCTATTAAAAAG‑CCTGTTTTTATCCCCTTT‑ATCC‑AAAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                       <  2:264305/145‑1 (MQ=255)
                                                                                                        atccttggtttggttTCCCC‑ATAAAAAAC‑CCGGGTTTTACCCCTTTT‑TTCC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                    <  2:337366/143‑1 (MQ=255)
                                                                                                         acccTGGTTGGGCAACACC‑TTTAATCGG‑CCATGTTTTATCCCTTTT‑ATCC‑ATAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                 <  1:191228/148‑1 (MQ=255)
                                                                                                          tcttgttttGCCTCCCCC‑ATTAAAAAG‑CCGTGTTTAGCCCCTTTT‑ATTC‑ACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                       <  2:351065/137‑1 (MQ=255)
                                                                                                           ccTGTTTTGTTTTCGCC‑ATAAATAAC‑CCAGTTTTTACCCCTTTT‑ATCC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaata                                <  1:853167/149‑1 (MQ=255)
                                                                                                           ccTGTTTTGGCTTCCCC‑ATAAAAAAG‑CCAGTTTTTACCCCCTTT‑TTCC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                              <  2:165160/151‑1 (MQ=255)
                                                                                                           ccTGTTTGGGCTCCCCC‑TTTAAAAGG‑CCGGTTTTTACCCCCTTT‑TTTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaata                                <  1:354917/149‑1 (MQ=255)
                                                                                                             tggtttggttCCGCA‑ATAAAAAAG‑CAGGTTTTTACCCCCTTA‑ATCC‑ACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                              <  1:914227/149‑1 (MQ=255)
                                                                                                                tttGGTTTCGCC‑ATTAATAAC‑CCGGGTTTTATCCCCTTT‑ATCC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaag                                     <  1:895137/139‑1 (MQ=255)
                                                                                                                tttGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                         >  1:84438/1‑151 (MQ=255)
                                                                                                                   ggCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCt                      >  1:440263/1‑151 (MQ=255)
                                                                                                                   ggCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCt                      >  2:466942/1‑151 (MQ=255)
                                                                                                                   ggCTCCGCC‑ATTAAAAAC‑CCAGTTTTTATCCCCTTT‑ATCC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataa                               <  2:969522/142‑1 (MQ=255)
                                                                                                                    gtatcctcattaaataaA‑CCAGGTTTTATCA‑CTTT‑ATTCAAAAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                        <  1:133335/137‑1 (MQ=255)
                                                                                                                     gttccccctttAATAAG‑CCAGGTTTTATCCCTTTT‑ATCC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCt                      <  2:397389/141‑1 (MQ=255)
                                                                                                                      tccccc‑ATTAAAAAC‑CCAGTTTTTACCCCTTTT‑TTCC‑AAAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtt                     <  2:1144046/145‑1 (MQ=255)
                                                                                                                       tCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttttt                  >  2:358053/1‑151 (MQ=255)
                                                                                                                        cccc‑ATTAATAAG‑CCAGGTTTTATCCCCTTT‑ATCC‑AAAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt                   <  1:1160103/147‑1 (MQ=255)
                                                                                                                        cGCC‑ATAAAAAAG‑CCGGTTTTTATCCCTTTT‑TTCC‑AAAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                 <  2:413100/151‑1 (MQ=255)
                                                                                                                         gCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCAATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg                >  2:241912/1‑151 (MQ=255)
                                                                                                                         cgc‑ATTAAAAAA‑GCCGGTTTTTTCCCCTTTAATTC‑AAAT‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                 <  1:86164/149‑1 (MQ=255)
                                                                                                                           c‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc              >  1:999616/1‑151 (MQ=255)
                                                                                                                                aataaG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa           >  2:700436/1‑150 (MQ=255)
                                                                                                                                 aagaG‑CCGGGTTTTTTCCCCCTT‑CTTC‑CCAATCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa            <  2:143424/146‑1 (MQ=255)
                                                                                                                                    aG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt      >  2:50093/1‑151 (MQ=255)
                                                                                                                                    aG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt      >  1:334117/1‑151 (MQ=255)
                                                                                                                                       ccAGGTTTTTCCCCTTTT‑ATCC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg     <  2:78667/150‑1 (MQ=255)
                                                                                                                                        cAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg     >  1:1043660/1‑149 (MQ=255)
                                                                                                                                         aGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc  >  2:602612/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAG‑CCAGGTTTTATCACCTTT‑ATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGC  >  NC_000913/4537964‑4538249

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: ‑