Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 89,263 A→G intergenic (+231/‑371) cra → / → mraZ DNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091389,2630AG97.9% 146.5 / ‑4.2 47intergenic (+231/‑371)cra/mraZDNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (26/20);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

TCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTAACAAGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTTCGCACGTTGGACGTAAAATAAACA  >  NC_000913/89019‑89495
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                         
ccAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGTGCGTAGAGTACGGTAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACttt                                                                                                                                                                                                                                     <  2:107407/250‑1 (MQ=255)
  aGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTc                                                                                                                                                                                                                                   >  1:418095/1‑251 (MQ=255)
      gTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGTATTTCCCGCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTg                                                                                                                                                                                                                               <  2:376287/251‑1 (MQ=255)
         gCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCtt                                                                                                                                                                                                                            >  1:299418/1‑251 (MQ=255)
                          aaTACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGTCTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTgg                                                                                                                                                                                                           <  1:373756/251‑1 (MQ=255)
                             aCGCGCCAGGTGAATTCCCCTCTTGCTCGTAGAGTACGGGAAGTGACATCAATATCCTTAAAGTAAAAAAATCTATTCCTTACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCAt                                                                                                                                                                                                        <  2:138848/251‑1 (MQ=255)
                                        gAATTTCCCTCTGGCGGGTAGAGTGCGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTAGTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTGTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCcaca                                                                                                                                                                                             <  2:139863/251‑1 (MQ=255)
                                                ctctGGCGCGTAGAGTACGGGAATGGACCTCCATATCCTTAAAGTAAATAAGACTATTCCTGACTCTTATTAATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCACGCGCCACACAGACgc                                                                                                                                                                                      >  2:332361/1‑249 (MQ=255)
                                                         gTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACt                                                                                                                                                                           >  1:464507/1‑251 (MQ=255)
                                                             agTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGtt                                                                                                                                                                        >  1:363228/1‑250 (MQ=255)
                                                              gTACGGGACGTGACAACAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATGAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTAACAAGGACTGTTaa                                                                                                                                                                      <  2:422246/251‑1 (MQ=255)
                                                                          acctCAATATCCTTAAAGTAAATAAGACTATTCCTGACTCTTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGCTTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCTCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGa                                                                                                                                                               >  2:78999/4‑246 (MQ=255)
                                                                                   tccTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAAccc                                                                                                                                                                                                                <  2:655908/187‑1 (MQ=255)
                                                                                   tccTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAAccc                                                                                                                                                                                                                >  1:655908/3‑189 (MQ=255)
                                                                                     cTTAAAGTAAATAAGACTATTCCTGAATATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGCATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGAAGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAAt                                                                                                                                               >  2:82690/1‑251 (MQ=255)
                                                                                               aTAAGACTATTCCTGTATATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGGTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCGTTGCTTTTGCGCTAACCCTCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAAGGCCGCTcgcg                                                                                                                                     <  2:661336/251‑1 (MQ=255)
                                                                                                  aGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTACCAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACCGGGAATGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCg                                                                                                                                  >  2:669736/1‑251 (MQ=255)
                                                                                                          ccTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTACTTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTCCATGCAGAGATGAATTGCCGATTTAACAAGAACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCACAAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGACGCTCGCGTCgaagactg                                                                                                                           >  2:359670/1‑243 (MQ=255)
                                                                                                          ccTGACTATTATTGATAAGTTCTTTTAAACACGCCCGTTAATTAACTCCCCAGCTGAAATTCACAATAATTAAGTGATATCGACAGAGAGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGACCTAAAATGCCGCTCGCGTCGCAAACTg                                                                                                                           >  1:596835/1‑250 (MQ=255)
                                                                                                            tGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGACACACAGACGAAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGac                                                                                                                         >  1:605733/1‑250 (MQ=255)
                                                                                                              aCTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACt                                                                                                                      <  1:134990/251‑1 (MQ=255)
                                                                                                              aCTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCGGAAATTCACAATAATTAAGTGATATCGACAGCTCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCGAAAATGCCGCTCGCGTCGCAAACTGACACt                                                                                                                      <  2:462260/251‑1 (MQ=255)
                                                                                                                      tGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCAACAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAACTGCCGCTCGCGTCGCAAACTGACACTTTAAAttt                                                                                                              >  2:328901/1‑251 (MQ=255)
                                                                                                                                     aaCCCGGCCGTTAATTAACCAACCAGCTGAAATTCACAATAATTAAGTGATATCGCCAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGt                                                                                               <  2:194244/251‑1 (MQ=255)
                                                                                                                                        ccgcccgTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAg                                                                                            <  1:461080/251‑1 (MQ=255)
                                                                                                                                                               cTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACg                                                                      >  1:161124/1‑250 (MQ=255)
                                                                                                                                                                   aaTTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAg                                                                 <  1:61732/251‑1 (MQ=255)
                                                                                                                                                                       cacaATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGAt                                                             >  1:282088/1‑251 (MQ=255)
                                                                                                                                                                          aataatTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGTTGAGGGAtttt                                                          >  2:360482/1‑251 (MQ=255)
                                                                                                                                                                           ataatTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGAttttt                                                         <  2:161124/251‑1 (MQ=255)
                                                                                                                                                                            taatTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACCAGCACTTTTCTTTGCTTTTCCGCAAACCCGCTGGCATCAAGCGACACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTCAAACGGAGATGACGATGAGGGAtttttt                                                        >  1:677646/1‑251 (MQ=255)
                                                                                                                                                                                   gTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTg                                                                                                                           <  2:431518/177‑1 (MQ=255)
                                                                                                                                                                                   gTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTg                                                                                                                           >  1:431518/1‑177 (MQ=255)
                                                                                                                                                                                        atCGACAGCTCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCCTTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTAt                                             <  2:418095/250‑1 (MQ=255)
                                                                                                                                                                                           gACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCAt                                         >  1:367687/1‑251 (MQ=255)
                                                                                                                                                                                             cAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATaa                                       <  2:510380/251‑1 (MQ=255)
                                                                                                                                                                                                  cgTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGTTGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGtt                                   >  2:432902/1‑250 (MQ=255)
                                                                                                                                                                                                   gTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGAAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCAGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTa                                  >  1:336939/1‑250 (MQ=255)
                                                                                                                                                                                                                 ttattacctgcagagATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCACCAGGGACTGTTAACCGGGGCTGATATGTCCTAAAATGCCGCTCGCGACGCAAACTGCCACTTTCTATTTGCTGAGGATAATCGCGAGTCCTTTTAAACCGGCGCTGACGATGTGGGATTTTTTCTTACAGCTCTTCATAACGTTAATTTACTGCGCACGt                   >  2:277474/15‑251 (MQ=255)
                                                                                                                                                                                                                     tACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTgg                <  2:683494/250‑1 (MQ=255)
                                                                                                                                                                                                                       catgcagagATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACg             >  1:160753/9‑251 (MQ=255)
                                                                                                                                                                                                                        atgcagagATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGt                                                                                               >  2:629106/8‑168 (MQ=255)
                                                                                                                                                                                                                        atgcagagATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGt                                                                                               <  1:629106/161‑1 (MQ=255)
                                                                                                                                                                                                                          gcagagATGAATTGCCGATTGAATAAGCACGTTTCTTTGCTTGTGCGCAAACCAGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTGAACCGGGGAAGATATGTCCTAAAATTCCGCTCGCGTCGCAAACTGACTCTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGGTAATTTGCTGCGCCCGTTGGACGTa           <  1:445714/245‑1 (MQ=255)
                                                                                                                                                                                                                            agagATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAACACGGTGATGACGATGAGGGCTTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTaaa         >  2:350930/4‑250 (MQ=255)
                                                                                                                                                                                                                             gagATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAAt       <  2:299418/249‑1 (MQ=255)
                                                                                                                                                                                                                                   aaTTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAaaca  >  1:592363/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                         
TCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTAACAAGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTTCGCACGTTGGACGTAAAATAAACA  >  NC_000913/89019‑89495

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