Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,795,336 A→C I136V (ATA→GTC waaF → ADP‑heptose‑‑LPS heptosyltransferase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,795,3360AC98.5% 229.8 / ‑5.6 69I136V (ATA→GTCwaaFADP‑heptose‑‑LPS heptosyltransferase 2
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (31/36);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

CTATGCCTCTCGGTCACGGAGCGCTGGAAATCGGCGAACGCCGCAAACTGGGTCATAGCCTGCGTGAAAAGCGCTACGACCGCGCCTACGTCTTACCCAACTCCTTCAAATCTGCATTAGTGCCTTTCTTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGAGTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGG  >  NC_000913/3795095‑3795574
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                              
cTATGCCTCTCGGTCACGGAGCGCTGGAAATCGGCGAACGCCGCAAACTGGGTCATAGCCTGCGTGAAAAGCGCTACGACCGCGCCTACGTCTTACCAAACTCCTTCAAATCTGCATTAGTGCCTTTCTTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGcc                                                                                                                                                                                                                                       >  1:118059/1‑251 (MQ=255)
  aTGCCTCTCGGTCACGGAGCGCTGGAAATCGGCGAACGCCGCAAACTGGGTCATAGCCTGCGTGAAAAGCGCTACGACCGCGCCTACGTCTTACCAAACTCCTTCAAATCTGCATTAGTGCCTTTCTTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGcc                                                                                                                                                                                                                                       <  1:93161/249‑1 (MQ=255)
    gCCTCTCGGTCACGGAGGGATGGAAAGCGGCGACCGCCGCAAACTGGGTCATAGCCTGCGGGAAAAGCGCTACGACCGCGCCTACGTCTTACCACACTCCTTCAAATCTGCATGAGTGCCTGTCTGCGCGGGTATTCCTCAGCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGACCGCTATGTCGCGCTGGCCTAt                                                                                                                                                                                                                                    <  2:4935/250‑1 (MQ=255)
          cGGTCACGGAGCGCTGGAAATCGGCGAACGCCGCAAACTGGGTCATAGCCTGCGTGAAAAGCGCTACGACCGCGCCTACGTCTTACCAAACTCCTTCAAATCTGCATTAGTGCCTTTCTTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACaa                                                                                                                                                                                                                               >  1:16873/1‑249 (MQ=255)
                cGGAGCGCTGGAAATCGGCGAACGCCGCAAACTGGGTCATAGCCTGCGTGAAAAGCGCTACGACCGCGCCTACGTCTTACCAAACTCCTTCAAATCTGCATTAGTGCCTTTCTTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCAtt                                                                                                                                                                                                                        >  1:411531/1‑250 (MQ=255)
                 ggAGCGCTGGAAATCGGCGAACGCCGCAAACTGGGTCATAGCCTGCGTGAAAAGCGCTACGACCGCGCCTACGTCTTACCAAACTCCTTCAAATCTGCATTAGTGCCTTTCTTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTAt                                                                                                                                                                                                                      <  2:436866/251‑1 (MQ=255)
                                    aaCGCCGCAAACTGGTTCATAGCCTGCGTGAAAAGCGCTACGGCCGCGCCTACGTCTTACCAAACTCCTTCAAATCTGCATTAGTGCCTTTCTTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTg                                                                                                                                                                                                   <  2:97619/251‑1 (MQ=255)
                                             aacgtggccATAGCCTGCGTGAAAAGCGCTACGACCGCGCCTACGTCTTACCAAACTCCTTCAAATCTGCATTAGTGCCTTTCGTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTTCGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCa                                                                                                                                                                                          <  2:245662/243‑1 (MQ=255)
                                                         gcttGCGTGAAAAGCGCTACGGCCGCGCCTACGTCTTACCAAACTCCTTCAAATCTGCATGAGTGCCTTTCTTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTAGGACAAAGGCAGTATGCGTCCCGCACAAGATCTGCCGCAGCCATTGTTTTGGCCg                                                                                                                                                                              <  2:54079/248‑1 (MQ=255)
                                                          ccTGCGTGAAAAGCGCTACGACCGCGCCTACGTCTTACCAAACTCCTTCAAATCTGCATTAGTGCCTTTCTTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTCTGTCGCGCTGGCCTATGACAAAGGCATTATGAGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGc                                                                                                                                                                             >  1:214657/1‑251 (MQ=255)
                                                                tGAAAAGCGCTACGACCGCGCCTACGTCTTACCAAACTCCTTCAAATCTGCATTAGTGCCTTTCTTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGc                                                                                                                                                                       >  1:153735/1‑251 (MQ=255)
                                                                     aGCGCTACGACCGCGCCTACGTCTTACCAAACTCCTTCAAATCTGCATTAGTGCCTTTCTTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATTGCCGCAGTTGCAGGTg                                                                                                                                                                  >  1:333793/1‑251 (MQ=255)
                                                                        gcTACGACCGCGCCTACGTCTTACCAAACTCCTTCAAATCTGCATTAGTGCCTTTCTTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGc                                                                                                                                                               <  1:98043/251‑1 (MQ=255)
                                                                             gACCGCGCCTACGTCTTACCAAACTCCTTCAAATCTGCATTAGTGCCTTTCTTCGCGGGTATTCCTCATCGCCCCGGCTGGCGCGGCGAGATGCGCTACGGTTTCCTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAgg                                                                                                                                                          <  1:329633/251‑1 (MQ=255)
                                                                              aCCGCGCCTACGTCTTACCAAACTCCTTCAAAACTGCATTAGTGCCTTTCTTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACGCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGt                                                                                                                                                         <  1:17416/251‑1 (MQ=255)
                                                                                 acgcCTAAGTCTTACCACACTCCTTCAAATCTGCATTAGTGCCTTTCTTCGCGGGTATTCCTCATCGCCCCGGCTGGCTCGGCGCGATGCGCTCCGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTAGCGCTTGCCTATGACACAGGCATTATGCGTACCGCACCAGATCTGCCTCAGCCATTGTTATGGACGAAGGTGCAGGTGAGCGCAGGTgc                                                                                                                                                       >  2:281396/2‑249 (MQ=255)
                                                                                       aCGTCTTACCAAACTCCTTCAAATCTGCATTAGTGCCTTTCTTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACACAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAACAATc                                                                                                                                                 >  2:109066/1‑250 (MQ=255)
                                                                                          tCTTACCAAACTCCTTCAAATCTGCATTAGTGCCTTTCTTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCatat                                                                                                                                             >  2:313095/1‑251 (MQ=255)
                                                                                             tACCAAACGCCTTCAAATCTGCATTAGTGCCTTTCTTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGACGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCTAAGGTGAAAAATCAtata                                                                                                                                            <  1:73252/249‑1 (MQ=255)
                                                                                              accaaaCTCCGTCAAATCAGCATTAGTGCCTTGCTTCGCGGGTATTCCTCAGCGCCCCGGCTTGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCTTGCTCGATAAAGAAGCCTGGCCGCTAATGGGGGAACGCTATGTCGCGCTGGCCTATTACAAAGGCATGATGCTTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCt                                                                                                                                         <  1:268485/247‑1 (MQ=255)
                                                                                                  aaCTCCTTCAAATCGGCAGTAGTGCCTTTCGTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCGGGCCTATGACAAAGGCATTATGCGTACCGCACAAGAGCTGCCGCAGCCATTGTTATGGCCGCAGGTGCAGGTGAGCGAAGGTGAAAAAGCATAGACCTGTaa                                                                                                                                     <  2:417410/251‑1 (MQ=255)
                                                                                                   aCTCCTTCAAATCTGCATTAGTGCCTTTCTTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGGGAGCGAAGGTGAAAAATCATATACCTGTAAt                                                                                                                                    >  2:109710/1‑251 (MQ=255)
                                                                                                              tCTGCATTAGTGCCTTTCTTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCt                                                                                                                         >  2:384137/1‑251 (MQ=255)
                                                                                                              tCTGCATTAGTGCCTTTCTTAGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGAACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCt                                                                                                                         >  2:271173/1‑251 (MQ=255)
                                                                                                                 gCATTAGTGCCTTTCTTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtc                                                                                                                      >  2:106513/1‑251 (MQ=255)
                                                                                                                   aTTAGTGCCTTTCTTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAACGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtcat                                                                                                                    >  1:220690/1‑251 (MQ=255)
                                                                                                                      aGTGCCTTTCTTCGCGGGTATTCCTCAGCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtcatca                                                                                                                  <  2:333793/250‑1 (MQ=255)
                                                                                                                         gCCTTTCTTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGaa                                                                                                               <  2:220569/250‑1 (MQ=255)
                                                                                                                         gCCTTTCTTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAAc                                                                                                              >  1:247075/1‑251 (MQ=255)
                                                                                                                             ttcttcGCGGGGATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTCTGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGAAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTcc                                                                                                          >  2:176747/1‑251 (MQ=255)
                                                                                                                               cttcGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCg                                                                                                         >  1:299368/1‑250 (MQ=255)
                                                                                                                               cctcGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTGCATCAGAACGTCCg                                                                                                         <  2:411531/248‑1 (MQ=255)
                                                                                                                                 tcGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCg                                                                                                         <  1:135663/248‑1 (MQ=255)
                                                                                                                                     gggTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTCCTCAACGATGTACGCGTGCTCGATAAAGACGCCTGGCCGCTACTGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTCCCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAACTAATATACATGTCATCACTTTTCGCTTTCATCAGAACGTCCGATGATTg                                                                                                  >  1:46438/1‑251 (MQ=255)
                                                                                                                                            ccTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAATCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGCTTTTg                                                                                           <  1:112738/251‑1 (MQ=255)
                                                                                                                                                    cACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGCTTTTGCCCGGGTg                                                                                   >  1:61997/1‑251 (MQ=255)
                                                                                                                                                            ggCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGCTTTTGCCCGGGTGCGGAGt                                                                             <  2:17583/249‑1 (MQ=255)
                                                                                                                                                                  gcgAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGCTTTTGCCCGGGTGCGGAGTTTGGTCCg                                                                     <  1:384137/251‑1 (MQ=255)
                                                                                                                                                                      gaTGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGCTTTTGCCCGGGTGCGGAGTTTGGTCCGGCa                                                                  >  2:256423/1‑250 (MQ=255)
                                                                                                                                                                       aTGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAACGATGCCTGGCCGCTAATGGTGGAAAGCTATGTCGCGCTGGCCTATGCCAAAGGCATTATGCGTACCGCACACGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGCGCGAAGGAGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGCACGTCCGATGATTGGCTTTTGCCCGGGTGCGGAGTTTGGTCCGGcca                                                                 >  1:420421/1‑248 (MQ=255)
                                                                                                                                                                       aTGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAATCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGCTTTTGCCCGGGTGCGGAGTTTGGTCCGGCaa                                                                 >  1:123211/1‑250 (MQ=255)
                                                                                                                                                                              aCGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGCTTTTGCCCGGGTGCGGAGTTTGGTCCGGCAAAACGCTg                                                          <  2:117406/250‑1 (MQ=255)
                                                                                                                                                                                ggTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAg                                                                                                                                                                >  1:419405/1‑146 (MQ=255)
                                                                                                                                                                                ggTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAg                                                                                                                                                                <  2:419405/146‑1 (MQ=255)
                                                                                                                                                                                ggTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGCTTTTGCCCGGGTGCGGAGTTTGGTCCGGCAAAACGCTGGcc                                                       >  1:106691/1‑251 (MQ=255)
                                                                                                                                                                                     aCTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGCTTTTGCCCGGGTGCGGAGTTTGGTCCGGCAAAACGCTGGCCACACt                                                  <  1:361470/251‑1 (MQ=255)
                                                                                                                                                                                     aCTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCAGCAGAACGTCCGATGATTGGCTTTTGCCCGGGTGCGGAGTTTGGTCCGGCAAAACGCTGGCCACACt                                                  <  2:375442/251‑1 (MQ=255)
                                                                                                                                                                                           cGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGCTTTTGCCCGGGTGCGGAGTTTGGTCCGGCAAAACGCTGGCCACACTACCACt                                            <  1:194710/251‑1 (MQ=255)
                                                                                                                                                                                            gATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGACAAATCATATACCTGTCATCAATTTTCGCTTTCATCAGACCGTCCGATGATTGGCTTTTGCCCGGGTGCGGAGTTTGGCCCGGCAAACCGCTGGCCACACTCCCACtt                                           >  2:157634/1‑250 (MQ=255)
                                                                                                                                                                                             atgtccgcgtgatgtaTAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGAGCGGGCCTATGACAAAGGCATTAGGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCACAGTTGCAGGTTAGCGAAGGTGAAAAATCATATACCGGTAATCAATTTTCGTTTGCATCAGAACGTCCGATGATTGGCGTTGGCCCGGGTGCGGATTTTGGTCCGGCAAAACGCTGGCCACACGAccccttt                                          <  2:306783/236‑6 (MQ=255)
                                                                                                                                                                                               gTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAAGGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTTAGCGAAGGTGAAAAATCATATACCTGTAAGCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGCTTTTGCCCGGGTGCGGAGTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATg                                         <  1:176747/250‑1 (MQ=255)
                                                                                                                                                                                                 acgcaagCTCGATCAAGAAGCCTGGCCGCGCATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCGTTGTTATGGCCGAAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGCTTTTGCCCGGGTGCGGAGTTTGGTACGGAAAAACGCTGGCCACAATACCACtcagcgg                                      >  2:336361/7‑245 (MQ=255)
                                                                                                                                                                                                 aCGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGCTTTTGCCCGGGTGCGGAGTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATGCgg                                      >  1:369772/1‑251 (MQ=255)
                                                                                                                                                                                                          cGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGCTTTTGCCCGGGTGCGGAGTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCa                              >  1:87539/1‑250 (MQ=255)
                                                                                                                                                                                                           gATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGCTTTTGCCCGGGTGCGGAGTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCaaa                            <  1:313095/251‑1 (MQ=255)
                                                                                                                                                                                                           gATAAAGAAGCCTGGCAGCTAATGGTGGAACGCTCTGTGGCGCGGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCGGCAGTTGCAGGTGAGCGAAGGTGAAAAATCAGATACCTGTAATCAATTTGCTCTTTCATCAGCACGGGCGATGACTGGCTTTTGCCCGGGTGCGGAGGTTGGTCCGGCAAAACGCTGGGCACGCGGCCACTATGCTGAGCTGGCaaa                            <  2:433206/251‑1 (MQ=255)
                                                                                                                                                                                                                 gaagCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGCTTTTGCCCGGGTGCGGAGTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCt                      >  2:68298/1‑251 (MQ=255)
                                                                                                                                                                                                                 gaagCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGCTTTTGCCCGGGTGCGGAGTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCt                      <  1:274038/251‑1 (MQ=255)
                                                                                                                                                                                                                  aagCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGCTTTTGCCCGGGTGCGGAGTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCtg                     <  1:106513/251‑1 (MQ=255)
                                                                                                                                                                                                                    gCCTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAATATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGGGAAAAATCATATACCTATAATCAATTTTCGCTTTCATCATAACGTCCGATGATTGGCTTTTGCCCGGGTGCGGAGTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATTCGGATCTGGCAAAGCTGCtgat                   <  2:197916/251‑1 (MQ=255)
                                                                                                                                                                                                                      cTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACACAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGaaa                                                                                                                                                     >  1:254973/1‑119 (MQ=255)
                                                                                                                                                                                                                      cTGGCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGTTGATCGAAGTTGaaa                                                                                                                                                     <  2:254973/119‑1 (MQ=255)
                                                                                                                                                                                                                        ggCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACg                                                                                                             <  1:97314/157‑1 (MQ=255)
                                                                                                                                                                                                                        ggCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACg                                                                                                             >  2:97314/1‑157 (MQ=255)
                                                                                                                                                                                                                        ggCCGCTAATGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGCTTTTGCCCGGGTGCGGAGTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCtgattgat               <  1:424770/251‑1 (MQ=255)
                                                                                                                                                                                                                                aTGGTGGAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGCGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGCTTTTGCCCGGGTGCGGAGTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGCAGGtt        >  2:311623/1‑250 (MQ=255)
                                                                                                                                                                                                                                aTGGTGGAACGCTATGTCGCGCTGGCCTATAACAAATGCTTTATTCGGACCGCAAAAGATCTGCCGCAGCCATTGTGATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGATCGTCCGAGGATTGGCTTTTGCCCGGTTTCGGAGTGTGGTCCGGCAATACGCTGTCCACACTACCACTATGCGGAGCTGTCAAAGCAGCTGACTGATGACGGTTa       <  2:420421/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ggAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTATGCGTACCGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGCTTTTGCCCGGGTGCGGAGTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAgg  <  1:30332/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ggAACGCTATGTCGCGCTGGCCTATGACAAAGGCATTAGGCGTACCGCACAAGATCTGCCGCAGCCATTGTTAGGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTCATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGCTTTTGCCCGGGTGCGGAGTTTGGTCCGGCAAGACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAg   <  2:373522/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                              
CTATGCCTCTCGGTCACGGAGCGCTGGAAATCGGCGAACGCCGCAAACTGGGTCATAGCCTGCGTGAAAAGCGCTACGACCGCGCCTACGTCTTACCCAACTCCTTCAAATCTGCATTAGTGCCTTTCTTCGCGGGTATTCCTCATCGCACCGGCTGGCGCGGCGAGATGCGCTACGGTTTACTCAACGATGTACGCGTGCTCGATAAAGAAGCCTGGCCGCTAATGGTGGAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGAGTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGG  >  NC_000913/3795095‑3795574

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: