Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,624,670 C→T Q40Q (CAG→CAA ytjB ← putative inner membrane protein Smp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,624,6700CT98.6% 227.1 / ‑4.6 71Q40Q (CAG→CAAytjBputative inner membrane protein Smp
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (36/34);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-01

CCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAAC  >  NC_000913/4624425‑4624903
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                         
ccgcacgCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATAGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTAGGTACGCATCAGCGGTGCCACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCAGGGGATTACGTTGTCg                                                                                                                                                                                                                                      >  2:430795/7‑251 (MQ=255)
  gcccgCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTTCCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCt                                                                                                                                                                                                                                    <  1:301412/247‑1 (MQ=255)
        ttttttACCGTCGAGCGGCAGCCGGTCGCGCACTTCGCCGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGTATGCGATTTTCATCCGGCGAGTCGGGACGCCTCAGCGGTGAACCGTTCAGCGTTACCTGACGAGCCTAGGTGCGAGCAAGTTCTTCCAGCTGTGGTTTTCGTTGTCGCtggtggg                                                                                                                                                                                                                              >  2:263503/1‑250 (MQ=255)
             tACCGTCGGGCGCCAGCCGGGCGCGCACTTCGTCGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGAGCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGttttt                                                                                                                                                                                                                          <  2:424130/250‑2 (MQ=255)
                        gcCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGACCGGGCGATAAGATAGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCCGTTCTTCCCGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATg                                                                                                                                                                                                               >  1:375644/1‑250 (MQ=255)
                          cAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATAAGGCTTTCATACGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTAGATGGAGGTTTTGACTAAACCATGa                                                                                                                                                                                                              >  1:320636/1‑249 (MQ=255)
                              cGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTcc                                                                                                                                                                                                        <  1:439504/251‑1 (MQ=255)
                              cGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGGTTTGACTAAACCATGACGCTcc                                                                                                                                                                                                        >  2:228040/1‑251 (MQ=255)
                                    cgcACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCAt                                                                                                                                                                                                  <  1:240000/251‑1 (MQ=255)
                                       aCTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAg                                                                                                                                                                                               <  1:276184/251‑1 (MQ=255)
                                            gcagCTTTCGCCAGAACGGGCGATAAGATCGCCGTGTTCGTCAGTCACGCCCGCGTCGAGGATACGGTTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCa                                                                                                                                                                                          <  2:198195/248‑1 (MQ=255)
                                               gCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCAAGGGTGCGGGCCAGCTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGCCGCTACCTGCATCAGCGCCACt                                                                                                                                                                                        >  1:454325/1‑250 (MQ=255)
                                                cTTTCGCCAGAACGGGCGCTAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTACACCATGACGCTCCCTGCATCAGCGCCACTac                                                                                                                                                                                      >  2:314930/1‑250 (MQ=255)
                                                       ccgaacaggCGATAAGATCGCCTTTTTCCTCATCCATGCCCGCGTCGAGGATACGGCTTTCATCCGTTATCTGATCGAGGATCGCTTGTATGTGTTGTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTGTCCTGTCGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGGCGCTGGTGGTTTTGACGAAACCATGACGCTCCCTTCATCAGCGCCACTAACAAGGCa                                                                                                                                                                               <  2:171556/243‑1 (MQ=255)
                                                          aaCGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGa                                                                                                                                                                            <  1:241672/251‑1 (MQ=255)
                                                              ggCGATAAGATCGCCTTGGTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGa                                                                                                                                                                        <  1:227361/250‑1 (MQ=255)
                                                                 gATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGGTACCTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAAACATGTCGCTCCCTGCATCAGCGCACCTAACAAGGCTAGaacgaac                                                                                                                                                                      >  2:384575/1‑244 (MQ=255)
                                                                           ccTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACCATccc                                                                                                                                                           >  2:266524/1‑249 (MQ=255)
                                                                              tGTTCGTCATACACGCCCGCGGCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTTCGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACt                                                                                                                                                          <  2:341866/249‑1 (MQ=255)
                                                                                  cGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGCGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCg                                                                                                                                                     >  1:431250/1‑250 (MQ=255)
                                                                                   gTCATACACGCACGCGTCGAGGATAAGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGGGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCAAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCCTGACGCACCCTGCATCAACGCCACTAACAAGGCAAGACAGAACAGGACAATCACAGCCCGAt                                                                                                                                                   >  1:262884/1‑251 (MQ=255)
                                                                                     cATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGc                                                                                                                                                 >  1:290144/1‑251 (MQ=255)
                                                                                         cacgcccGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGcc                                                                                                                                             >  1:287289/5‑251 (MQ=255)
                                                                                                      gggATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGttt                                                                                                                                <  2:237649/250‑1 (MQ=255)
                                                                                                      aGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGttt                                                                                                                                >  1:352309/1‑251 (MQ=255)
                                                                                                       ggATACGGCGTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTGTCACCCGGCGAGTCGGTACGCAGCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGtttt                                                                                                                               <  2:431250/251‑1 (MQ=255)
                                                                                                         aTACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTtgt                                                                                                                             >  2:249962/1‑251 (MQ=255)
                                                                                                             ggCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGa                                                                                                                         <  2:352309/251‑1 (MQ=255)
                                                                                                               cTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATAACTGCCCGATGCAGCCGACATTTCAGTTTTGTGCGAGc                                                                                                                       >  1:221253/1‑251 (MQ=255)
                                                                                                                tttCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGcc                                                                                                                      >  1:247783/1‑251 (MQ=255)
                                                                                                                        gTTAACTGATCGAGGATCGCCTTAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCa                                                                                                              <  1:388097/251‑1 (MQ=255)
                                                                                                                           aaCTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGCGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGAc                                                                                                                                                                                                                        >  2:28187/1‑142 (MQ=255)
                                                                                                                           aaCTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGAc                                                                                                                                                                                                                        <  1:28187/142‑1 (MQ=255)
                                                                                                                               gATCGAGGAGCCTCTGAATGCGTTTTTCATCCGGCGAGCCGGTACGCATCAGCGGTGCAACGTTCAGCGTTCCCTGACGCGCGAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGGTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGTCGGAAGTTTAGGTTTGTGCGAGCCATCTTCCTCCCTTTg                                                                                                       <  1:405491/251‑1 (MQ=255)
                                                                                                                                   gagtgtCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCCGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGaaaa                                                                                                   <  2:433673/246‑1 (MQ=255)
                                                                                                                                   gAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGaaaa                                                                                                   <  1:191121/251‑1 (MQ=255)
                                                                                                                                    aGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCCCCCTTTGAAAAt                                                                                                  >  2:169848/1‑251 (MQ=255)
                                                                                                                                              gAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACtt                                                                                         >  2:197682/1‑250 (MQ=255)
                                                                                                                                              gAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGAAAATCAATGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCGTCTTCCACCCTTTGAAAATTTGAGCCtt                                                                                         >  2:82185/1‑250 (MQ=255)
                                                                                                                                                aTGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAAt                                                                                      <  1:249962/251‑1 (MQ=255)
                                                                                                                                                       tttCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGcc                                                                                <  1:159467/250‑1 (MQ=255)
                                                                                                                                                        ttCATCCGGCGAGTCGGGACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGTCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCTCCCGTTGAAAATTTGAGACTTAATGTTGCCa                                                                               <  1:106489/250‑1 (MQ=255)
                                                                                                                                                         tCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGcc                                                                                >  1:7914/1‑248 (MQ=255)
                                                                                                                                                              cggcgattCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTgag                                                                                             <  2:329487/223‑1 (MQ=255)
                                                                                                                                                              cggcgAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTgag                                                                                             >  1:329487/5‑230 (MQ=255)
                                                                                                                                                               ggcgAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCaa                                                                                                                                                                              >  1:407624/4‑148 (MQ=255)
                                                                                                                                                               ggcgAGTCGGGACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCaa                                                                                                                                                                              <  2:407624/145‑1 (MQ=255)
                                                                                                                                                                   aGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAg                                                                                                                        <  2:237378/198‑1 (MQ=255)
                                                                                                                                                                   aGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAg                                                                                                                        >  1:237378/1‑198 (MQ=255)
                                                                                                                                                                     gcGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCGGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACGAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATAc                                                                 <  2:437045/250‑1 (MQ=255)
                                                                                                                                                                      cGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTCACAAGGCAAGACAGAACAGGACCATCCCTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGATAATTTGAGACTTAATGTTGCCAGAAGCAATGGAtcc                                                                 >  2:125658/1‑248 (MQ=255)
                                                                                                                                                                      cGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCAGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACACGGCAAGACAGAACAGGACCATCACTGCCCGATGCAGCCGGAATTTCAGTTTCGTGCGAGCCATCTTCCACCCTTTGAAACTTTGAGACTTAATGTTGCCAGAAGCAAAGGATAcc                                                                >  2:71185/1‑250 (MQ=255)
                                                                                                                                                                          aCGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAgg                                                             <  1:108130/250‑1 (MQ=255)
                                                                                                                                                                            gccacAGCGGTGCAACGGTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAg                                                          <  2:117005/247‑1 (MQ=255)
                                                                                                                                                                            gCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGACGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAGTTTGAGACTTAATATTGACAGAAGCAATGGATACAAgga                                                            >  2:263466/1‑248 (MQ=255)
                                                                                                                                                                            gCATCAGCGGTGCAACGTTCAGCGTTACCCGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGGTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCGTGGAAAATTTGAGACTTAATGTTGCCAGAAGCAAGGGATACAAGGTAg                                                          <  2:454325/251‑1 (MQ=255)
                                                                                                                                                                                 aGCGGTGCAACGTTCAGCGTTCCCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCa                                                     <  1:116975/251‑1 (MQ=255)
                                                                                                                                                                                   cGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATg                                                   <  1:247142/251‑1 (MQ=255)
                                                                                                                                                                                    ggTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGCTTTGACTAAACCATGCCGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGATACTTAATGTTGCCAGAAGCAATGGATCAAAGGTAGCCTCATGc                                                  >  2:97224/1‑251 (MQ=255)
                                                                                                                                                                                        cAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGACGCAATGGATACAAGGTAGCCTCATGCGtt                                               >  1:250214/1‑250 (MQ=255)
                                                                                                                                                                                            gTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTAtttt                                          <  1:125658/251‑1 (MQ=255)
                                                                                                                                                                                            gTTCACCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTAtttt                                          >  1:130108/1‑251 (MQ=255)
                                                                                                                                                                                                      aCCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCg                                <  1:197682/251‑1 (MQ=255)
                                                                                                                                                                                                             gcgtctcggttcttgCCAGTTCTTCCGTCTGTGGATTAATTTTTCGCTGGTGGTTGTGACTAAACCTTGACGCTCCCCGCATCAGCTCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCGCCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTATCCTCATGCTTTATTTTCCCTGCTTCGAACGAtt                         <  2:396750/237‑1 (MQ=255)
                                                                                                                                                                                                                  aGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACa                    <  1:418253/251‑1 (MQ=255)
                                                                                                                                                                                                                         gggCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGCCTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGGGCCt             >  2:382108/1‑251 (MQ=255)
                                                                                                                                                                                                                         gggCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCt             >  1:87125/1‑251 (MQ=255)
                                                                                                                                                                                                                          agCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTACCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGAGAGCCATCTTCCCCCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTTTTTTCCCTGCTTCGAACGATTTTACAgggacctt            >  2:147676/2‑245 (MQ=255)
                                                                                                                                                                                                                                    tcCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAACTTTGAGACTTCATGTTGCCAGAAGCAATGGCTACAAGGAAGCCTCTTGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAAc  >  1:453313/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tcCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTaa   <  1:228040/250‑1 (MQ=255)
                                                                                                                                                                                                                                    tcCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAAc  >  1:68551/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                         
CCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAAC  >  NC_000913/4624425‑4624903

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