Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,027,819 A→G C318C (TGT→TGC fadA ← 3‑ketoacyl‑CoA thiolase (thiolase I)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,027,8190AG97.7% 303.4 / ‑3.6 86C318C (TGT→TGCfadA3‑ketoacyl‑CoA thiolase (thiolase I)
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (43/41);  total (44/42)
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.56e-01

CAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCACAACCAACGACCGCCATC  >  NC_000913/4027671‑4027966
                                                                                                                                                    |                                                                                                                                                    
cAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCa                                                                                                                                                    <  1:1494507/151‑1 (MQ=255)
      gAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCa                                                                                                                                               <  1:606702/150‑1 (MQ=255)
       aaCGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAgg                                                                                                                                             <  1:1574821/151‑1 (MQ=255)
          gTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATc                                                                                                                                          >  1:277920/1‑151 (MQ=255)
            cTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTg                                                                                                                                        >  2:486257/1‑151 (MQ=255)
             tttGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTgc                                                                                                                                       >  1:418015/1‑151 (MQ=255)
              ttGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTgcg                                                                                                                                      >  1:101625/1‑151 (MQ=255)
              ttGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTg                                                                                                                                        <  2:268286/149‑1 (MQ=255)
                  gTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGc                                                                                                                                  <  2:1078000/151‑1 (MQ=255)
                  gTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGc                                                                                                                                  >  1:498499/1‑151 (MQ=255)
                   ttCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGc                                                                                                                                  >  1:1184870/1‑150 (MQ=255)
                   ttCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGc                                                                                                                                  >  2:874543/1‑150 (MQ=255)
                    tCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCa                                                                                                                                 >  2:1388877/1‑150 (MQ=255)
                     ccATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCaaa                                                                                                                               <  1:766454/151‑1 (MQ=255)
                     ccATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCaaa                                                                                                                               <  2:1703198/151‑1 (MQ=255)
                      cATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCaa                                                                                                                                >  1:629010/1‑149 (MQ=255)
                       aTCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGGTCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAgg                                                                                                                             >  2:1315744/1‑151 (MQ=255)
                       aTCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAgg                                                                                                                             >  2:1446280/1‑151 (MQ=255)
                           gATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTc                                                                                                                         >  2:1536371/1‑151 (MQ=255)
                           gATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTc                                                                                                                         <  2:1172359/151‑1 (MQ=255)
                           gATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTc                                                                                                                         <  1:445177/151‑1 (MQ=255)
                             ttCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGACCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGt                                                                                                                      >  1:751434/1‑151 (MQ=255)
                              tCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGtt                                                                                                                     >  2:1132474/1‑151 (MQ=255)
                              tCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGt                                                                                                                      <  2:1630988/150‑1 (MQ=255)
                              tCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGt                                                                                                                      <  1:509389/150‑1 (MQ=255)
                               cagcagCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTc                                                                                                                    >  1:856129/1‑151 (MQ=255)
                                 gcagcGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCAt                                                                                                                  <  2:742559/151‑1 (MQ=255)
                                  cagcGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCa                                                                                                                   <  1:1320504/149‑1 (MQ=255)
                                  cagcGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCa                                                                                                                   >  2:1320599/1‑149 (MQ=255)
                                    gcGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGAt                                                                                                                                           >  2:1527509/1‑124 (MQ=255)
                                    gcGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGAt                                                                                                                                           <  1:1527412/124‑1 (MQ=255)
                                      gtggtgCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCa                                                                                                              >  1:1288475/1‑150 (MQ=255)
                                        ggtgCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAac                                                                                                           <  2:1567671/151‑1 (MQ=255)
                                        ggtgCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAac                                                                                                           <  1:455774/151‑1 (MQ=255)
                                               aTACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGa                                                                                                    <  1:1634665/151‑1 (MQ=255)
                                                 aCGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGat                                                                                                   >  2:1511254/1‑150 (MQ=255)
                                                   gcgcACCGGACCACCCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATc                                                                                                <  1:1021260/151‑1 (MQ=255)
                                                   gcgcACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATc                                                                                                >  1:268313/1‑151 (MQ=255)
                                                    cgcACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCg                                                                                               >  2:30357/1‑151 (MQ=255)
                                                     gcACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCg                                                                                               <  2:176341/150‑1 (MQ=255)
                                                      cACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCt                                                                                             >  1:893216/1‑151 (MQ=255)
                                                      cACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCt                                                                                             >  1:593342/1‑151 (MQ=255)
                                                      cACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCt                                                                                             >  1:34802/1‑151 (MQ=255)
                                                        ccGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCAGGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTg                                                                                            <  2:41939/150‑1 (MQ=255)
                                                            aacaacCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAg                                                                                        <  1:1647145/150‑1 (MQ=255)
                                                             acaacCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGaa                                                                                      <  2:1623411/151‑1 (MQ=255)
                                                               aacCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGaaa                                                                                     <  1:758598/150‑1 (MQ=255)
                                                                 cccAGGGGATGACCCAGGGCGATCGCGCCACGGTTAAGGTTGATTTTCTCGTCAATCTGTTCAATTAGTCCCAGGTCTTTAATGCTTGGCAGAATGAGCGCGGCAAAGGTTCGTTTCATTCCAAACACGCCGATATCGCTGGCAGAAAGc                                                                                   <  1:1673402/150‑1 (MQ=255)
                                                                 cccAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGcc                                                                                  <  2:631205/151‑1 (MQ=255)
                                                                      cGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCt                                                                             <  2:33348/151‑1 (MQ=255)
                                                                      cGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCt                                                                             >  2:1570030/1‑151 (MQ=255)
                                                                      cGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGCCCCAGATCTTTAATGCAGGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCg                                                                               <  2:855025/149‑1 (MQ=255)
                                                                         aTGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCttt                                                                           <  2:156393/150‑1 (MQ=255)
                                                                           gACCCAGCGCGATCGCGCCCCCTTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGTATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCtttttt                                                                        <  2:377048/151‑1 (MQ=255)
                                                                            aCCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTc                                                                       >  2:809319/1‑151 (MQ=255)
                                                                             ccaaGCGCGATCGCGACACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCa                                                                      <  1:1221359/148‑1 (MQ=255)
                                                                              ccAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTc                                                                       >  2:449323/1‑149 (MQ=255)
                                                                              ccAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTACAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAg                                                                     <  1:442659/151‑1 (MQ=255)
                                                                                      accgcgCCACCGTTGAGTTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCTTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGtt                                                             <  2:1518602/149‑1 (MQ=255)
                                                                                        cgcgCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGttt                                                            <  1:579302/150‑1 (MQ=255)
                                                                                        cgcgCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGttt                                                            <  1:694047/150‑1 (MQ=255)
                                                                                          cgcCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCg                                                          >  2:1246131/1‑150 (MQ=255)
                                                                                          cgcCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGa                                                         <  1:1049171/151‑1 (MQ=255)
                                                                                            ccACCGTTGAGGTTGATTTTCCCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGCCGGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAACCACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTc                                                           <  1:912932/147‑1 (MQ=255)
                                                                                                cGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCgg                                                   >  2:1301703/1‑151 (MQ=255)
                                                                                                 gTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGa                                                  <  2:781105/151‑1 (MQ=255)
                                                                                                   tGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAAc                                                >  1:84589/1‑151 (MQ=255)
                                                                                                       gTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCggg                                            >  2:1010821/1‑151 (MQ=255)
                                                                                                       gTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCggg                                            >  2:965293/1‑151 (MQ=255)
                                                                                                          gATCTTCTCGTCAATCGGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCGGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGATGCCGGAACCGGGcc                                          <  1:883680/150‑1 (MQ=255)
                                                                                                           aTCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGt                                        >  2:1497483/1‑151 (MQ=255)
                                                                                                                 tcGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAAccc                                   >  1:1619985/1‑150 (MQ=255)
                                                                                                                      aaTCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATc                             >  2:1382997/1‑151 (MQ=255)
                                                                                                                         cTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGAt                          >  2:63787/1‑151 (MQ=255)
                                                                                                                           gCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATg                         <  1:934399/150‑1 (MQ=255)
                                                                                                                            cTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATggg                       >  2:969926/1‑151 (MQ=255)
                                                                                                                               aaTTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTc                     >  1:826517/1‑150 (MQ=255)
                                                                                                                               aaTTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTc                     >  1:1215426/1‑150 (MQ=255)
                                                                                                                                     tCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCAcaacc               >  1:885631/1‑150 (MQ=255)
                                                                                                                                      cccAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCAcaaccaa             <  2:751495/151‑1 (MQ=255)
                                                                                                                                      cccAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCAcaaccaa             <  1:1446182/151‑1 (MQ=255)
                                                                                                                                         aGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCACAACCAACg           >  1:1079267/1‑150 (MQ=255)
                                                                                                                                            tCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCACAACCAACGACCg       >  1:1234006/1‑151 (MQ=255)
                                                                                                                                               ttAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCACAACCAACGACCGcc     >  2:80710/1‑150 (MQ=255)
                                                                                                                                               ttAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCACAACCAACGACCGc      <  1:809253/149‑1 (MQ=255)
                                                                                                                                                 aaTGCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCACAACCAACGACCGCCATc  >  2:1190633/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCACAACCAACGACCGCCATC  >  NC_000913/4027671‑4027966

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