Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,192,392 C→T A82T (GCC→ACC)  thiG ← 1‑deoxy‑D‑xylulose 5‑phosphate:thiol sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,192,3920CT92.1% 161.3 / NA 66A82T (GCC→ACC) thiG1‑deoxy‑D‑xylulose 5‑phosphate:thiol sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (34/24);  minor base G (1/4);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGTTTCCAGTCCCTGATTCGAGCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGCGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACAGATGTGAATCAAACGTTTTGTCCGCAATACGTAACATT  >  NC_000913/4192145‑4192636
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
cgGGTTTCCAGTCCCTGATTCGAGCCAATCGGCGCGCCGAGCGGCCTCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGggtggc                                                                                                                                                                                                                                                     >  1:211752/1‑251 (MQ=255)
          gtccccgataagAGCCAATCGGCGCGCAGAGGGGCCTCACCGATGAACAGCCGACTTCTTCAAGACGTTTAAACCATACCGGATCGGCCCCGCAGTAAGGCAGCACGACCAATCCATGTTGTACCCGCGTTTCGGCGGCTTTACGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGCGAATCTCTACTTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCCGCTGGGTGGCGACAATGGct                                                                                                                                                                                                                                           >  2:167324/12‑251 (MQ=255)
          gccccTAATTCGAGCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAGCCGACTGCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGct                                                                                                                                                                                                                                           <  2:228535/249‑1 (MQ=255)
            cccTGATTCGAGCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGcttc                                                                                                                                                                                                                                         <  1:99677/251‑1 (MQ=255)
            cccTGAGTCGAGCCAAGCGGCGCGCCGACCGGCAGCACCGCTGCACAGCCGACTGCTTCCATACGGTTACACAATACGGGATCGTCCCCGCAGTAAGGCAGCACGACAAATCCCTGGTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCTTGCGTCAGGGTGAAGCGCTAGTTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGTAAATGGcttc                                                                                                                                                                                                                                         <  2:245298/251‑1 (MQ=255)
                gATTCGAGCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTcc                                                                                                                                                                                                                                     <  1:130207/251‑1 (MQ=255)
                gATTCGAGCCAATCATCGGGCCGAGTGGCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCCGTAAGGCAGCACGACAAATCCCTGTTGTAACAGCGTTTCGTCGGCTTTCAGGGTTTCGATGGGATCGGGCACAAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAAAGCT‑TCACGAGCCAGATGGGCGTCGAAAAGGGCTTTTTcc                                                                                                                                                                                                                                     >  2:75243/1‑251 (MQ=255)
                     gAGCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGt                                                                                                                                                                                                                                >  1:81341/1‑251 (MQ=255)
                     gAGCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGt                                                                                                                                                                                                                                >  1:339764/1‑251 (MQ=255)
                         cAATCGGCGCGCCGAGCGGCATCACCGCTGCACAGCCGACTTCTTCCTGACGTTTAAACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑c                                                                                                                                                                                                                           <  2:84909/251‑1 (MQ=255)
                              ggCGCGCCGAGCGGCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATCCCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTACCCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGCAAATGGCTTCTTCCGCTGTTTT‑CGcccc                                                                                                                                                                                                                      >  2:351758/1‑251 (MQ=255)
                              ggCGCCCCGAGCGGCATCACAGCTGCACACCCGCCTTCCTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCCCGACACATCCCTGTTGCACCAGCGTATCGGCGGCTTTCAGGGTTTCGATGGGATCGGCCAACAGCCAGCGGGCGTCCGGGAGACTCTCTAATCTTAACCCGTTTGTGCCTACCGCT‑CCACGAGCCAGATGGGGCGCGCAAAAGGCTTCTTCCGCTGTTTT‑CAcccc                                                                                                                                                                                                                      >  2:193164/1‑251 (MQ=255)
                               gcgcgcCGAGCGGCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCg                                                                                                                                                                                                                     >  2:272748/1‑251 (MQ=255)
                                 gcgcCGAGCGGCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGa                                                                                                                                                                                                                   <  1:406545/251‑1 (MQ=255)
                                     cgggCGGCATCACCGCGGCCCAGCCGACTTCTTCCAGACGGTTACGCAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCGGTTGTACCAGCGTGTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTa                                                                                                                                                                                                              <  2:232438/248‑1 (MQ=255)
                                                     tGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTc                                                                                                                                                                                               >  1:10516/1‑250 (MQ=255)
                                                      gCACAGCCGACTTCTTCCAGACGTGTACACAAGACCGGAGCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTGCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTc                                                                                                                                                                                               <  2:50689/249‑1 (MQ=255)
                                                      gCACAGCCGACTTCTTCCAGACGGTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTcac                                                                                                                                                                                             <  1:15043/251‑1 (MQ=255)
                                                       cacaGCCGACTTCTTCCAGACGTTTACACAACACCGGATCGGCCCCGCAGTAAGTCAGCACGTCAAATCCCTGTTGTACCGGCGTTTCGGCGGCTTTCATGGTTTCGTTGGGATCGGGTAACATCCAGCGGGCGTCAGGGTGAAGCGCTAATGTTATCCAGTTTTTGCCTAACGCT‑GCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCTCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTcaca                                                                                                                                                                                            <  1:56992/251‑1 (MQ=255)
                                                        acaGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCacac                                                                                                                                                                                           <  1:272748/251‑1 (MQ=255)
                                                             cGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACccgcc                                                                                                                                                                                      >  2:311428/1‑251 (MQ=255)
                                                                                acacAATACCGGATCGGCACCGCAGTAAGGCAGAACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCTGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGAGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGACAGATGGGAGGCGAACATGGCTTCTTCCGATGTTTT‑CGCACAGGAT‑GTATTTGGCAGCAGGGGCACACCCGCCGCGATAAGCGGTTCGAgg                                                                                                                                                                    >  2:324694/1‑250 (MQ=255)
                                                                                   caATACCGGATCGGCCCCGCAGTAAGGCCGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGAt                                                                                                                                                                  >  1:132205/1‑249 (MQ=255)
                                                                                    aacccgggccccgCCCCGCATGAAGGCAGCACGACAAATCCTGGTTGGACCAGCGTTCCGGCGTCTTTCAGGGTTTGGAGGGGACCGGCCAACACCCGGCGGGCGGCAGGGTGAATCTCTAATTTTAACCAGTTTGCGCCTACCGCT‑TCACGAGCCAGAGGGGTGGCAAAAATGGCTTCTCCCGCTGTTTT‑CGCCCGGGTT‑GTATGTGGCATAAGGTTCACTCCCGCCGCGATAAGCGGTTCGAGGATAGc                                                                                                                                                               <  2:64554/239‑1 (MQ=255)
                                                                                     aTACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGAATCGGGCAACAGCCAGCGGGCGTCAGGGTGACTCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGAGGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGAAGCAGGGTCACACCCGCCCCGATAAGCGGTTCTGGGATAg                                                                                                                                                                >  2:187239/1‑249 (MQ=255)
                                                                                         cGGAGCGGCCCCGCAGTAGGGCAGCACGACAAATCCCTGTTGTACCGGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCAAAAATGGCTTCTTCCGCTGTTTT‑CGCCCGGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGcgtcgt                                                                                                                                                          <  1:331659/251‑1 (MQ=255)
                                                                                          ggATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGCT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGtt                                                                                                                                                         >  2:329283/1‑251 (MQ=255)
                                                                                               ggccccgaagaagggcagagcGACAATTCCCGTTTGTACCACCGTTTCGGCGGCTTTCAGGTTTTCGGTGGGAGCGGGAAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCATTTTGTGCCTACCGCT‑TCCCGAGCCGGATGGGTGGCGAAAAGGCCTTCTCCCGCTTTTTC‑CGCCGGGGAT‑GTATTTGGCGGCAGGTTCACCCCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCt                                                                                                                                                    <  2:336360/231‑1 (MQ=255)
                                                                                               ggCCCCGAAGTAAGGCAGCACTACAAACCCCTGTTGTCCCGGCGTTTCGGCGGCTTGCAGGGTTCCGTTGGGATCGGGAAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTACCCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTTTTTT‑CGCCCCGGAT‑GTTTTTGGCAGCAGGTTCACACCCGCCGCGATAAGCGGTGCGAGGATAGCGGCGTGGTGCt                                                                                                                                                    <  1:341831/251‑1 (MQ=255)
                                                                                                  cccGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAAAAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGGTCGAGGAAAGCGTCGTTGTGCTGgc                                                                                                                                                 >  2:276515/1‑251 (MQ=255)
                                                                                                       aGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCa                                                                                                                                              >  2:97336/1‑249 (MQ=255)
                                                                                                         tAAGGCAGCAAGACAAATCCATGTTGTACCAGCGTTTCGGCGACTTTCAGGGTTTCGACGGGAACGGGCAACCGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGACACTGGCTTCTTCCGCTGTTTT‑CGCCCCGGCT‑GTCTTCGGCAGCCGGGTCACACCCGCCGCGATAAGCGGTTCGAGGGTAGCGTCGTTGTGCTGGCCCAAGTc                                                                                                                                          >  2:57857/1‑251 (MQ=255)
                                                                                                          aaGGCAGTATGACAATTCCCGGTTGTACCGGCGTTCCGGCGGTTTACAGGGTGTCGATGGGTGCGGCAAACAGCCATCGGGCGGAGGGGTTAATCTTTAAGTTTATCCAGTTTGTGCCTA‑CCCTGCCACGAGCCATATGGGGGGCGAAAAGGGCTTCTCCGGCTTTTTG‑CGCCC‑GGATTGTATGTGGGAGGGGGATCCCACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTc                                                                                                                                          <  1:66041/250‑1 (MQ=255)
                                                                                                           aGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCCAGTCg                                                                                                                                         >  1:312373/1‑250 (MQ=255)
                                                                                                                     aCAAATCCCTGTTGTAACAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCAt                                                                                                                              >  1:316573/1‑251 (MQ=255)
                                                                                                                               gttgtACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCACGGCGCCACGTCTCATCGCCAGTGTc                                                                                                                    >  2:234163/1‑251 (MQ=255)
                                                                                                                                         ggTTTCGGCGGCTTTCAGGGTTTCGATGGGATGGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTGGTGCCTAACGCT‑TCCCGAGCCAGATGGGGGGCGAAAATGGCTTCTTCCGCTGTTTG‑CGCCCGGGAG‑GTATTGGGCAGCGGGGCCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAgctggc                                                                                                          <  1:182239/250‑1 (MQ=255)
                                                                                                                                                  ggcTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGa                                                                                                   >  1:242293/1‑249 (MQ=255)
                                                                                                                                                      ttCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGACTCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGcgcg                                                                                             >  2:152577/1‑251 (MQ=255)
                                                                                                                                                        cAGGGTTTGGAGGGGAGCGGACAAAAGGCAGGGGGGGTCGGGGTGAATCTATAATTTTAACCAGTTTGTCCCTAACGCT‑CCAGGAGCAAGAGGGGTGGCGAAATTGGCTTCTTCCGCTGTTTT‑CGCCCGGGGT‑GTATGGGGCGGCGGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCgg                                                                                            <  2:339764/250‑1 (MQ=255)
                                                                                                                                                          ggtttgcgaggggATCGGGTAACAGCTAGCGGGCTTCGGGGTGAATCGCTAATTTTACCCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATTGGGGGCTAAATTGGCTTCTGCCGCTGTTTT‑CGCCCCGGAT‑GTTGTTGGCAGCGGGGTCTCGCCCGCCGCGATAAGCGGTTCGAGTATAGCTTCGTTGTGCTGGCGCAAGTCGACACGTTGCAGCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATc                                                                                         <  1:34522/241‑1 (MQ=255)
                                                                                                                                                          gggTTTCGATGGGATCGGGCAACAGCCAGCGGGCTTCGGGGTAAATCTCTAATTTTACCCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTCCCGCTGTTTC‑CGCCCGGGAT‑GTATTGGGCAGCAGGGTCACCCCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATc                                                                                         <  1:120779/251‑1 (MQ=255)
                                                                                                                                                              ttCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCt                                                                                     >  1:319374/1‑251 (MQ=255)
                                                                                                                                                               ccgaggggttCGGGCAACAGCCGGGGGGCGTCAGGGTGAATCTATAATTTTAACCAGTTTGTGCCAAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTCCCTCTGTTTT‑CGCCCCGGAT‑GTATTGGGCAGCAGGTTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCt                                                                                     <  1:301496/241‑1 (MQ=255)
                                                                                                                                                                 gATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCCGATTGTGCCTAACGCT‑TCACGATACAGATGGGTGGCGAAAATGGCTTCTTACGATGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTAGAGGATAGCGTCGTTGCGCTGGCGCACGTCGACACGTTTAATCGCCAGTGTCACCAGCTGGCTGCCGGACGCGCGGATCGCCTccc                                                                                  >  2:250649/1‑250 (MQ=255)
                                                                                                                                                                     ggATCGGGCAACAGCCGGCGGGCGTCGGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTACCTCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTTTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCAt                                                                              <  1:280110/251‑1 (MQ=255)
                                                                                                                                                                                        gggCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAAt                                                           >  1:177852/1‑251 (MQ=255)
                                                                                                                                                                                        gggCGTCAGGGCGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAAt                                                           >  1:360739/1‑251 (MQ=255)
                                                                                                                                                                                         ggCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAAtt                                                          >  2:145355/1‑251 (MQ=255)
                                                                                                                                                                                         gacGTCAGGGTGAAGCTCTATTTTAAACCAGGTTGTGCCTAACGCT‑TCCCGCGCCAGAGGGGTGGCGAAAATGGTTGCTGCCGCTGTTTG‑CGCCGGTGGG‑GTATTGGTCAGAAGGGTCACCCCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGTGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAAtt                                                          <  1:265509/249‑1 (MQ=255)
                                                                                                                                                                                              cAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGAGCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCt                                                     >  2:17122/1‑251 (MQ=255)
                                                                                                                                                                                                 ggTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCtgtg                                                  >  1:73975/1‑251 (MQ=255)
                                                                                                                                                                                                    gAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCg                                               >  1:348276/1‑251 (MQ=255)
                                                                                                                                                                                                    gAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAACTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTCCCg                                               >  2:136833/1‑251 (MQ=255)
                                                                                                                                                                                                            aatttaaaccgttgtgtgCTTAACGCT‑TCACGAGCCAGATGGGGGGCGAAATTGGTTTCTTCCGCTGTTTT‑CCCCCGGGGT‑GTATTTGGCAGCGGGGTCACCCCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCTAATTTCCCTGTGCCGGTaaa                                          <  1:162180/234‑1 (MQ=255)
                                                                                                                                                                                                              ttttAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAAc                                         >  1:212943/1‑247 (MQ=255)
                                                                                                                                                                                                                    ccAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACAGAtgtgcat                               >  2:104056/1‑248 (MQ=255)
                                                                                                                                                                                                                     cAGTTTGTGACTAACGCT‑TCACGAGCCAGATTGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCCGGGTCACACCCGCCGCGATACGCGGTTCGAGGATAGCGTCGTTTTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCCGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCCTCAGTTGTGAAGAAGCGAATTTCCCTTTGCCGGTACACCTATGTGAATc                              >  2:40514/1‑251 (MQ=255)
                                                                                                                                                                                                                       gTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTTCCTGTGCCGGTAAACAGATGTGAATCa                             >  1:395613/1‑250 (MQ=255)
                                                                                                                                                                                                                           gtgCCTACCGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CCCCCCGGAT‑GTATTTGGCAGCAGGGTCCCACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACAGATGTGAATCAAACGt                        <  1:351758/251‑1 (MQ=255)
                                                                                                                                                                                                                              ccTAACGCT‑TCACGAGCAAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCGGGAT‑GTTTTTGGCAGCGGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTACCCTGTGCCGGTAAACAGATGTGAATCAAACGtttt                     <  1:112381/251‑1 (MQ=255)
                                                                                                                                                                                                                               cTAACGCT‑TCCCGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTAATCGCCAGTGTCACCAGCTGGCTGCCGGCAGCGCGGATCGCCTCCACCATCCGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACAGAGGGGAATCAAACGTTTTg                    >  2:375571/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aaCGCT‑TCACGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACAGATGTGAATCAAACGTTTTGTc                  >  1:254908/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cT‑TCACGAGCCAGATGGGAGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGCATTTCCCTGTGCCGGTAAACAGATGTGAATCAAACGTTTTGTCCGc               >  2:377339/1‑250 (MQ=255)
                                                                                                                                                                                                                                          ccGAGCCAGAGGGGTGGCAAAAATGGCTTCTCCCGC‑GTTTTGCGCCCCGGTT‑TTATTGGGCAGCGGGGTCCCCCCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACAGATGTGAATCAAACGTTTTGTCCGCAATAc          <  1:272384/250‑1 (MQ=255)
                                                                                                                                                                                                                                          aCGAGCCAGATGGGTGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACAGATGTGAATCAAACGTTTTGTCCGCAATAc          >  1:426982/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  gATGGGTGGCAAAAATGGCTTCTTCCGCTTTTTC‑CGCCCCGGAT‑GTTTTTGGCAGCAGGGTCACCCCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACAGATGTGAATCAAACGTTTTGTCCGCAATACGTAACAtt  <  1:279259/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
CGGGTTTCCAGTCCCTGATTCGAGCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAGGGTGAATCTCTAATTTTAACCAGTTTGTGCCTAACGCT‑TCACGAGCCAGATGGGCGGCGAAAATGGCTTCTTCCGCTGTTTT‑CGCCCCGGAT‑GTATTTGGCAGCAGGGTCACACCCGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACAGATGTGAATCAAACGTTTTGTCCGCAATACGTAACATT  >  NC_000913/4192145‑4192636

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