Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 438,740 T→G 5.0% N480H (AAC→CAC)  dxs ← 1‑deoxyxylulose‑5‑phosphate synthase, thiamine triphosphate‑binding, FAD‑requiring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913438,7400TG 5.0% 669.6 / 20.1 194N480H (AAC→CAC) dxs1‑deoxyxylulose‑5‑phosphate synthase, thiamine triphosphate‑binding, FAD‑requiring
Reads supporting (aligned to +/- strand):  ref base T (89/95);  new base G (5/5);  total (94/100)
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.65e-01

Only 100 of 199 total aligned reads displayed.

AGCGTGGCGTTCAGCGATTCGGCGACTTTCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAAGAGAGATCAAAAGCACCCTGATGG  >  NC_000913/438592‑438888
                                                                                                                                                    |                                                                                                                                                    
aGCGTGGCGTTCAGCGATTCGGCGACTTTCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGtt                                                                                                                                                      <  4:1534984/149‑1 (MQ=255)
aGCGTGGCGTTCAGCGATTCGGCGACTTTCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGtt                                                                                                                                                      <  1:740796/149‑1 (MQ=255)
 gCGTGGCGTTCAGCGATTCGGCGACTTTCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTg                                                                                                                                                     <  1:1722532/149‑1 (MQ=255)
        gTTCAGCGATTCGGCGACTTTCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCAcgcg                                                                                                                                              >  4:1703660/1‑149 (MQ=255)
              cGATTCGGCGACTTTCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAgc                                                                                                                                        <  4:1318102/149‑1 (MQ=255)
               gATTCGGCGACTTTCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAgcg                                                                                                                                       >  2:654219/1‑149 (MQ=255)
                aTTCGGCGACTTTCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAgcgc                                                                                                                                      >  2:437485/1‑149 (MQ=255)
                          tttCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACg                                                                                                                            >  3:1150453/1‑149 (MQ=255)
                           ttCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACgg                                                                                                                           <  1:1742291/149‑1 (MQ=255)
                           ttCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACgg                                                                                                                           <  3:86351/149‑1 (MQ=255)
                                 gcTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATc                                                                                                                     <  2:1722076/149‑1 (MQ=255)
                                 gcTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATc                                                                                                                     <  4:810501/149‑1 (MQ=255)
                                 gcTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATc                                                                                                                     <  1:905413/149‑1 (MQ=255)
                                 gcTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATc                                                                                                                     <  1:365090/149‑1 (MQ=255)
                                     tgcataccggaaatggtcattatgaccccgagcgatgaaaacgaatgtcgccagatgctctataccggctatcactataacgatggccCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTta                                                                                                                 >  4:296763‑M1/89‑149 (MQ=255)
                                       ggCATCAGCGTACCACAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGTTAGTTTTTCCAGCGGCGTCAGTTCCACGCTGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTtata                                                                                                               <  4:1347104/149‑1 (MQ=255)
                                       ggCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGACTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTtata                                                                                                               <  1:14164/149‑1 (MQ=255)
                                          aTCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTg                                                                                                            >  3:328221/1‑149 (MQ=255)
                                          aTCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTg                                                                                                            >  2:1632286/1‑149 (MQ=255)
                                             aGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATa                                                                                                         <  2:1555333/149‑1 (MQ=255)
                                                  aCCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCgg                                                                                                    <  2:609933/149‑1 (MQ=255)
                                                   ccAACGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTa                                                                                                                                          >  4:570134/1‑110 (MQ=255)
                                                   ccAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTa                                                                                                                                          <  3:570134/110‑1 (MQ=255)
                                                   ccAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGt                                                                                                   <  1:208236/149‑1 (MQ=255)
                                                   ccAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGt                                                                                                   >  3:1117821/1‑149 (MQ=255)
                                                    cAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGta                                                                                                  >  2:21565/1‑149 (MQ=255)
                                                    cAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGta                                                                                                  >  4:1149221/1‑149 (MQ=255)
                                                     aaaGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGtat                                                                                                 >  2:1741415/1‑149 (MQ=255)
                                                      aaGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGtata                                                                                                <  1:160549/149‑1 (MQ=255)
                                                      aaGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGtata                                                                                                >  1:1087785/1‑149 (MQ=255)
                                                        gTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATaga                                                                                              <  1:1444837/149‑1 (MQ=255)
                                                        gTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATaga                                                                                              <  2:1087785/149‑1 (MQ=255)
                                                         ttAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATagag                                                                                             >  4:142928/1‑149 (MQ=255)
                                                         ttAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATagag                                                                                             >  4:343808/1‑149 (MQ=255)
                                                          tAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGc                                                                                            <  3:1631763/149‑1 (MQ=255)
                                                             ggATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATc                                                                                         <  4:1117821/149‑1 (MQ=255)
                                                                   ccAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGa                                                                                   >  2:1293317/1‑149 (MQ=255)
                                                                    cAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGAc                                                                                  >  1:1105689/1‑149 (MQ=255)
                                                                     aGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACa                                                                                 <  4:1489122/149‑1 (MQ=255)
                                                                      gTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACAt                                                                                <  1:580084/149‑1 (MQ=255)
                                                                       tttCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGTCCATCGTTATGGTGATAGCCGGTATGGAGCATCTg                                                                                       >  3:183793/1‑141 (MQ=255)
                                                                       tttCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACAtt                                                                               >  3:548226/1‑149 (MQ=255)
                                                                       tttCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACAtt                                                                               >  1:1117042/1‑149 (MQ=255)
                                                                          ctcGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGt                                                                            <  4:444737/149‑1 (MQ=255)
                                                                             gCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGtttt                                                                         <  2:1226392/149‑1 (MQ=255)
                                                                              ccACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTc                                                                        <  2:1459329/149‑1 (MQ=255)
                                                                              ccACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTc                                                                        <  1:1574359/149‑1 (MQ=255)
                                                                               cACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTta                                                                                                                 <  4:1630008/107‑1 (MQ=255)
                                                                               cACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCa                                                                       >  1:1662156/1‑149 (MQ=255)
                                                                               cACGACGCTGCACCACGCCTTTGGCGATTGGTAGTGTGTCCAGCGGCGTCAGTTGCACGTCGATCGAGTTGCCACGCGGGTAGCGCACCGCGGGCGGGCCGTTGTTATAGTAAAACCCTGTATAGAGCATCTGGCGACATGCGTTCTCa                                                                       <  4:361462/149‑1 (MQ=255)
                                                                                  gacgCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCg                                                                    >  3:752369/1‑149 (MQ=255)
                                                                                       ttCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCgg                                                               >  1:618419/1‑149 (MQ=255)
                                                                                       ttCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCgg                                                               <  3:1607852/149‑1 (MQ=255)
                                                                                       ttCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCgg                                                               >  3:963466/1‑149 (MQ=255)
                                                                                       ttCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCgg                                                               <  4:18761/149‑1 (MQ=255)
                                                                                         cacaATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACAGCGTTGCCACGCGGGTAGCGCACCGCTGACgg                                                                                                                           >  4:713267/1‑87 (MQ=255)
                                                                                         cacaATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACAGCGTTGCCACGCGGGTAGCGCACCGCTGACgg                                                                                                                           <  3:713267/87‑1 (MQ=255)
                                                                                          acaATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGt                                                            >  2:1429224/1‑149 (MQ=255)
                                                                                            aaTGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGc                                                                                            <  3:1338882/115‑1 (MQ=255)
                                                                                            aaTGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCa                                                          >  3:722986/1‑149 (MQ=255)
                                                                                            aaTGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCa                                                          >  3:1073968/1‑149 (MQ=255)
                                                                                            aaTGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCa                                                          >  1:804033/1‑149 (MQ=255)
                                                                                              tGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGAt                                                                                                          <  4:145823/99‑1 (MQ=255)
                                                                                              tGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGAt                                                                                                          <  4:480910/99‑1 (MQ=255)
                                                                                               gCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATaa                                                       <  3:45442/149‑1 (MQ=255)
                                                                                               gCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATaa                                                       <  4:963466/149‑1 (MQ=255)
                                                                                                 cTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGt                                                                                                                   <  4:684726/87‑1 (MQ=255)
                                                                                                  tttGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGa                                                    <  3:142928/149‑1 (MQ=255)
                                                                                                   ttGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGAc                                                   >  1:316695/1‑149 (MQ=255)
                                                                                                   ttGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGAc                                                   >  2:1312985/1‑149 (MQ=255)
                                                                                                   ttGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGAc                                                   <  4:334163/149‑1 (MQ=255)
                                                                                                    tGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGAcc                                                  <  3:1465948/149‑1 (MQ=255)
                                                                                                      ccAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCAt                                                <  1:437485/149‑1 (MQ=255)
                                                                                                       cAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCggg                                                              >  2:366323/1‑134 (MQ=255)
                                                                                                       cAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCAtt                                               >  1:476327/1‑149 (MQ=255)
                                                                                                         aTTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTc                                             >  1:1178764/1‑149 (MQ=255)
                                                                                                             gTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGt                                         <  1:37463/149‑1 (MQ=255)
                                                                                                              tAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTa                                        >  4:211086/1‑149 (MQ=255)
                                                                                                              tAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTa                                        <  1:1541289/149‑1 (MQ=255)
                                                                                                                gTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATg                                      >  2:1246912/1‑149 (MQ=255)
                                                                                                                   tttCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAg                                   <  1:1293317/149‑1 (MQ=255)
                                                                                                                      ccAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAgcgc                                >  1:682284/1‑149 (MQ=255)
                                                                                                                        aGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAg                              <  4:13665/149‑1 (MQ=255)
                                                                                                                          cggcgTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGt                            <  3:211086/149‑1 (MQ=255)
                                                                                                                            gcgTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTaa                          <  2:80067/149‑1 (MQ=255)
                                                                                                                            gcgTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTaa                          <  3:629260/149‑1 (MQ=255)
                                                                                                                              GTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAaga                        >  4:274332‑M1/1‑149 (MQ=255)
                                                                                                                              GTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAaga                        >  4:1041076‑M1/1‑149 (MQ=255)
                                                                                                                               TCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAagag                       <  1:409391‑M1/149‑1 (MQ=255)
                                                                                                                                  GTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAAgagaga                    >  3:945734‑M1/1‑149 (MQ=255)
                                                                                                                                   ttCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAAGAGAGAt                   <  4:860332/149‑1 (MQ=255)
                                                                                                                                     ccACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAgg                             >  2:1167480/1‑137 (MQ=255)
                                                                                                                                       aCGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAAGAGAGATCaaa               >  2:1485805/1‑149 (MQ=255)
                                                                                                                                       aCGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAAGAGAGATCaaa               >  3:103027/1‑149 (MQ=255)
                                                                                                                                        cGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAAGAGAGATCaaaa              <  4:1073968/149‑1 (MQ=255)
                                                                                                                                        cGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAAGAGAGATCaaaa              >  2:931580/1‑149 (MQ=255)
                                                                                                                                           tacccGCGTGGCAACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAagagag                     >  2:251747/4‑139 (MQ=255)
                                                                                                                                           tacccGCGTGGCAACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAagagag                     <  1:251747/136‑1 (MQ=255)
                                                                                                                                            acccGCGTGGCAACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAagagag                     <  2:838769/136‑1 (MQ=255)
                                                                                                                                                    tGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAAGAGAGATCAAAAGCACCCTGATgg  <  4:722986/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGCGTGGCGTTCAGCGATTCGGCGACTTTCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAAGAGAGATCAAAAGCACCCTGATGG  >  NC_000913/438592‑438888

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