Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,903,952 C→T 6.1% G549D (GGT→GAT)  bglF ← beta‑glucoside specific PTS enzyme II/BglG kinase/BglG phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,903,9520CT 6.1% 305.9 / 6.4 98G549D (GGT→GAT) bglFbeta‑glucoside specific PTS enzyme II/BglG kinase/BglG phosphatase
Reads supporting (aligned to +/- strand):  ref base C (50/41);  new base T (4/2);  total (55/43)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.81e-01

ACTGATTAATACCGGCGTCGTCAGATCAAATCCGGCCTCGCGAATAGCAGGGATATCAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACT  >  NC_000913/3903804‑3904100
                                                                                                                                                    |                                                                                                                                                    
aCTGATTAATACCGGCGTCGTCAGATCAAATCCGGCCTCGCGAATAGCAGGGATATCAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATAcc                                                                                                                                                     >  2:147258/1‑150 (MQ=255)
 gtGATTAATACAGGAGTCGTCAGATCAATTCCGGACTAGCGAATAGCAGGGATATAAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCg                                                                                                                                                    <  2:630569/149‑1 (MQ=255)
        aTACCGGCGTCGTCAGATCAAATCCGGCCTCGCGAATAGCAGGGATATCAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGaa                                                                                                                                             >  2:1103723/1‑150 (MQ=255)
         tACCGGCGTCGTCAGATCAAATCCGGCCTCGCGAATAGCAGGGATATCAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGGTATCGACaccgta                                                                                                                                            <  2:516155/150‑6 (MQ=255)
         tACCGGCGTCGTCAGATCAAATCCGGCCTCGCGAATAGCAGGGATATCAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAAt                                                                                                                                            >  2:1328774/1‑150 (MQ=255)
           ccGGCGTCGTCAGATCAAATACTGCCTCGAGAATATCAGGGATATCAAAAGAAATCATACGATCGCCTTTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCTACATGAATCa                                                                                                                                          >  2:607895/1‑150 (MQ=255)
            cGGCGTCGTCAGATCAAATCCGGCCTCGCGAATAGCAGGGATATCAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAg                                                                                                                                         >  1:1238402/1‑150 (MQ=255)
            cGGCGTCGTCAGATCAAATCCGGCCTCGCGAATAGCAGGGATATCAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTCGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAg                                                                                                                                         <  2:1352884/150‑1 (MQ=255)
             ggCGTCGTCAGATCAAATCCGGCCTCGCGAATAGCAGGGATATCAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAgg                                                                                                                                        >  2:1135131/1‑150 (MQ=255)
              gCGTCGTCAGATCAAATCCGGCCTCGCGAATAGCAGGGATATCAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGa                                                                                                                                       <  1:538550/150‑1 (MQ=255)
                  cgtcAGATCAAATCCGGCCTCGCGAATAGCAGGGATATCAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGAtctc                                                                                                                                   >  2:49804/1‑150 (MQ=255)
                   gtcAGATCAAATCCGGCCTCGCGAATAGCAGGGATATCAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTcc                                                                                                                                  <  2:957771/150‑1 (MQ=255)
                     cAGATCAAATCCGGCCTCGCGAATAGCAGGGATATCAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCcac                                                                                                                                >  2:480729/1‑150 (MQ=255)
                      aGATCAAATCCGGCCTCGCGAATAGCAGGGATATCAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCcaca                                                                                                                               >  1:730471/1‑150 (MQ=255)
                      aGATCAAATCCGGCCTCGCGAATAGCAGGGATATCAAAAGAAATCAGCCGATCGCATCTATTGACCTTGTCACCCACGATGACGTGAGCGGAAAAGAATTTGCCGTCCAGATTTACTGTGTCGATAAATAACTGAATCAGTATCTCcacc                                                                                                                               >  2:728613/1‑149 (MQ=255)
                       gATCAAATCCGGCCTCGCGAATAGCAGGGATATCAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCacac                                                                                                                              >  1:1913/1‑150 (MQ=255)
                         tCAAATCCGGCCTCGCGAATAGCAGGGATATCAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACca                                                                                                                            <  1:1135131/150‑1 (MQ=255)
                          cAAATCCGGCCTCGCGAATAGCAGGGATATCAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCTATACCGACATGAATCATGATCTCCACACcat                                                                                                                           <  2:636591/150‑1 (MQ=255)
                           aaaTCCGGCCTCGCGAATAGCAGGGATATCAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACcatc                                                                                                                          >  1:757970/1‑150 (MQ=255)
                            aaTCCGGCCTCGCGAATAGCAGGGATATCAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGGCGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACcatca                                                                                                                         >  1:1029459/1‑150 (MQ=255)
                               ccGGCCTCGCGAATAGCAGGGATATCAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCt                                                                                                                      >  2:47883/1‑150 (MQ=255)
                                               cAGGGATATCAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATgg                                                                                                      >  2:159920/1‑150 (MQ=255)
                                                 gggATATCAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCg                                                                                                    <  1:1207176/150‑1 (MQ=255)
                                                      atCAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTaa                                                                                               >  1:1016537/1‑150 (MQ=255)
                                                        cAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAAtg                                                                                             >  2:588143/1‑150 (MQ=255)
                                                         aaaaGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAAtgt                                                                                            >  2:1511791/1‑150 (MQ=255)
                                                         aaaaGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAAtgt                                                                                            <  1:43646/150‑1 (MQ=255)
                                                          aaaGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAAtgtg                                                                                           >  2:1524445/1‑150 (MQ=255)
                                                          aaaGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAAtgtg                                                                                           >  2:1530537/1‑150 (MQ=255)
                                                           aaGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCTGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTgg                                                                                          >  2:509571/1‑150 (MQ=255)
                                                            tcaatacAGGCGATCGGCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGc                                                                                         >  2:211982/7‑150 (MQ=255)
                                                            aGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGc                                                                                         <  1:286998/150‑1 (MQ=255)
                                                             gAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCg                                                                                        >  2:178025/1‑150 (MQ=255)
                                                                aTCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGaac                                                                                     >  1:1475083/1‑150 (MQ=255)
                                                                  cAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGaacaa                                                                                   >  2:1589276/1‑150 (MQ=255)
                                                                   aGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATATCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGaacaac                                                                                  >  2:1570877/1‑150 (MQ=255)
                                                                    gCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACg                                                                                 >  2:63286/1‑150 (MQ=255)
                                                                       gATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAg                                                                              <  2:1238402/150‑1 (MQ=255)
                                                                       gATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAg                                                                              <  2:1565639/150‑1 (MQ=255)
                                                                       gATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAg                                                                              >  1:190712/1‑150 (MQ=255)
                                                                                 aTTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGAcc                                                                    <  1:1589276/150‑1 (MQ=255)
                                                                                  ttGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGAccc                                                                   <  1:159920/150‑1 (MQ=255)
                                                                                     aCCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATTAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCa                                                                <  1:597009/150‑1 (MQ=255)
                                                                                     aCCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCa                                                                >  1:168825/1‑150 (MQ=255)
                                                                                      ccTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCaa                                                               >  1:348968/1‑150 (MQ=255)
                                                                                       cTTGTCACACACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGCTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAAc                                                              >  2:882180/1‑150 (MQ=255)
                                                                                        ttGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAAcc                                                             <  1:47883/150‑1 (MQ=255)
                                                                                          gggccccACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCgg                                                           <  2:918074/146‑1 (MQ=255)
                                                                                           ccgctcACGTCAACGTGGGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGga                                                          >  2:913131/6‑150 (MQ=255)
                                                                                            cACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTTCGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCACCCGgag                                                         >  2:726731/1‑150 (MQ=255)
                                                                                              cccACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGaa                                                       >  1:340340/1‑150 (MQ=255)
                                                                                               ccACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAAc                                                      >  1:1182132/1‑150 (MQ=255)
                                                                                                cACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACg                                                     >  2:1615366/1‑150 (MQ=255)
                                                                                                   gTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCAc                                                  <  2:1475083/150‑1 (MQ=255)
                                                                                                    ttGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACt                                                 <  1:1084840/150‑1 (MQ=255)
                                                                                                    ttGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACt                                                 <  1:1346283/150‑1 (MQ=255)
                                                                                                     tGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACtt                                                <  1:1548210/150‑1 (MQ=255)
                                                                                                       aCGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCa                                              <  2:168825/150‑1 (MQ=255)
                                                                                                        cGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCAc                                             >  1:1045143/1‑150 (MQ=255)
                                                                                                        cGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCAc                                             <  1:211982/150‑1 (MQ=255)
                                                                                                             gCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAAcc                                        >  1:164225/1‑150 (MQ=255)
                                                                                                             gCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAAcc                                        >  1:1471137/1‑150 (MQ=255)
                                                                                                             gCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAAcc                                        >  1:405388/1‑150 (MQ=255)
                                                                                                               ggAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGa                                      >  1:799810/1‑150 (MQ=255)
                                                                                                                gAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAg                                     >  2:1336310/1‑150 (MQ=255)
                                                                                                                gAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAg                                     >  2:272236/1‑150 (MQ=255)
                                                                                                                 aaaaGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTTACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAGGCAATTCGACCCGCAACCGGAGAACGCACTTCACCGACCGAgg                                    <  2:757970/150‑1 (MQ=255)
                                                                                                                  aaaGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAggg                                   >  2:931939/1‑150 (MQ=255)
                                                                                                                  aaaGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAggg                                   >  2:1132617/1‑150 (MQ=255)
                                                                                                                   aaGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGc                                  <  2:1016537/150‑1 (MQ=255)
                                                                                                                     gAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAg                                <  1:1103723/150‑1 (MQ=255)
                                                                                                                       aTTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCGTCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGaa                              <  2:730471/150‑1 (MQ=255)
                                                                                                                       aTTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGaa                              >  1:490316/1‑150 (MQ=255)
                                                                                                                         ttGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATg                            <  1:480729/150‑1 (MQ=255)
                                                                                                                           gCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGc                          <  1:1328774/150‑1 (MQ=255)
                                                                                                                            ccGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCa                         <  1:1222067/150‑1 (MQ=255)
                                                                                                                            ccGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCa                         <  2:1471137/150‑1 (MQ=255)
                                                                                                                            ccGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCa                         <  1:147258/150‑1 (MQ=255)
                                                                                                                              gTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAAt                       >  2:766290/1‑150 (MQ=255)
                                                                                                                               tCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATa                      <  2:1182132/150‑1 (MQ=255)
                                                                                                                                      gattcatgtCGGTATCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTAc               >  1:1505402/7‑150 (MQ=255)
                                                                                                                                      gattcatgtCGGTATCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTAc               >  1:1505403/7‑150 (MQ=255)
                                                                                                                                      gattcatgtCGGTATCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTAc               >  1:229265/7‑150 (MQ=255)
                                                                                                                                        aCGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTAccc             <  1:178025/150‑1 (MQ=255)
                                                                                                                                           gtgtCGATACCGACATGAATCAGTATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAAc          <  1:728613/150‑1 (MQ=255)
                                                                                                                                           gtgtCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAAc          <  1:412097/150‑1 (MQ=255)
                                                                                                                                           gtgtCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAAc          <  1:1106226/150‑1 (MQ=255)
                                                                                                                                           gtgtCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAAc          <  2:190712/150‑1 (MQ=255)
                                                                                                                                            tgtCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACa         >  1:481954/1‑150 (MQ=255)
                                                                                                                                              tCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGc                          >  1:1125753/1‑131 (MQ=255)
                                                                                                                                              tCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGc                          <  2:1125753/131‑1 (MQ=255)
                                                                                                                                              tCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAgg       >  2:19822/1‑150 (MQ=255)
                                                                                                                                               cGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGc      >  1:1601777/1‑150 (MQ=255)
                                                                                                                                                gATATCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGcc     <  1:1570877/150‑1 (MQ=255)
                                                                                                                                                gATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGcc     <  1:1132617/150‑1 (MQ=255)
                                                                                                                                                gATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGcc     >  1:193230/1‑150 (MQ=255)
                                                                                                                                                  tACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGCCCCGCAACCGGAGAACGCCCTTCACCAACCGCGGGCGGAATGGCAATACCTTTACCCAACCGGcccc   >  2:1444376/1‑148 (MQ=255)
                                                                                                                                                  tACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCAc   <  1:63286/150‑1 (MQ=255)
                                                                                                                                                   aCCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACt  <  1:1615366/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACTGATTAATACCGGCGTCGTCAGATCAAATCCGGCCTCGCGAATAGCAGGGATATCAAAAGAAATCAGCCGATCGCCTGTATTGACCTTGTCACCCACGTTGACGTGAGCGGAAAAGAATTTGCCGTCCAGTTTTACGGTGTCGATACCGACATGAATCAGGATCTCCACACCATCATCTGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACT  >  NC_000913/3903804‑3904100

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