Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,289,022 C→A 15.8% Q71K (CAA→AAA)  yraJ → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,289,0220CA15.8% 160.1 / 3.8 65Q71K (CAA→AAA) yraJputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (39/15);  new base A (2/9);  total (41/24)
Fisher's exact test for biased strand distribution p-value = 1.29e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

AAAATATAAAAACACAGGTCATCAGGGAATGCCACAACGACACCACCAGGGACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCA  >  NC_000913/3288784‑3289259
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
aaaaTATAAAAAACCTGGGCCTCAGGGAATGCCACAACGACACCACCAGGGACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGcac                                                                                                                                                                                                                                                   <  1:492841/239‑1 (MQ=255)
  aaTATAAAAACACAGGTCATCAGGGAATGCCACAACGACACCACCAGGGACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGCATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑ACAAAAATATGCAC‑aa                                                                                                                                                                                                                                                >  2:213802/1‑239 (MQ=255)
       aaaaaCACAGGTCATCAGGGAATGCCACAACGACACCACCAGGGACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTAcc                                                                                                                                                                                                                                           <  2:620754/239‑1 (MQ=255)
           acacAGGTCCTCAGGGAATGCCACAACGACCCCACCAGGGACATTAACGCACCCCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGa                                                                                                                                                                                                                                       <  1:54798/239‑1 (MQ=255)
                ggTCATCAGGGAATGCCACAACGACACCACCAGGGACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAACTATGCCC‑AACTACCAGGAACttt                                                                                                                                                                                                                                  >  2:646387/1‑238 (MQ=255)
                      caGGGAATGCCACAACGACACCACCAGGGACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCCCAAAAA‑GAGGAATGC‑AACAAACCGATCTTTCGCGTTTTTC‑AAAAAAATATGAAC‑AACTACCAGGAACTTATCAGGt                                                                                                                                                                                                                            >  2:72309/1‑239 (MQ=255)
                      caGGGAATGCCACAACGACACCACCAGGGACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGAATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAATAATCGAC‑AACTACCAGAAACTTATCAGGt                                                                                                                                                                                                                            >  2:935230/1‑239 (MQ=255)
                           aaTGCCACAACGACACCACCAGGGACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAA‑AACTACCAGGAACTTATCAGGTTGa‑ta                                                                                                                                                                                                                      >  2:264411/1‑239 (MQ=255)
                             tGCCCCAACCAACCCCCCAGGGAAATAAACGCCAATCGAAACAGTTGGGGGTCAGCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGATTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATc                                                                                                                                                                                                                    <  1:213802/239‑1 (MQ=255)
                             tGCCACAACGACACCACCAGGGACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AACCAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTCCCAGGAACTTATCAGGTTGa‑ta                                                                                                                                                                                                                      >  2:778898/1‑237 (MQ=255)
                              gCCACAACGACACCACCAGGGACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCCC‑AACTa                                                                                                                                                                                                                                             >  2:514764/1‑214 (MQ=255)
                                cacaACGACACCACCAGGGACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAATAGTATTATTTCGACCCCATTATCCTGGAAACCACAAAAA‑GTGGGATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAAAATGCAC‑AACTACCAGGAACTTTCCAGGttatatacggg                                                                                                                                                                                                                  >  1:797066/1‑230 (MQ=255)
                                cacaACGACACCACCAGGGACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGAGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATAGGAAC‑AACTACAAGGAACTTATAAGGTTGT‑AATCTgc                                                                                                                                                                                                                 >  2:943402/1‑238 (MQ=255)
                                 ccacagaaacccccAGGGACATAAACGCACTCCGAAACAGTGGGCGCTCATCATTAAACGTTGTTTGCCGAAGGGGCTCACTGGCAGCGTCATGGTTTGCACTACCGCTAACGCCGAATAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGc                                                                                                                                                                                                                <  1:684841/231‑1 (MQ=255)
                                  caaccccccccccaGGGACATAAACGCACCCCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCt                                                                                                                                                                                                               <  1:943402/228‑1 (MQ=255)
                                   aaCGACACCACCAGGGACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGcac                                                                                                                                                                                                                                                   <  1:219005/204‑1 (MQ=255)
                                   aaCGACACCACCAGGGACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGcac                                                                                                                                                                                                                                                   >  2:219005/1‑204 (MQ=255)
                                   aaCGACACCACCAGGGACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAgg                                                                                                                                                                                                                                        >  2:907270/1‑214 (MQ=255)
                                           ccccaGGGACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATaaaaa                                                                                                                                                                                                      <  1:941619/237‑1 (MQ=255)
                                                gggACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGATTAAAAAGAAgg                                                                                                                                                                                                 >  2:362315/1‑239 (MQ=255)
                                                  gACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑CAAAAAATATGCAC‑AACTACCAGGAAATTATAAAGGAAA‑AATCTAGCGGAAAAAAAAAAAGGtt                                                                                                                                                                                               >  1:429841/1‑239 (MQ=255)
                                                  gACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGCGCTCACAGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATAATTTCGACCCCATTATGCTGGAAACCACAAAAA‑AGGGTATCA‑AAACAACCGCTCTGTCCCGTTTTCC‑ATAAAAAAAAGACA‑AACACCCGGGAACTTACCAGATTAA‑TACCTgcccaataagaaaaaaggg                                                                                                                                                                                                >  1:899336/1‑220 (MQ=255)
                                                                 cGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAAAATGCAC‑AACTACCAGGAACTATTCAGGTTGC‑TATCTGGCTGAATAAAAAAAAGGTTTCAAAAAAAAAAA‑tt                                                                                                                                                                               >  1:52819/1‑239 (MQ=255)
                                                                           ggCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑Acc                                                                                                                                                                     >  2:777798/1‑238 (MQ=255)
                                                                           ggCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GGGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGT‑TATCTGGCTGAATAAAAAAAAGGGTTCACAGAAAAAAA‑TTAAAATT‑ACCg                                                                                                                                                                    >  1:944467/1‑239 (MQ=255)
                                                                                 catcatcaAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGCATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCG‑TTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGCTGA‑TATCTGGCTGAATAAAAAAAAGGTTTACCAGAAAAAAA‑TTAACTTT‑ACTGCCAATga                                                                                                                                                             >  2:365813/1‑238 (MQ=255)
                                                                                                  tGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCa                                                                                                                                              >  2:337395/1‑238 (MQ=255)
                                                                                                  tGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGCTGA‑TATCCGGCTGAATAAAACGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCAAATGCCGAGCAACTTTTGCa                                                                                                                                              >  1:523858/1‑238 (MQ=255)
                                                                                                          gTGCTCACTGGCAGCGGCATTGTTTGCACTCCCGGAAACGCCGAAGGGTATTATTTCGACCCCATTATGCTGGAAACCAAAAAAA‑GTGGTATGA‑AAACAACGGATCTGTCGCGTTTTTC‑AAAAAAATATGCCA‑AACTACCAGGAATTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTACAAAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGtt                                                                                                                                     <  2:429841/239‑1 (MQ=255)
                                                                                                                  tGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGCAGAGTATTATTTCGACCCCATTATGCTGGAAACCACCAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAg                                                                                                                             >  2:97708/1‑239 (MQ=255)
                                                                                                                     cAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCCATGCAGATCAACTTCTGCAGCCCCAATTTACgggaaaaa                                                                                                                          >  1:962291/1‑233 (MQ=255)
                                                                                                                          gCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGAAACTTCTGCAGCCCAAGCTTACGGTAGAAAAACt                                                                                                                      >  2:338782/1‑238 (MQ=255)
                                                                                                                            ttttttggccctcccGCTAAGCCCGAAGGGTTTTATTCCCCCCCCTTTTGTCTGAAAACCACAAAAA‑GTGGAATCC‑AACCAACCGATCTGCCGCGGTTTTA‑AAAAAAATATGCAA‑AACTCCCAGGAACTTATCAGGTTTA‑TATCTGGTTAAATAAAAAGAAGTTTGCAAAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTAc                                                                                                                    <  1:97708/225‑1 (MQ=255)
                                                                                                                                 ttGCACTACCTCTAACGACGAAAAGTATTATTTCGCCCCCATTATGCTGGAACCCAAAAAAA‑GTGGTATGC‑AAACAACCGTTCTGGCGCGTTTTTC‑AAAAAAAAATGCAC‑AACTACAAGGACTTTATCAGGTTGA‑TATCTGGCTGAAAAAAAAGAAGTTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAg                                                                                                               <  2:375817/238‑1 (MQ=255)
                                                                                                                                  tccACTCCCGCTAACGCAGAAGAGTTTTATTTCGCCCCAATTATGCTGAAACCCAAAAAAA‑GGGGTATGC‑AAACAACCGATCTGTCGGGTTTTTC‑AAAAAAATAGGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCt                                                                                                             <  2:48949/237‑1 (MQ=255)
                                                                                                                                  tGCCCTCCCGCAAACGCCAAAGAGTTTTATTCCGCCCCAATTAGGCTGAAACCCAAAAAAA‑GTGGAATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAGAATATGCAA‑ACCTACCAGAAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCt                                                                                                             <  2:155895/239‑1 (MQ=255)
                                                                                                                                            cTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAg                                                                                                   >  2:992353/1‑239 (MQ=255)
                                                                                                                                            cTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGAAACTTCTGAAGCCACAGTTTACGGTAGAACAAATACGGCAGCTGGGATTTAAg                                                                                                   >  2:44795/1‑239 (MQ=255)
                                                                                                                                               aaagagAAGGGTATTATTCCCACCCCATTATGCT‑GAAAACACAAAAAGGTGGTAGGC‑AAACAACCGATTTGTCGCGTTTTAA‑AAAAAAATAGGAAC‑AACTACAAGGAATTTAGCAGTTGGA‑TATTTGGCTGAAAAAAAAGAATGTTCAAAAAAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGt                                                                                                 <  2:293802/233‑1 (MQ=255)
                                                                                                                                                 gCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTA‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGACCACCTTCTGCAGCCACAGTTTACGGTAAAAAAACAACGTGAACTGGGAATAAAGGGGGa                                                                                              >  1:490410/1‑239 (MQ=255)
                                                                                                                                                            ttattTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTCCGTGAGCTGGGAATTAAGGTGGTTGAAATCCCgg                                                                                   >  2:370585/1‑239 (MQ=255)
                                                                                                                                                              agtcgcaCCCAATTACGCTGAAAACCAAAAAAA‑GGGGTATGA‑AAACAACCGTTCTGGCGGGTTTTTA‑AAAAAAATAGGCAA‑ACCTACCAGGAACTTATCAGGTTGA‑TATCTGGCTAAATAAAAAGAAGGTTT‑ACAAAAAAAAATTTACTTTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGgcg                                                                                 <  2:797066/233‑1 (MQ=255)
                                                                                                                                                                 tCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGACAAACAACGTGAGCTGGGTATTAAGGGGGATGAAAACCCGGCgctg                                                                              >  1:467896/1‑239 (MQ=255)
                                                                                                                                                                  cGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAAAAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCgctg                                                                              >  2:418738/1‑238 (MQ=255)
                                                                                                                                                                    aCCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTTTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTAACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGTTCAAACAAACACGTGAGCTGGGGCTTAAGGCGGATGAAATccctgtcccttaa                                                                           >  1:293802/1‑229 (MQ=255)
                                                                                                                                                                    aCCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCgctg                                                                              >  2:713532/1‑236 (MQ=255)
                                                                                                                                                                    aCCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCcctggct                                                                           >  2:646135/1‑239 (MQ=255)
                                                                                                                                                                    aCCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGAGAGCTGGGGATTAAGGGGGAGGAAAACCGGGCgcgggcg                                                                           >  1:997828/1‑238 (MQ=255)
                                                                                                                                                                    aCCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTAACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCCCAGTTTACGGTGGAACAACTACGTGAGTTGGGTATTAAGGGGGAGGAAAACCCgagacgggcg                                                                           >  1:682070/1‑230 (MQ=255)
                                                                                                                                                                       ccATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGACTTAAGGGGGATGAAATCCCGGGGCTGGCTGaa                                                                        >  1:509067/1‑239 (MQ=255)
                                                                                                                                                                        cATTATGCTGGAAACCAAAAAAA‑GTGTTATGC‑AAACAACCGATCGGTCGCTTTTTTA‑AAAAAAATATGCAC‑AACTACCAGGAATTTACAAGGTTGA‑TATCTGGCTGAAAAAAAAGAAGGTTTCACAAAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGaaa                                                                       <  2:509067/239‑1 (MQ=255)
                                                                                                                                                                             ttgcggAAACCCAAAAAA‑GGGGGAAGGAAAACAACCGATCGTTCGGTTTTTTC‑AAAAAAATATGACA‑AACTCCCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGTTTTCAAAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGa                                                                   <  2:52819/235‑1 (MQ=255)
                                                                                                                                                                                 ggAAACCAAAAAAA‑GGGGAAACA‑AAAAAACCGATCTGTCGCTTTTTTA‑AAAAAAAATGGCAA‑AACTACAAGAAACTTAAAGGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCAAAAAAAAAAA‑TTAAATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACg                                                              <  2:944467/239‑1 (MQ=255)
                                                                                                                                                                                  gAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGGAATCCCGGCGCTGGCTGAAAAAGATGACGa                                                             >  2:72452/1‑239 (MQ=255)
                                                                                                                                                                                  gAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGa                                                             >  2:758648/1‑239 (MQ=255)
                                                                                                                                                                                       cacaAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGAATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCg                                                        >  2:627871/1‑239 (MQ=255)
                                                                                                                                                                                         aaaaat‑GTGGTAGGA‑AAACAACCGATTTGTCGCTTTTTTC‑AAAAAAATATGCCC‑AACTACCAGGAACTTATCAGGTGTA‑TATCTGGATAAAAAAAAAAAAGGTTTCAAAGAAAAAAA‑TTACTTTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTg                                                      <  2:523858/238‑1 (MQ=255)
                                                                                                                                                                                              a‑GTTGTATGC‑AAACCCCCTGTCTTGCGCTTTTTTA‑AAAAAAAGAGTCCA‑AACCACCAGCAAATTAACAGGTTAA‑TTTCTGGCTAAATAAAAAGAAGGTTTCAAAGAAAAAAT‑TTACATTG‑ACCGCCAATGCAGAGCAGCTTCTGCAGCCACAGGTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCaa                                                 <  1:264411/239‑1 (MQ=255)
                                                                                                                                                                                                  ggTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTca                                                                                                                                                                                            >  1:470717/1‑99 (MQ=255)
                                                                                                                                                                                                  ggTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTca                                                                                                                                                                                            <  2:470717/99‑1 (MQ=255)
                                                                                                                                                                                                          aaaCAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCCCTTGac                                       >  2:399379/1‑238 (MQ=255)
                                                                                                                                                                                                           aacaccCGATCTGCCGCGTTTTAC‑AAAAAAATATGCAC‑AACTACCAGAAACTTATCAGGTTGA‑TATCTGGCTAAATAAAAAGAAGGTTTCAAAAAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAAc                                      <  2:467896/239‑1 (MQ=255)
                                                                                                                                                                                                               caccgtatctgccctttttCAAAAAAAATACGCCA‑ACCTCCAAGGAACTTATATGGTGTA‑TATTTGGCTGAAAAAAAAAAAGTTTCCACAAAAAAAAA‑TTACTTTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACaaa                                   <  2:609345/225‑1 (MQ=255)
                                                                                                                                                                                                                   gtttttttgctttttt‑AAAAAAATTAATACAAACTTCCAGAAAGTTAACAGGTTGT‑TATGTGGCTAAAAAAAAAGAAGCTTCAACAGAAAATAA‑TTACTTTT‑CCCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCAt                               <  1:362950/229‑1 (MQ=255)
                                                                                                                                                                                                                   aTCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAAt                                                                                        >  2:459428/1‑181 (MQ=255)
                                                                                                                                                                                                                   aTCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCAt                               >  2:274088/1‑238 (MQ=255)
                                                                                                                                                                                                                      tGTCGCGTTTTTA‑AAAAAAATGGGCAA‑AACTTCCAGGACCTTATCAGGTTGA‑TTTTAGGCTGAAAAAAAGGAAGGTTTCAAAAAAAAAAA‑TTACATTA‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTccc                           <  2:899336/239‑1 (MQ=255)
                                                                                                                                                                                                                           tataaaaaacaaaaaaTACGAAC‑AA‑TACCCAGAACTTTACTAGTTTATTATTTGGTTAAAAAAAAAAAAGG‑TTCACAAAAAAAAA‑TTTCCTTTAACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTa                       <  1:928614/229‑1 (MQ=255)
                                                                                                                                                                                                                                               cac‑aACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATT‑CCGGTACAGCTGCTGAATTTGATTTCa  >  2:248469/1‑239 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
AAAATATAAAAACACAGGTCATCAGGGAATGCCACAACGACACCACCAGGGACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAA‑GTGGTATGC‑AAACAACCGATCTGTCGCGTTTTTC‑AAAAAAATATGCAC‑AACTACCAGGAACTTATCAGGTTGA‑TATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAA‑TTACATTT‑ACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCA  >  NC_000913/3288784‑3289259

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