Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,328,622 T→C D496D (GAT→GAC rdcA → regulator of diguanylate cyclase RdcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,328,6220TC91.2% 177.0 / 7.3 68D496D (GAT→GACrdcAregulator of diguanylate cyclase RdcA
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base C (28/34);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-01

CTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGTT  >  NC_000913/4328376‑4328868
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
cTCGATGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGTTGAGCGGCGAGATTAAACATGAAATCGATCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACgcgc                                                                                                                                                                                                                                                     >  2:151048/1‑251 (MQ=255)
     tGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGa                                                                                                                                                                                                                                                  >  1:120570/1‑249 (MQ=255)
        aTGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGt                                                                                                                                                                                                                                              <  1:329727/250‑1 (MQ=255)
               gCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGt                                                                                                                                                                                                                                         <  2:272206/248‑1 (MQ=255)
                      gCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAgcc                                                                                                                                                                                                                                    <  1:103262/246‑1 (MQ=255)
                            cATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTc                                                                                                                                                                                                                          >  1:239422/1‑250 (MQ=255)
                            cATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTc                                                                                                                                                                                                                          >  2:478688/1‑250 (MQ=255)
                            cATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAAAAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAAAATGAAATAGAGCTGGCCCTGACCTCCACCAACCACTTTCTGCGTCAACAGCTAGACGCGCTGAAGGGACAGTAACACGCCTTGTTTc                                                                                                                                                                                                                          >  2:367009/1‑250 (MQ=255)
                              tcaaCTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCAGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAgg                                                                                                                                                                                                                       >  2:15801/3‑251 (MQ=255)
                                   tGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGgatgat                                                                                                                                                                                                                  >  2:110029/1‑251 (MQ=255)
                                    gTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGAtt                                                                                                                                                                                                                 <  1:238057/251‑1 (MQ=255)
                                     ttGAATTACGCGCAGCAGGTGCGGGAATACCTGGATTTTCTTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGAtt                                                                                                                                                                                                                 <  2:309755/250‑1 (MQ=255)
                                              gcgcAGCAGGCGCGGGAATTCCTTGAGTTTCTTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCCCCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAAGCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCTGCCTTGTTTCAGGATGATGCTGAGCCGtt                                                                                                                                                                                                       <  2:228270/251‑1 (MQ=255)
                                                gcagcaGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTa                                                                                                                                                                                                      >  1:368459/1‑250 (MQ=255)
                                                gcagcaGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTa                                                                                                                                                                                                      >  1:162038/1‑250 (MQ=255)
                                                        cgcgGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGATCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGAtt                                                                                                                                                                                             <  1:334311/251‑1 (MQ=255)
                                                        cgcgGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGGCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGAtt                                                                                                                                                                                             <  2:427745/251‑1 (MQ=255)
                                                          ggggAATTCCTGGATGTTCGGTCGCACGGGTGAAACGTCGCTTGTGAGAAAGTGCGGCAAAATAGCCACCAGATCGAAGAAAGTATGCAGTTATAGCAACTCAATCAGGCGCAGGGGAGCGGCGAGATTAAACATGAAATCGAGCTGGCTCTGACGGCCTCCAACCACGTTCTGCGTCAACTGCAAGACGCTCTGAAGGTGCAGGTAGCAGCCTTGTTACAGGAGGAGTCGGAGCCGTTGTTCGAGATTCg                                                                                                                                                                                           <  2:198653/250‑1 (MQ=255)
                                                              aaTACCTGGATTTTCGTGCGCACGGTTTAAACGTCGCTTGTGAACAAGTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTAcc                                                                                                                                                                                       <  2:102489/251‑1 (MQ=255)
                                                                                           aaCGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACACGCTGTGAGACACTGTTACAGACGgcgc                                                                                                                                                           >  1:445496/1‑250 (MQ=255)
                                                                                              ggcGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTGGCAGCTAGTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGaa                                                                                                                                                       <  2:459156/249‑1 (MQ=255)
                                                                                                  cTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGCCCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCCCAGCACAcc                                                                                                                                                   >  2:228917/1‑251 (MQ=255)
                                                                                                   ttGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACAcc                                                                                                                                                   >  1:233598/1‑250 (MQ=255)
                                                                                                        aaCAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGt                                                                                                                                             <  1:196325/251‑1 (MQ=255)
                                                                                                            aTTGCGGCAAACGATACACCAGACCGAAGAAAGTTTGCAGCTATTGCACCTCAATCAGGCGCAGGTGAGCGGCGAGCTTAAACATGAAATAGAGCTGGCCCTGCCCTCCGCCAACCAATTTCTGCGTCACCAGCAAGACGCGCTGAAGGTGCACTTCGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGCGCCGCTGTTCCAGACGGCGCCGGACACCCTCACTcgcg                                                                                                                                        >  2:27667/1‑251 (MQ=255)
                                                                                                                  gCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCGCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTa                                                                                                                                  <  1:309737/251‑1 (MQ=255)
                                                                                                                         atCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACACTCAGTCGAGACTTTACGCAGCg                                                                                                                           >  2:253266/1‑251 (MQ=255)
                                                                                                                           ccaccaGATCGAAGAAAGTTTGCAGCTAATGCAACTCAATAAGCCCCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGCCGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGtt                                                                                                                         >  1:170069/1‑251 (MQ=255)
                                                                                                                                    cgaggacaGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGAGAGATTAAACATGAACTCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCATTGCTTCAGAATGATTCTGAGCAGTTAAGCGAGATTCGTACCCACTGAGAGACCCTGTGACAGCCGTCGCAAAACACCCTCAGCCGCGCCTTTACGCTGCGTTTTGCCgcgc                                                                                                                >  2:138192/8‑248 (MQ=255)
                                                                                                                                    cGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGc                                                                                                                <  1:110029/251‑1 (MQ=255)
                                                                                                                                       agaaAGTTTGCAGATATTGCAACTCCATCCGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGACGGTGCAGTTAGCCGCCTTGTTTCAGGAGGATTCGGCGCCGTTAAGCGAGATTCGTACCAGCTGTGAGCCACTGTTACAGACGGCACAGATCACCATCGGTCTGCACTTATCGCTGCGTTTTgacaaacttg                                                                                                             >  2:195033/1‑242 (MQ=255)
                                                                                                                                           aGTTTGCAGTTATTGCAACTCAATCAGGCGCAGGTGAGCGGTGAGATTAAACATGAAATCGAGTTGGCCCTGACCTCCGCCAACCACTTTCTGCGTTAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGAtt                                                                                                                                                                                                                 <  1:243999/148‑1 (MQ=255)
                                                                                                                                           aGTTTGCAGCTATTGCAGCTCAATCAGGCGCATGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGCAGGTGCAGTTAGCCGCCTTGTTTCAGGATGAtt                                                                                                                                                                                                                 >  2:243999/1‑148 (MQ=255)
                                                                                                                                                 cAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCt                                                                                                   >  1:474619/1‑251 (MQ=255)
                                                                                                                                                   gCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTg                                                                                                                                                                        <  1:183058/181‑1 (MQ=255)
                                                                                                                                                   gCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTg                                                                                                                                                                        >  2:183058/1‑181 (MQ=255)
                                                                                                                                                   gCTATTGCAACTCAATCAGGCGCAGGGGAGCGGCGAGATTAAACAGGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGGGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACGTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCttt                                                                                                 <  2:368459/251‑1 (MQ=255)
                                                                                                                                                               cAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGtt                                                                                       <  1:369709/249‑1 (MQ=255)
                                                                                                                                                                 aTCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAAGCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGGCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCGTTTGCCGGGTGTTAAcc                                                                                   <  2:185579/251‑1 (MQ=255)
                                                                                                                                                                   cAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGa                                                                                 <  2:41152/251‑1 (MQ=255)
                                                                                                                                                                          aGGTGAGCGGCGAGATTAAACAGGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGCCGCGCTGAAGGTGCAGTTATCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACATACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTAtt                                                                          <  1:259595/251‑1 (MQ=255)
                                                                                                                                                                           gggggacggcgAGATTAAACATGAAATCGAGCTGGCCCTGGCCTCCGCCAACCACCTTCTGCGTCAACAGCAATACGCGCTGAAGGTGCAGTTAGCCGCCTTGGTTCAGGATGATGCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTc                                                                         <  2:473615/245‑1 (MQ=255)
                                                                                                                                                                            gTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACATCCGCCAACCACTTTCTGCGTCAACAGCACGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGAGCAGAACCCCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCg                                                                        >  2:208867/1‑251 (MQ=255)
                                                                                                                                                                                   gcgAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCAtt                                                                 <  2:239422/251‑1 (MQ=255)
                                                                                                                                                                                    cgAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTg                                                                <  1:295189/251‑1 (MQ=255)
                                                                                                                                                                                    cgAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGGTTTGCCGAGCTTGAATCCACCCTATGCCGGGTGTTAACCGATGTTATTCGCCCCAtgg                                                                >  1:369582/1‑249 (MQ=255)
                                                                                                                                                                                      ggaTTAAACATGAAACCGAGCGGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTTCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAg                                                              <  2:278982/250‑1 (MQ=255)
                                                                                                                                                                                                  aaaactaGCTGGCCCTGATCTCCGCCATACACTTTCTGCGTCAACATCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGTATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTGACAGACGGCGCAGAACACCATCAGTAGCGATTTTACGCTGCGTTTTGCCGAGCTTGAATCCCCCCTTTTCCGGGTGTTAACCTATGTTATTCGCCCCATTGAGCTACAAGTCaaa                                                  <  1:131874/245‑1 (MQ=255)
                                                                                                                                                                                                       gAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGa                                             <  1:428730/251‑1 (MQ=255)
                                                                                                                                                                                                        aGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGGCTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTCTTCGCCCCTTTGAGCACCGAGTCAAAATGGa                                             >  2:205757/1‑250 (MQ=255)
                                                                                                                                                                                                               ccTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGCGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGACACAATGGAATTggg                                       >  2:86640/1‑247 (MQ=255)
                                                                                                                                                                                                               ccTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGTCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCg                                     >  2:137061/1‑251 (MQ=255)
                                                                                                                                                                                                                    ccgccGCCATCCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTTAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCGTGAATCCTCCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCa                                <  1:448940/248‑1 (MQ=255)
                                                                                                                                                                                                                         gCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGtt                           >  2:168212/1‑251 (MQ=255)
                                                                                                                                                                                                                          ccAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGca                                                            >  1:13507/1‑217 (MQ=255)
                                                                                                                                                                                                                          ccAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGttt                          <  1:28857/251‑1 (MQ=255)
                                                                                                                                                                                                                          ccAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGGGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGca                                                            <  2:13507/217‑1 (MQ=255)
                                                                                                                                                                                                                             accccTTTCTGGGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGTGTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACGTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGt                       <  2:445496/247‑1 (MQ=255)
                                                                                                                                                                                                                                    tCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTggg                 <  1:478688/250‑1 (MQ=255)
                                                                                                                                                                                                                                        cGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGCTGCTTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGCGAGCCACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCACCAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTc             >  1:465286/1‑250 (MQ=255)
                                                                                                                                                                                                                                           cAACAGCAAGACGCGTTGGTGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCAtt          <  1:483100/250‑1 (MQ=255)
                                                                                                                                                                                                                                            aaCAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCAtttt        <  1:450141/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aaCAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGGGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCCAAATGGCATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCGtttt        >  2:137112/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 cAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGCGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCCTTGAGCAACAAGTCAAAATGGAATTGCGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGt   >  2:310758/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  aaGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGtt  <  2:233598/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  aaGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTgg                  >  2:291045/1‑235 (MQ=255)
                                                                                                                                                                                                                                                  aaGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTgg                  <  1:291045/235‑1 (MQ=255)
                                                                                                                                                                                                                                                  aaGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGtt  <  1:232435/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
CTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGTT  >  NC_000913/4328376‑4328868

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