Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,446,482 C→T L791L (CTG→TTG)  tamB → translocation and assembly module subunit TamB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,446,4820CT98.3% 182.2 / NA 59L791L (CTG→TTG) tamBtranslocation and assembly module subunit TamB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (25/32);  minor base G (0/1);  total (26/33)
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.90e-01

CCCGAATGCGGAACTGTGCGTGCCGCAAACTATCGATGCGGGGGCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGCTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCGTGGGGAAAAAGCAGCGGGGATGGTGAGTG  >  NC_000913/4446235‑4446713
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        
cccGAATGCGGAACTGTGCGTGCCGCAAACTATCGATGCGGGGGCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGc                                                                                                                                                                                                                                       >  1:125754/1‑251 (MQ=255)
       gCGGAACTGTGCGTGCCGCAAACTATCGATGCGGGGGCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGAAACCACGAAAGAAGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGCCGTAAAGGGCAGGTAACGCAAAACGTCAACGATGCGGCGTTGACCGTGGc                                                                                                                                                                                                                                >  2:340049/1‑251 (MQ=255)
        cGGAACTGTGCGTGTCTCAAACTATCGATGTGGTGGCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGTCCAGCTGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACTTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCg                                                                                                                                                                                                                               <  1:415305/251‑1 (MQ=255)
               gtgcgtgcCGCAAACTATCGATGCGGGGGCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGGGGCGTTTCAg                                                                                                                                                                                                                         <  1:135007/250‑1 (MQ=255)
                        gCAAACTATCGATGCGGGGGCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATc                                                                                                                                                                                                               >  1:378314/1‑251 (MQ=255)
                          aaaCTATCGATGCGGGGGCCGAAGGGCGTGCGGTAGTGGATCTCAGCCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGGTGCGGCGTTGCCCGTGGCGTTTCAGCCACTGAATCTg                                                                                                                                                                                                             >  1:163501/1‑251 (MQ=255)
                               aTCGATGCGGGGGCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGCCGTTGCCTGGGACACCACGAAAGAGGGTCTGCCGCAGGGCGGTATCACCCTTTCGGGGCGGAAAGTGAAGGGAACGCAAACCGTCAACGATGCGGCGAGGCCCGTGGCGTTTCACACACTGGATCTggaggg                                                                                                                                                                                                        >  2:48142/1‑246 (MQ=255)
                                   aTGCGGGGGCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGa                                                                                                                                                                                                     >  1:258154/1‑250 (MQ=255)
                                      ggggggCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTg                                                                                                                                                                                                 <  1:192083/250‑1 (MQ=255)
                                          ggCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTa                                                                                                                                                                                             <  1:459689/251‑1 (MQ=255)
                                                     gggagatggtgAACCGCAACCGCTTCGTATGCGCCATGCTGAAACCGTTGATGCCAGAAACCACGCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCGGCCGCAGGGCAGTATCTCCCGTTCGGTGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTTCCCGTGGCGTTTCAGACACTGAATCTGACGGCTGAATTGCGTTACAACCTTGcc                                                                                                                                                                                  <  2:34580/245‑1 (MQ=255)
                                                         ggtggtGAATCTCAACAGCTTCGACCTAGCCATGATGAAACCGTTTATGCCAGAAACCACTCAGGCCCGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGAGACACCACGCAAAAGAGACAGCCGCAGAGCAGTATCACACTTTCGGGGCGTAACGTGCAAGTAAAGACAACCGTCAACGATGCGGCGTAGACCGGAGCATTGACGACACTGAGTCTGACGGCGGAATTGCATAACAAACGTTCCGAAt                                                                                                                                                                              >  2:297950/1‑251 (MQ=255)
                                                             gtgAATCTCAACCGTTTCGTCACCGCCAGGCTGAAACCGTGTAGTCCAGAAACCACTCAGGCCAGCGGTATCTTCTCGGGTAAAGCGGACGTTGCCTGGGACCCCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCGTGCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGAGGCGTTTCAGACACTTATTCTGACGGCGGAATTGCGTAACAACCGTTCCGAATTggg                                                                                                                                                                          <  2:450762/251‑1 (MQ=255)
                                                               gAATCTCAACCGCTTCGCCCTCGCCATGCTGAACCCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCt                                                                                                                                                                        <  2:339489/251‑1 (MQ=255)
                                                                          gCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGc                                                                                                                                                             >  2:175405/1‑251 (MQ=255)
                                                                          gCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGGCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAAACGTGCAGACTAGGGCTGGGCCAACCGc                                                                                                                                                             >  2:264620/1‑251 (MQ=255)
                                                                            ttCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCt                                                                                                                                                           >  1:368057/1‑251 (MQ=255)
                                                                              cGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGCATTGGGCTGGACCATCCGACTgc                                                                                                                                                         >  1:469278/1‑250 (MQ=255)
                                                                                 ccTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACTATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCTTGCCGAATTGGGCTGGACCATCCGCCTGACCa                                                                                                                                                      <  1:404795/251‑1 (MQ=255)
                                                                                 ccTCGCCATGCTGAAACCGGTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTACCTGGGACCCCACGAACGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGCGCAGGTAACGCAAACCGTCAACGATACGCCGTTGCCCGTGGCGTTTCAGCCACTGAAAATGCCGGCGGCATTGCGTAACAACCGAGCCGACTTGGGCTGGACCATCCGCCCGAca                                                                                                                                                       >  2:275221/1‑249 (MQ=255)
                                                                                   ccGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAAt                                                                                                                                                    <  1:175405/250‑1 (MQ=255)
                                                                                    cGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATa                                                                                                                                                   <  1:339185/251‑1 (MQ=255)
                                                                                      ccATGCTGAAACCGTTTATGCCAGAAACCACGCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAAc                                                                                                                                                 <  2:415253/251‑1 (MQ=255)
                                                                                       cATGCTGAAACCGTTTATGCCCGAAACCAATCAGACCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCCGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTACCGAATTGGGCTGGCCCCTCCGCCTGACCAATAACg                                                                                                                                                >  2:66195/1‑251 (MQ=255)
                                                                                       cATGCTGAAACCGTTTAGGCCAGAAACCCCTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACg                                                                                                                                                <  2:400468/251‑1 (MQ=255)
                                                                                             gAAACCGTTTATGCCAGAATCCACTCAGGCCGGCGTTATCTTCACGGGTAAAGCTGACGTTGCCTTGGACACCACGAAAGAGGGGCTGCTGCAGGGCAGTATCCCCCTTTCGGGGCGTAACGTGCAGGTATCGCAAACCGTCAACGAGGCGGCGGTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAAGTTCGTAACAACCGTGCCGAATTGGGCTGTACCATCCGCCTGACCAATAACGGccc                                                                                                                                            <  2:452017/249‑2 (MQ=255)
                                                                                             gAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCCTCCGCCTGCCCAATa                                                                                                                                                   >  2:53113/1‑242 (MQ=255)
                                                                                             gAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATa                                                                                                                                                   <  1:53113/242‑1 (MQ=255)
                                                                                                    tttATGCCAGAAACCCCTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGAt                                                                                                                                     <  2:475426/249‑1 (MQ=255)
                                                                                                      tATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGa                                                                                                                                  <  1:371812/250‑1 (MQ=255)
                                                                                                                       aGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTTCAGGGGACCGAt                                                                                                                >  1:420191/1‑251 (MQ=255)
                                                                                                                         gCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATcc                                                                                                              <  2:378314/251‑1 (MQ=255)
                                                                                                                                gTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGCCAGGTGCAGGTGACCGATCCGCCAGgc                                                                                                       >  1:3031/1‑251 (MQ=255)
                                                                                                                                      tCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGgccgcc                                                                                                   <  1:265310/249‑1 (MQ=255)
                                                                                                                                      tCACGGGGAAAGCGGCCCTTGCCTGGGACCCCACGTAAGAGGGGAGGCCGCAGGGCGGTATCACCCGTTCGGGGCGTAACTTGCAGGTAACGCAAACCGTCAACGAGGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGccgccg                                                                                                  <  2:163501/250‑1 (MQ=255)
                                                                                                                                              aaaGCGGATGTTGCCTGGGACTCCCCGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGTGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGTCTTTTCAGACACTGAATCGGACGGCGGAATTGCGTAACAACCGTTCCGAATTGGGCTGGACCATCCGCCTGGCCAATAACGTCCAGTTTGATGTACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTtgg                                                                                         <  1:458228/251‑1 (MQ=255)
                                                                                                                                               aaGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGAAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGAACCGCAAGACCGACGTAATCTtggg                                                                                        >  1:322532/1‑250 (MQ=255)
                                                                                                                                                    gACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGa                                                                                                                                                                     <  1:19318/169‑1 (MQ=255)
                                                                                                                                                    gACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGAGCAGGCAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGa                                                                                                                                                                     >  2:19318/1‑169 (MQ=255)
                                                                                                                                                       gTTGCCTGGGACACCACGAACGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGGGCAGGTGACCGTTCCGCAAGGCCGCCGTAATCTTGGTGGCACCGt                                                                                >  2:472200/1‑251 (MQ=255)
                                                                                                                                                          gCCTGGGACCCCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTGCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGTTGGTCCATCCGCCTGACCAATAACGGCCAGTTGGAGGGACAGGGGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCa                                                                              <  1:316007/250‑1 (MQ=255)
                                                                                                                                                              gggACACCACGAAAGAGGGGCTGGCGCAGGGCAGTATCACCCGTTCGGGGCGTAACGTGCAGGTAACGCAACCCGTCAACGATGGGGCTTTGCCCGTGGCGTTTCAGACACTGAATCTGTCGGCGGAATTGCGTAACAACCGGTCCGAATTGGGCTGGTCCATCCGCGTGACCAATAAGGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAatat                                                                         <  1:313026/250‑1 (MQ=255)
                                                                                                                                                                         aaaGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTACAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATACGTAACTTCa                                                              >  1:339941/1‑250 (MQ=255)
                                                                                                                                                                          agggggggCTGCCGCGGGGCAGGTCCCCCCTGTGGGGGCGTAACGTGCAGGTAACGAAACCCGTCAACGGTGGGGGGTTGCCCGGGGCGTTTCAGACACGGAATCGGACGGCGGAATGCCGTAACAACCGTGCCGATTTGGTCGGGACAACCCGCCTGCCCAATAACGGCCAGTTTGATGGAAAGGTGCAGGTTACCGATCCTCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAAc                                                            <  2:222089/247‑1 (MQ=255)
                                                                                                                                                                           agagGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCACCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAAcc                                                           <  1:466506/251‑1 (MQ=255)
                                                                                                                                                                             gggggCTGCCGCAGGGCAGTATCCCCCTTGCGGGGCGTAACGTGCAGGTAACTCAACCCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAAGAACGGCCAGTTTTTTGGACAGGTGCAGGTTCCCGATCTGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCtt                                                         <  1:230049/250‑1 (MQ=255)
                                                                                                                                                                                ggCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCg                                                      <  1:43170/251‑1 (MQ=255)
                                                                                                                                                                                ggCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCAAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCg                                                      <  1:495240/251‑1 (MQ=255)
                                                                                                                                                                                        cAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGCATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGgccgcc                                                                                                   >  2:408397/1‑199 (MQ=255)
                                                                                                                                                                                        cAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGgccgcc                                                                                                   <  1:408397/199‑1 (MQ=255)
                                                                                                                                                                                         aGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAGCGTCAATATCCGTAACTTCCACCTTGCGATGCTAAAc                                             >  2:116657/1‑251 (MQ=255)
                                                                                                                                                                                             cAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAAcccc                                          >  1:301703/1‑250 (MQ=255)
                                                                                                                                                                                                 aTCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCtt                                     <  1:367030/251‑1 (MQ=255)
                                                                                                                                                                                                   cacccttttggtggtaaaggTGCAGGTACCGTAAACCGTAAACGATGCGGCGTTGCCGGGGGCTTTCCAGACACTGAAGCTGACGTCGGAATTGCGTAACAACCGTGCCGATTTGGGCTGGACCATCCGCGCGACCAATAACGGATAGTTTGATGGACAGGTGCAGGTGCCCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTGCAACCTTGCGATGATAAACCCCATCTTTa                                   <  1:243663/232‑1 (MQ=255)
                                                                                                                                                                                                         gtgggggCTTAACGGGCAGGTACCGCAAACCGCCAACGATGCGGGGTTGCCCGTGGCGTTTCAGACATGGAATCGGACGGAGGAATTGGTTAACAACCGTGCGGAATTGGGCTTGCCCACCCGCCTGCCAAATAACGGCCAGTTGGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCg                               <  2:372227/246‑1 (MQ=255)
                                                                                                                                                                                                                      gTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGAGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCCTCCGCCTGACCAATCACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCACATACCGTAACTTCACCCTTGCGATGATAAACCCCATATTTACCCGTGGGGacaaaactgc                >  2:429227/1‑246 (MQ=255)
                                                                                                                                                                                                                          aGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTCACAACCGTGCCGAATTGGGCTGGCCCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGAACGATCCGCAAGGCCGCCGTAATCTTGGTGGCCACGTCACTATCCGTACATTCAACCTTGCGAAGATAAACCCCATCTTTACCCAGGGGGAAAAAGCAGAgggg            >  2:160772/1‑251 (MQ=255)
                                                                                                                                                                                                                           ggTAACGCAAACCGTAAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCGTGGGGAAAAAGTAGCGGGGa           <  1:441324/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aaCCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCGTGGGGAAAAAGCAGCGGGGATGGTGAGTg  >  1:231909/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        
CCCGAATGCGGAACTGTGCGTGCCGCAAACTATCGATGCGGGGGCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGCTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCGTGGGGAAAAAGCAGCGGGGATGGTGAGTG  >  NC_000913/4446235‑4446713

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