Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,571,975:1 +G 7.1% coding (71/1401 nt) gadB ← glutamate decarboxylase B, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,571,9731.G 7.1% 223.1 / 12.0 56coding (73/1401 nt)gadBglutamate decarboxylase B, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base . (23/29);  new base G (2/2);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATT  >  NC_000913/1571689‑1572256
                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                           
aaCGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTggta                                                                                                                                                                                                                                                                                         >  1:221446/1‑287 (MQ=255)
aaCGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGaa                                                                                                                                                                                                                                                                                         >  1:571840/1‑289 (MQ=255)
              ttttCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCAcc                                                                                                                                                                                                                                                                           >  2:417750/1‑289 (MQ=255)
              ttttCGGCCCAGGCGCATGCAACAGATCGGCAACCAGATTTACGCAACGCAGGTCGATGGCTGCGGATTGCAGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCAcc                                                                                                                                                                                                                                                                           <  1:617048/289‑1 (MQ=255)
                           cgcaTGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGaa                                                                                                                                                                                                                                                                 <  1:466128/286‑1 (MQ=255)
                                                      tACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACt                                                                                                                                                                                                                                   <  2:74314/289‑1 (MQ=255)
                                                             cgcaGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGcttct                                                                                                                                                                                                                            >  2:51325/1‑289 (MQ=255)
                                                                  gTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTAt                                                                                                                                                                                                                         >  2:10750/1‑288 (MQ=255)
                                                                  gTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATcc                                                                                                                                                                                                                       >  1:605369/1‑289 (MQ=255)
                                                                   tCGTTGGCTGCGGATGTCGTATATTCTTCTTGGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTAt                                                                                                                                                                                                                         <  2:705231/286‑1 (MQ=255)
                                                                     gATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCAt                                                                                                                                                                                                                     >  1:203807/1‑289 (MQ=255)
                                                                           tgcggattgcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTaaa                                                                                                                                                                                                               <  2:605369/288‑1 (MQ=255)
                                                                            gcggattgcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAAATTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACt                                                                                                                                                                                                                                   >  2:209736/1‑267 (MQ=255)
                                                                            gcggattgcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAAATTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACt                                                                                                                                                                                                                                   <  1:209736/267‑1 (MQ=255)
                                                                                 ttgcggatATTCTTCTTTGTCGATCCAGTTTTTTTTAATGGATAAATCCATCAATTTGTGGCCATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTa                                                                                                                                                                                                        <  2:418213/289‑1 (MQ=255)
                                                                                  tggggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGGTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTaa                                                                                                                                                                                                       <  2:620719/286‑1 (MQ=255)
                                                                                     ggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTaaaa                                                                                                                                                                                                     <  2:508762/288‑1 (MQ=255)
                                                                                              tcttTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATc                                                                                                                                                                                           >  1:346678/1‑289 (MQ=255)
                                                                                                  tGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCAt                                                                                                                                                                                       >  1:25980/1‑289 (MQ=255)
                                                                                                     cGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTaa                                                                                                                                                                                    <  2:25980/289‑1 (MQ=255)
                                                                                                     cGATCCAGTTTTGGTTAATGGATAAATCCATCAATTTGTGGACATTGTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTaa                                                                                                                                                                                    <  2:141153/289‑1 (MQ=255)
                                                                                                       aTCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTaaaa                                                                                                                                                                                  >  2:641454/1‑289 (MQ=255)
                                                                                                                ttGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGAtt                                                                                                                                                                        <  1:275794/289‑1 (MQ=255)
                                                                                                                       tGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGa                                                                                                                                                                  >  2:450319/1‑289 (MQ=255)
                                                                                                                           tAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATaaa                                                                                                                                                              <  1:417750/289‑1 (MQ=255)
                                                                                                                                 cATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGc                                                                                                                                                        <  1:51325/289‑1 (MQ=255)
                                                                                                                                   tCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCaa                                                                                                                                                      <  1:450319/289‑1 (MQ=255)
                                                                                                                                   tCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCaa                                                                                                                                                      >  1:399390/1‑289 (MQ=255)
                                                                                                                                        ttGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCa                                                                                                                                                 <  1:80992/289‑1 (MQ=255)
                                                                                                                                               cATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACt                                                                                                                                          <  2:346678/289‑1 (MQ=255)
                                                                                                                                               cATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACt                                                                                                                                          <  2:28974/289‑1 (MQ=255)
                                                                                                                                                              aGGACTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACa                                                                                                                           <  2:212206/289‑1 (MQ=255)
                                                                                                                                                                 tCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACa                                                                                                                           <  2:203807/287‑1 (MQ=255)
                                                                                                                                                                     gCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCt                                                                                                                    >  2:316285/1‑289 (MQ=255)
                                                                                                                                                                               gCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGttattt                                                                                                          >  2:96448/1‑289 (MQ=255)
                                                                                                                                                                                      tCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTg                                                                                                    <  1:10750/289‑1 (MQ=255)
                                                                                                                                                                                              gCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTa                                                                                            >  2:682782/1‑288 (MQ=255)
                                                                                                                                                                                               cGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACa                                                                                                                           <  1:84346/256‑1 (MQ=255)
                                                                                                                                                                                               cGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACa                                                                                                                           >  2:84346/1‑256 (MQ=255)
                                                                                                                                                                                               cGTTGCCATCAAGATATAAGTGGTCATTGAGAATCTGGAATGCGACAGCGTCGCGCATTTCGTGCAGCGGGAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTATTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATAGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTAc                                                                                           <  2:251675/288‑1 (MQ=255)
                                                                                                                                                                                                   gCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGaa                                                                                      <  1:316285/289‑1 (MQ=255)
                                                                                                                                                                                                         aaGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATAGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTaa                                                                                                                                 <  2:101848/240‑1 (MQ=255)
                                                                                                                                                                                                         aaGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATAGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTaa                                                                                                                                 >  1:101848/1‑240 (MQ=255)
                                                                                                                                                                                                                aaTTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCa                                                                         >  2:28305/1‑289 (MQ=255)
                                                                                                                                                                                                                      tCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGaaa                                                                   >  1:197976/1‑289 (MQ=255)
                                                                                                                                                                                                                       cATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGaaaa                                                                  >  1:369979/1‑289 (MQ=255)
                                                                                                                                                                                                                            aTAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTaa                                                                                                                                 >  2:372051/1‑221 (MQ=255)
                                                                                                                                                                                                                            aTAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTaa                                                                                                                                 <  1:372051/221‑1 (MQ=255)
                                                                                                                                                                                                                                    gAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGaa                                                                                                                                                                                                                                                                 <  2:408534/85‑1 (MQ=255)
                                                                                                                                                                                                                                    gAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGaa                                                                                                                                                                                                                                                                 >  1:408534/1‑85 (MQ=255)
                                                                                                                                                                                                                                      aTGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGaa                                                   >  1:299358/1‑289 (MQ=255)
                                                                                                                                                                                                                                            cATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACa                                                                                                                           <  2:104787/211‑1 (MQ=255)
                                                                                                                                                                                                                                            cATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACa                                                                                                                           >  1:104787/1‑211 (MQ=255)
                                                                                                                                                                                                                                                       aTTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTa                                  <  1:641454/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      cGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACt                                                                                                                                          <  2:577945/154‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      cGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACt                                                                                                                                          >  1:577945/1‑154 (MQ=255)
                                                                                                                                                                                                                                                                                       gATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCAtt  <  1:28305/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                           
AACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATT  >  NC_000913/1571689‑1572256

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: