Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,389,434 T→C 15.4% A107A (GCT→GCC yfbP → TPR‑like repeats‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,389,4340TC15.4% 163.2 / 2.3 75A107A (GCT→GCCyfbPTPR‑like repeats‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (37/24);  new base C (10/2);  total (49/26)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGCATATGGTAACGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCAGCGAATATTGGTTAGCACTGGCAATCCAGTATCGCCGATGTAAAAAAACAGAAGCTGCAGCACAGGCTGCGCTTAATGCTTATCTTGGCAATTGGGCATTTGGTGTTCCTGACAATAAAGTCATTCATTTGCTGTCACAGGCCGCAGA  >  NC_000913/2389210‑2389661
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                   
gcagcaTATGGTAACGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGACGATCAAGACTTTGCAGCAACTCTCTTCCGTGTTGCCCGGa                                                                                                                                                                                                                                 >  1:509957/1‑229 (MQ=255)
gcagcaTATGGTAACGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGGTCGCAGAATGGCTGGAAGTGAATGATGATTATCCCCATGAAAATGGCATCCGCTTCGCAGCTCAAGACTTTTCAGCAAATCCCTTCCGCGTTACCCGGa                                                                                                                                                                                                                                 >  2:328929/1‑229 (MQ=255)
  gggattgttgaagggggagtgttttgtttTTGACTTGCCTTTTGGTCAAGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTGCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGaaa                                                                                                                                                                                                                               <  2:324584/210‑1 (MQ=255)
              cGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGGCTTTCAACg                                                                                                                                                                                                                   >  1:378473/1‑229 (MQ=255)
              cGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGCCTTTCAACg                                                                                                                                                                                                                   >  1:626138/1‑229 (MQ=255)
              cGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGATGCAAATCTCTTCCGTGTTGCTCGGAAAAGACTTTCaagg                                                                                                                                                                                                                   >  2:348514/1‑227 (MQ=255)
                gagaCTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGGAGTTCGAAGAATGGCTGGAAGTGAATGATGATGGTCCCCATGAAAATGGCATCAGCATCGAGGATCAAGACTTTGACGCAAATCTCTTCCCTGTTGGTCGTAAAAGTATTTCacacgg                                                                                                                                                                                                                 >  2:325658/1‑224 (MQ=255)
                 ggacttgtttggTTTGTACTGGCCTTTGGGGCCGGAGCAAAAACTACCTATAGTGTGCGAAACCTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGa                                                                                                                                                                                                                <  2:418180/219‑1 (MQ=255)
                  gaCTGCTGCGGTTTGCACTCGCTTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATTATGATGATCCGCATGAAAATGGCATCAGCAGCGAAGATCAAGACTTTGCAGCAAATCCCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGa                                                                                                                                                                                                                <  1:325658/228‑1 (MQ=255)
                   aCTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGATCGAAGAATGGCTGGAAGGTAATGATGATGATCCCCATGAAAATGTCATCAGCATCGAAGATCAAGACTTTGCAGCAAAGCCTGTCCGTGTTGCCCTGAAATGCTTTTCAAAggagag                                                                                                                                                                                                              >  1:724774/1‑225 (MQ=255)
                            atgtgtatcagagacagggTCAGGAGCTAAAACTTCCTATAGTGTTCGAAACGTACCATGACGAATGGCGAATCGGACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTtgatg                                                                                                                                                                                                      <  1:660172/211‑1 (MQ=255)
                                gtgcgtgctgatgggTCAGGGGCATAAACTACCTTTAGTGTGGTAAACGTACCATGACGAAGGGTGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCa                                                                                                                                                                                                 <  2:303755/216‑1 (MQ=255)
                                      gCCTTTTGTGCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATTACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCa                                                                                                                                                                                           <  1:280019/229‑1 (MQ=255)
                                      gCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGTTCTAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCa                                                                                                                                                                                           >  1:406111/1‑229 (MQ=255)
                                           ttGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGCAGTGAATGATGATGATCCCCATGACAACGGCATCAGCATCGAAGCTCAAGACTTTGCAGCAAATCTCTTCCGTGTTCCACGGAAATGGCTTTCAACGGGAAAGCTGGATCACGCCCTGCCACGGtt                                                                                                                                                                                      >  2:713260/1‑229 (MQ=255)
                                                cAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGCGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCACATCCTTTCCGGGTTGCTCGGAAAGGGCTTTCAACGGGAAGGGTGGATGACCCCCCGCCCCTGTAACAGc                                                                                                                                                                                 >  2:539288/1‑229 (MQ=255)
                                                      cATAAGCTACGTGTTGTGTGCGAAACCTACCATGGCGAATGGCGAATCGTTACTGCCTTTTACAGCATTAAGAAGTTCGAAGAATGGCTTGAAGGTAATTATGATGATCCCCATGAACATGGCATCAGCATCGAAGAGCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATATCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCaa                                                                                                                                                                           <  1:333424/229‑1 (MQ=255)
                                                        tAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGAATTTGCAGCAAATCTCTTCAGTGTTGCCCGGAAATGACATTAAACGGGAAGGCTTACTGACGCACTGccccggtgtcttatcaccatc                                                                                                                                                                         >  2:660171/1‑210 (MQ=255)
                                                           aCTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATTATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCGGTGTTGCCCGGAAATGCCTTTCAACGGGAAGGCTTGTTGACCCCCTGCCACTGTTACCTCGAGCAACAGaa                                                                                                                                                                      >  2:753845/1‑229 (MQ=255)
                                                                    aGTGTGCGAAACGTACCATGACGAATGGCGACTCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGATGGAAGTGAATGATGATGATCACCATGAAAATGGCATCAGCATGGAAGATCAAGACTTTGCAGCAAATTTCTTCCTCGTTGCTCGGAAAAGGCTTTCAACGGGAAGGCTTAATACACCACTGCCACGGTAACAGCTGCCAACAAAAAAAtcacca                                                                                                                                                             >  2:30314/1‑224 (MQ=255)
                                                                      tgtgCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGa                                                                                                                                                           >  1:238725/1‑229 (MQ=255)
                                                                       gtgCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGCATGGCTGGAAGTGAATGATGATGATCCCCATGAACATGGCATCAGCATCGAAGATCAAAACTTTTCAGCAAATCTCTTCAGTGCTGCCCGGCAATGTTTTTCAAAGGGAAGGCTTTATTCCCACCTTCCACTGTTAACCGGAGCAACAGAAAAATTGCCTGaa                                                                                                                                                          >  2:292565/1‑229 (MQ=255)
                                                                       gtgCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGaa                                                                                                                                                          >  1:1620/1‑229 (MQ=255)
                                                                       gtgCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGaa                                                                                                                                                          >  1:606226/1‑229 (MQ=255)
                                                                       gtgCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGaa                                                                                                                                                          >  1:527360/1‑229 (MQ=255)
                                                                       gtgCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCACGGAAATGCATTTAAACGGGAAGGCTTGATGACGCACTGCCCCTGCTACAGCGAGCAACAGAACAATTGcacgaa                                                                                                                                                          >  1:288065/1‑224 (MQ=255)
                                                                       gtgCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCCTTTCCGTGTTGCCCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCCCTGCCACTGTTACAGCGAGCAACAAACCAATTCCCTGaa                                                                                                                                                          >  2:110403/1‑229 (MQ=255)
                                                                       gtgCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCCCTTCCGAGTTGCTCGGAAATGCCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAAAAAAACAAATGCCTGaa                                                                                                                                                          >  2:230476/1‑229 (MQ=255)
                                                                       gtgCGAAACGTACCATGACGAATGGAGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGACGGATGGCTGGAACTGAATGATGATGATCCCCATGAAAATGGCATAAGCATCGAAGATCAAAAATTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGCTGAAGCACTGCCACTGTTACAGCGAGCACCAGAACAATTGCCTGaa                                                                                                                                                          >  1:265875/1‑229 (MQ=255)
                                                                        tgCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGCTGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCCCGGAAATGTCTTTCAACGGGAAGGCTTGTTGACGCCTTGCCCCTGTTACAGCGACCAACGTACAAATTGCCTGagg                                                                                                                                                         >  2:39089/1‑227 (MQ=255)
                                                                        tgCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAATGGATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATAAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGTTGACGCACTGCCACGTTTAAAGCAAGCCACAAAACAAATGCCTGAAg                                                                                                                                                         >  2:311406/1‑229 (MQ=255)
                                                                         tcgaaactgtccATGACGAATGGCGGATCGTACCTTCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGt                                                                                                                                                        <  2:527360/219‑1 (MQ=255)
                                                                         gCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATTACGCACTGCCACTTTTACAGCGAGCAACAGAACAATTCCCTGAAGt                                                                                                                                                        >  1:544179/1‑229 (MQ=255)
                                                                         gCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCAATGCCCCTGTGACACGCACCAACAGTAAAATTTACTGAAgc                                                                                                                                                        >  2:343147/1‑228 (MQ=255)
                                                                          cGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGGGCCACAGAACAATTGACTGAAgca                                                                                                                                                       >  1:363927/1‑227 (MQ=255)
                                                                          cGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTc                                                                                                                                                       >  1:561037/1‑229 (MQ=255)
                                                                           gaaactgtgcATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCa                                                                                                                                                      <  1:713261/220‑1 (MQ=255)
                                                                           gaaactgaacATGCCGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCa                                                                                                                                                      <  2:559446/220‑1 (MQ=255)
                                                                           gAAACGTACCATGACGAATTGCGAATCGTACCTGCCTTTTCCAGCAGTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAAGGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCa                                                                                                                                                      <  2:527503/229‑1 (MQ=255)
                                                                           gAAACGTAACATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCa                                                                                                                                                      <  1:96564/229‑1 (MQ=255)
                                                                             aaCGTACCATGACGAATGGCGAATGGGACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAAGGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCAGc                                                                                                                                                    <  1:46663/229‑1 (MQ=255)
                                                                             aaCGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGGTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCAGc                                                                                                                                                    <  2:1620/229‑1 (MQ=255)
                                                                             aaCGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCa                                                                                                                                                      >  1:559447/1‑227 (MQ=255)
                                                                             aaCGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAAGGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGACAAATTGCCTGAAGTCAGc                                                                                                                                                    >  2:484289/1‑229 (MQ=255)
                                                                             aaCGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCAACAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCCCGGAAATGCCTTTCAACGGGAAGGCTTGATTACGCACTGCCACTGTTACAGGGAGCAACAGAAAAAATGCCTAAAGTccgg                                                                                                                                                    >  2:157997/1‑226 (MQ=255)
                                                                                  aCCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCCCTGTTACACCGAGCACCAAAAAAATTGCCGGAAGTAAggcaaaa                                                                                                                                               >  2:139824/1‑223 (MQ=255)
                                                                                  aCCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTCCAACTGTTCCAGCGAGCACAAAAAAAATTGCCTGAAGCCAGCGaca                                                                                                                                                >  2:134299/1‑226 (MQ=255)
                                                                                       gACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGGGAATGATGATGATCACCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCCCTTCCGTGTTGCTCGGAAATGTCTTTACACGGGAAGGCTTGATGACGCACGGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTAAGCgcgaattgcg                                                                                                                                          >  1:62895/1‑220 (MQ=255)
                                                                                            aTGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGCCAGCGAATATTGGTTTGCa                                                                                                                                     >  1:594840/1‑229 (MQ=255)
                                                                                                         ccTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGTTGGAGGTGAATGATGCTAATCCCAACAAAAATGGCATCAGCATCGAAGAGCAAGACTTTGAAGCAAATCTCTACCGGTTTGCTCGGAAATGTCTTTCAACGTGAAGCCTTGTTTACGACCTGCCACTGTTACATCGAGTCACAGAACAACTGCCTGAAGTCCGCGAAAATCGTTTacaccagccaaacaata                                                                                                                        >  2:327180/1‑213 (MQ=255)
                                                                                                                         ttAAGAAGTTAGAAGAAGGGCTGGAATTGAATGATGATGATCCCCATGAAAATGGCATCAGCATGGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCAGCGAATATTGGTTAGCACTGGCAATCCAGTATCGCCGATGTaaaaa                                                                                                        <  1:328929/229‑1 (MQ=255)
                                                                                                                                 ttCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCCCGGAAATGTCTTTCAACGGGAAGGCTTGCTGACGCACTGCCCCTGTTACAGCGAGCAACAGAACAATTGCCTGCAGTCAGCGAAAATTGGTTAGCACTGCCAATCCAGTAACCCCGCAGTAAAAAAAAAGaaa                                                                                                >  2:195255/1‑228 (MQ=255)
                                                                                                                                   cGAAGAATGGCTGGAAGTGAATGATGATGATCCTCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCAGCGAATATTGGTTAGCACTGGCAATCCAGTATCGCCGATGTAAAAAAACAGAAGCt                                                                                              <  1:157997/229‑1 (MQ=255)
                                                                                                                                         attggttgAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCAGCGAATATTGGTTAGCACTGGCAATCCAGTATCGCCGATGTAAAAAAACAGAAGCTgcagca                                                                                        <  2:378473/222‑1 (MQ=255)
                                                                                                                                         agggCTGGAAGTGAATGAGGATGCGCCACAGTAAAAGGGCCTCAGGATGGAAGATCAAGACTTTGCAGCAAAGCTCTTCCGGGTTGCTCGGAAATGTCTTTCAACGGGAAGGTGTGATGAAGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCGGAAGTCAGCGAATATTGGTTAGCACTGGCAATCCAGTATCGCCGATGTAAAAAAACAGAAGCTgcagca                                                                                        <  2:129415/227‑1 (MQ=255)
                                                                                                                                              tGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGCCTTTGCAGCAAATCACTACCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCAGCGAATATTGGTTAGCACTGGCAATCCAGTATCGCCGATGTAAAAAAACAGAAGCTGCAGCACAgg                                                                                    <  1:715439/228‑1 (MQ=255)
                                                                                                                                                  aGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCACATCTCTTCCGTGTTGCTCGGAAATGTCTTTTAAAGGGACGGCTTGATGAGCCACTGCCACTGTTACAGCGATCAACATACAAATTGCCTGAAGTCAGCGAATATAGATTAGCACTGTCAATCACGTATCCCCCATGTAAAAAAACAGaccgtgaaacacagacgacc                                                                               >  2:229828/1‑210 (MQ=255)
                                                                                                                                                  aGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCAGCGAATATTGGTTAGCCCTGGCAATCCAGTATCGCCGATGTAAAAAAAAAGAAGGTGCGGCAAAgggcggg                                                                               >  2:263572/1‑224 (MQ=255)
                                                                                                                                                  aGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCAGCGAATATTGGTTAGCACTGGCAATCCAGTATCACCGATGTAAAAAAACAGAAGATGCAGCACAGGCTgcg                                                                               >  1:607850/1‑229 (MQ=255)
                                                                                                                                                    tGAATGATGATGACCCCCAGGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCAGCGAATATTGGTTAGCACTGGCAATCCAGTATCGCCGATGTAAAAAAACAGAAGCTGCAGCACAGGCTGCGCt                                                                             <  1:39089/229‑1 (MQ=255)
                                                                                                                                                    tGAATGATGAAGATCCCCATGAAAATGGCATCAGCAGCGAAGATCAAGGCTTTGCAGCAAATCTCTTCCGTGTTGCCCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCAGCGAATATTGGTTAGCACTGGCAATCCAGTATCGCCGATGTAAAAAAACAGAAGCTGCAGCACAGGCTGCGCt                                                                             <  1:134299/229‑1 (MQ=255)
                                                                                                                                                      aaTGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGATTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCAGCGCATCTTGGTTAGCACTGCCAATCCAGAATGGCCGAGGAAAGAAACAAGAAGCTGCAGCACAGGCTGGGCTTa                                                                           >  1:349128/1‑229 (MQ=255)
                                                                                                                                                                     cATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACACTTGCCTGAAGTCAGCGAATATTGGTTAGCACTGGCAATCCAGTATCGCCGATGTAAAAACACAGAAGCTGCAGCAAAGGCTGCCCTTAATGCTTATCTTGACa                                                            >  2:441198/1‑229 (MQ=255)
                                                                                                                                                                             tGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCAGCGAATATTGGTTAGCACTGGCAATCCAGTATCGCCTATGTAAAAAAACCGAAGCTGCAGCACAGGCTGCGCTTAATCCTTATTTTGGCAATTGGGc                                                     >  1:363385/1‑228 (MQ=255)
                                                                                                                                                                                cATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCAGCGAATATTGGTTAGCACTGGCAATCCAGTATCGCCGATGTAAAAAAACAGAAGCTGCAGCACAGGCTGAGCTTAATGCTTATCTTGGCAATTGGGCAttt                                                 >  1:557856/1‑229 (MQ=255)
                                                                                                                                                                                 aTCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCCCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCAGCGAATATTGGTTAGCACTGGCAATCCAGTATCGCCGATGTAAAAAAACAGAAGCTGCAGCACAGGCTGCGCTTAATGCTTATCTTGGCAATTGGGCATTTg                                                <  1:186293/229‑1 (MQ=255)
                                                                                                                                                                                          gAAGATCAAGACTTTGCAGAAAATCTCTTCCGTGGTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCAGCGAATATTGGTTAGCACTGGCAATCCAGTATCGCCGATGTAAAAAAACAGAAGCTGCAGCACAGGCTGCGCTTAATGCTTATCTTGGCAATTGGGCATTTGGTGTTCCt                                        <  1:292565/228‑1 (MQ=255)
                                                                                                                                                                                                 aaGTCTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCAGCGAATATTGGTTAGCACTGGCAATCCAGTATCGCCGATGTAAAAAAACAGAAGCTGCAGCACAGGCTGCGCTTAATGCTTATCTTGGCAATTGGGCATTTGGTGTTCCTGACAATa                                 <  1:230476/228‑1 (MQ=255)
                                                                                                                                                                                                 aaGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCAGCGAATATTGGTTAGCACTGGCAATCCAGTATCGCCGATGTAAAAAAACAGAAGCTGCAGCACAGGCTGCGCTTAATGCTTATCTTGGCAATTGGGCATTTGGTGTTCCTGACAATa                                 <  1:205699/228‑1 (MQ=255)
                                                                                                                                                                                                        tGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCAGCGAATATTgg                                                                                                                                           <  1:24262/115‑1 (MQ=255)
                                                                                                                                                                                                        tGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCAGCGAATATTgg                                                                                                                                           >  2:24262/1‑115 (MQ=255)
                                                                                                                                                                                                        tGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCAGCGAATATTGGTTAGCACTGGCAATCCAGTATCGCCGATGTAAAAAAACAGAAGCTGCAGCACAGGCTGCGCTTAATGCTTATCTTGGCACTTGGGCATTTGGTGTTCCTGACAATAAAg                              >  1:186563/1‑224 (MQ=255)
                                                                                                                                                                                                             caAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGACTGAAGTCAGCGAATATTGGTTAGCACTGGCAATCCAGTATCGCCGAGGTAAAAAAACAGAAGCTGCAGCACAGGCTGGGCTTAATGCTTATCTTGGCAATTGGGCATTTGGGGTTCCTGACAATAAAGTCATTCAttt                    >  1:466950/1‑229 (MQ=255)
                                                                                                                                                                                                                               cTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCAGCGAATATTGGTTAGCACTGGCAATCCAGTATCGCCGATGTAAAAAAACAGAAGCTGCAGCACAGGCTGCGCTTAATGCTTATCTTGGCAATTGGGCATTTGGTGTTCCTGACAATAAAGTCATTCATTTGCTGTCACAGACCGCAGa  >  1:191930/1‑229 (MQ=255)
                                                                                                                                                                                                                               cTCGGAAATGTCCTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCAGCGAATATTGGTTACCACTGGCAATCCAGTATCGCCCATGTAAAAAAACAGAAGCTGCAGCACAGGCTCCGCTTAATGCTTATTTTGGCAATTGGGCATTTGGTGTTCCTAACAATAAAATCACTCATTTGCCGTCACAGGCCCCAg   >  2:4980/1‑228 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                   
GCAGCATATGGTAACGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACGCACTGCCACTGTTACAGCGAGCAACAGAACAATTGCCTGAAGTCAGCGAATATTGGTTAGCACTGGCAATCCAGTATCGCCGATGTAAAAAAACAGAAGCTGCAGCACAGGCTGCGCTTAATGCTTATCTTGGCAATTGGGCATTTGGTGTTCCTGACAATAAAGTCATTCATTTGCTGTCACAGGCCGCAGA  >  NC_000913/2389210‑2389661

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