Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,028,482 T→G 9.6% A97A (GCA→GCC) ‡ fadA ← 3‑ketoacyl‑CoA thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,028,4820TG 9.6% 264.7 / 19.7 94A97A (GCA→GCC) ‡fadA3‑ketoacyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base T (36/49);  new base G (4/5);  total (40/54)
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.22e-01

TGCCCGCCGCTTTGGCGACATTGCGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCG  >  NC_000913/4028337‑4028630
                                                                                                                                                 |                                                                                                                                                    
tGCCCGCCGCTTTGGCGACATTGCGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCt                                                                                                                                                  <  2:213599/150‑1 (MQ=255)
tGCCCGCCGCTTTGGCGACATTGCGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCt                                                                                                                                                  <  1:880/150‑1 (MQ=255)
  cccGCCGCTTTGGCGACATTGCGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGc                                                                                                                                                <  1:753019/150‑1 (MQ=255)
     gccgcTTTGGCGACATTGCGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATg                                                                                                                                             <  2:394317/150‑1 (MQ=255)
      ccgcTTTGGCGACATTGCGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACGAGACATGCCTGCGCATCGCCAGTGATGATCATTCGTGATGCGTCATGCAGTGCCTGCATgg                                                                                                                                            <  2:1190803/150‑1 (MQ=255)
       cgcTTTGGCGACATTGCGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGGCACTGCATgac                                                                                                                                           <  1:1566391/150‑3 (MQ=255)
       cgcTTTGGCGACATTGCGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCAGCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGGCACTGCATgac                                                                                                                                           <  1:1514706/150‑3 (MQ=255)
         cTTTGGCGACATTGCGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATg                                                                                                                                         <  1:506941/150‑1 (MQ=255)
          tttGGCGACATTGCGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGGCACTGCATgacgca                                                                                                                                        <  1:960416/150‑6 (MQ=255)
            tGGCGACATTGCGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAAc                                                                                                                                      <  2:209965/150‑1 (MQ=255)
            tGGCGACATTGCGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGGCACTGCATgacgcagc                                                                                                                                      <  2:174710/150‑8 (MQ=255)
               cGACATTGCGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACcac                                                                                                                                   <  1:122242/150‑1 (MQ=255)
               cGACATTGCGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACcac                                                                                                                                   <  2:843801/150‑1 (MQ=255)
               cGACATTGCGGCTCAGGCCGGAGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACcac                                                                                                                                   <  1:993343/150‑1 (MQ=255)
                gACATTGCGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCTTCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACcaca                                                                                                                                  >  2:104877/1‑150 (MQ=255)
                gACATTGCGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACcaca                                                                                                                                  >  1:182838/1‑150 (MQ=255)
                  cATTGCGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATg                                                                                                                                         >  1:1475949/1‑141 (MQ=255)
                  cATTGCGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATg                                                                                                                                         <  2:1475949/141‑1 (MQ=255)
                  cATTGCGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACcacaca                                                                                                                                >  2:437237/1‑150 (MQ=255)
                   aTTGCGGCTCAGGCCGGGGTGAAAATCGACGCCCTAACTAATCTGAACATGAACCATATGCTCCACCCCGCCAACCAGACATTCCAGCGCATCGCAAGTCATGATCATTCATGCAGAGTAATTCAGTTCCTCCATATATGAACCACACaa                                                                                                                               >  2:714336/1‑150 (MQ=255)
                    ttGCGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAg                                                                                                                              >  1:815445/1‑150 (MQ=255)
                    ttGCGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAg                                                                                                                              <  2:1029482/150‑1 (MQ=255)
                       cGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCTTCATGCAGTGCCTGCATTGATGAACCACACAAGCGa                                                                                                                           <  2:743586/150‑1 (MQ=255)
                       cGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGa                                                                                                                           >  2:513045/1‑150 (MQ=255)
                       cGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGa                                                                                                                           >  2:689207/1‑150 (MQ=255)
                       cGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGa                                                                                                                           <  2:1113876/150‑1 (MQ=255)
                       cGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGa                                                                                                                           >  1:1448990/1‑150 (MQ=255)
                        ggCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGAt                                                                                                                          >  2:1234650/1‑150 (MQ=255)
                        gccTCAGGCAGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCAAATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGAt                                                                                                                          >  2:1616037/3‑150 (MQ=255)
                          cTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTa                                                                                                                        <  2:557739/150‑1 (MQ=255)
                            cAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAAc                                                                                                                      <  1:428306/150‑1 (MQ=255)
                            cAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAAc                                                                                                                      >  2:786927/1‑150 (MQ=255)
                             aGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACg                                                                                                                     >  2:1030381/1‑150 (MQ=255)
                              ggCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACgg                                                                                                                    >  1:1356837/1‑150 (MQ=255)
                                 cGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCACGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTaa                                                                                                                 <  1:1470544/150‑1 (MQ=255)
                                   gggTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAAcc                                                                                                               >  2:407485/1‑150 (MQ=255)
                                   gggTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCCGCATGGATGAACCACACAAGCGATTAACGGTAAcc                                                                                                               >  2:667542/1‑150 (MQ=255)
                                    ggTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAAccg                                                                                                              >  1:1544890/1‑150 (MQ=255)
                                       gAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAAccgccg                                                                                                           >  1:615755/1‑150 (MQ=255)
                                       gAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAAccgccg                                                                                                           >  1:1537250/1‑150 (MQ=255)
                                                  ccGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAgtgt                                                                                                <  1:714336/150‑1 (MQ=255)
                                                     tGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGt                                                                                             <  1:217823/150‑1 (MQ=255)
                                                            tCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTg                                                                                      <  2:5230/150‑1 (MQ=255)
                                                             cGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGGAGTGCCTGCATGGATGAACCACACAAGCg                                                                                                                            <  1:1633709/111‑1 (MQ=255)
                                                             cGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCg                                                                                                                            >  2:1633709/1‑111 (MQ=255)
                                                              ggCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGcc                                                                                    <  2:525535/150‑1 (MQ=255)
                                                                cacaTGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCcag                                                                                  >  1:123091/1‑150 (MQ=255)
                                                                  caTGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCcagca                                                                                >  1:1531640/1‑150 (MQ=255)
                                                                    tGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTCCTGCCagcagc                                                                              <  1:1430575/150‑1 (MQ=255)
                                                                     ggCCCATATGCTCCACGCCGCCAACCCGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAgcg                                                                             >  2:1479392/1‑150 (MQ=255)
                                                                       cccATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGcgcc                                                                           >  1:723469/1‑150 (MQ=255)
                                                                          atatGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCcgcg                                                                        <  1:412319/150‑1 (MQ=255)
                                                                           tatGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGt                                                                       <  2:127424/150‑1 (MQ=255)
                                                                            atGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGtt                                                                      <  2:479614/150‑1 (MQ=255)
                                                                              gCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTAc                                                                    <  2:104060/150‑1 (MQ=255)
                                                                              gCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTAc                                                                    >  2:617857/1‑150 (MQ=255)
                                                                               cTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACg                                                                   >  2:909402/1‑150 (MQ=255)
                                                                                tCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACgg                                                                  <  1:20409/150‑1 (MQ=255)
                                                                                 ccACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACggg                                                                 <  2:123091/150‑1 (MQ=255)
                                                                                   aCGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCg                                                               <  2:1531640/150‑1 (MQ=255)
                                                                                   aCGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCg                                                               <  1:1616037/150‑1 (MQ=255)
                                                                                         ccAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTa                                                         <  2:815445/150‑1 (MQ=255)
                                                                                               aGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAAccc                                                   >  1:1547995/1‑150 (MQ=255)
                                                                                                        tGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCt                                               >  2:70361/1‑145 (MQ=255)
                                                                                                        tGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCt                                               <  1:70361/145‑1 (MQ=255)
                                                                                                         gcgcATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCg                                         <  1:1512390/150‑1 (MQ=255)
                                                                                                              tCGACAGTCATGATCATTCATACTGATTAATCCAGTACCTCAATGGATTAACCAAAAAATAAATTAAAAATAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTctgc                                    >  2:889954/1‑150 (MQ=255)
                                                                                                                 ccAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCtgctgc                                 <  2:210841/150‑1 (MQ=255)
                                                                                                                   aGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACGAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGcac                               <  2:1052229/150‑1 (MQ=255)
                                                                                                                   aGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGcac                               <  1:1587427/150‑1 (MQ=255)
                                                                                                                    gTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGcaca                              <  2:1448990/150‑1 (MQ=255)
                                                                                                                        tGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGCACACAAc                          >  2:728180/1‑150 (MQ=255)
                                                                                                                          aTCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGCACACAAccc                        <  1:437237/150‑1 (MQ=255)
                                                                                                                            cATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGCACACAACCCCa                      <  1:889752/150‑1 (MQ=255)
                                                                                                                              ttCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGCACACAACCCCAGt                    <  1:786927/150‑1 (MQ=255)
                                                                                                                               tCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGCACACAACCCAAGTa                   <  1:1175917/150‑1 (MQ=255)
                                                                                                                                cGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTaa                  <  2:952758/150‑1 (MQ=255)
                                                                                                                                 gTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATGTTAAAACCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTaaa                 <  1:104877/150‑1 (MQ=255)
                                                                                                                                 gTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTaaa                 >  2:1651105/1‑150 (MQ=255)
                                                                                                                                  tgctgcGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAat                >  2:123122/1‑150 (MQ=255)
                                                                                                                                   gttcatccATGCAGGCACTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAata               >  1:687098/8‑150 (MQ=255)
                                                                                                                                    ttcatccATGCAGGCACTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAatat              >  2:1227685/7‑150 (MQ=255)
                                                                                                                                    ctgcGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAatat              >  1:795000/1‑150 (MQ=255)
                                                                                                                                      gcGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATAtcg            <  1:1280319/150‑1 (MQ=255)
                                                                                                                                      catccATGCAGGCACTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGcac                               >  1:1046953/5‑131 (MQ=255)
                                                                                                                                      catccATGCAGGCACTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGcac                               <  2:1046953/127‑1 (MQ=255)
                                                                                                                                        gTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATAtcgtc          >  1:531026/1‑150 (MQ=255)
                                                                                                                                         tCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATAtcgtcg         <  2:182838/150‑1 (MQ=255)
                                                                                                                                         ccATGCAGGCACTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATAtcgtcg         >  2:462939/2‑150 (MQ=255)
                                                                                                                                          cATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGa        <  1:1234650/150‑1 (MQ=255)
                                                                                                                                            tGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAgg      <  1:909402/150‑1 (MQ=255)
                                                                                                                                            tGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAgg      >  1:672/1‑150 (MQ=255)
                                                                                                                                              cAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGgc    <  1:407485/150‑1 (MQ=255)
                                                                                                                                                gTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGgccg  >  2:897183/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TGCCCGCCGCTTTGGCGACATTGCGGCTCAGGCCGGGGTGAAAATCGACGCCGTGACTCATCGGCACATGGCCCATATGCTCCACGCCGCCAACCAGACATGCCTGCGCATCGCCAGTCATGATCATTCGTGCTGCGTCATGCAGTGCCTGCATGGATGAACCACACAAGCGATTAACGGTAACCGCCGGGACAGAGTGTGGTACTTCTGCCAGCAGCGCCGCGTTACGGGCGATATTAAAACCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCG  >  NC_000913/4028337‑4028630

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: ‑