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→TGC) fadA3‑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

CAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCACAACCAACGACCGCCATC  >  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)
                             ttCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGACCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑t                                                                                                                      >  1:751434/1‑151 (MQ=255)
                              tCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑tt                                                                                                                     >  2:1132474/1‑151 (MQ=255)
                              tCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑t                                                                                                                      <  2:1630988/150‑1 (MQ=255)
                              tCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑t                                                                                                                      <  1:509389/150‑1 (MQ=255)
                               cagcagCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTc                                                                                                                    >  1:856129/1‑151 (MQ=255)
                                 gcagcGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCAt                                                                                                                  <  2:742559/151‑1 (MQ=255)
                                  cagcGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCa                                                                                                                   <  1:1320504/149‑1 (MQ=255)
                                  cagcGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCa                                                                                                                   >  2:1320599/1‑149 (MQ=255)
                                    gcGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGAt                                                                                                                                           >  2:1527509/1‑124 (MQ=255)
                                    gcGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGAt                                                                                                                                           <  1:1527412/124‑1 (MQ=255)
                                      gtggtgCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCa                                                                                                              >  1:1288475/1‑150 (MQ=255)
                                        ggtgCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAac                                                                                                           <  2:1567671/151‑1 (MQ=255)
                                        ggtgCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAac                                                                                                           <  1:455774/151‑1 (MQ=255)
                                               aTACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGa                                                                                                    <  1:1634665/151‑1 (MQ=255)
                                                 aCGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGat                                                                                                   >  2:1511254/1‑150 (MQ=255)
                                                   gcgcACCGGACCACCCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATc                                                                                                <  1:1021260/151‑1 (MQ=255)
                                                   gcgcACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATc                                                                                                >  1:268313/1‑151 (MQ=255)
                                                    cgcACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCg                                                                                               >  2:30357/1‑151 (MQ=255)
                                                     gcACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCg                                                                                               <  2:176341/150‑1 (MQ=255)
                                                      cACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCt                                                                                             >  1:893216/1‑151 (MQ=255)
                                                      cACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCt                                                                                             >  1:593342/1‑151 (MQ=255)
                                                      cACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCt                                                                                             >  1:34802/1‑151 (MQ=255)
                                                        ccGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCAGGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTg                                                                                            <  2:41939/150‑1 (MQ=255)
                                                            aacaacCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAg                                                                                        <  1:1647145/150‑1 (MQ=255)
                                                             acaacCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGaa                                                                                      <  2:1623411/151‑1 (MQ=255)
                                                               aacCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGaaa                                                                                     <  1:758598/150‑1 (MQ=255)
                                                                 cccAGGGGATGACCCAGGGCGATCGCGCCACGGTTAAGGTTGATTTTCTCGTCAATCTGTTCAATTAGTCCCAGGTCTTTAATGCTTGGCAGAATGAGCGCGGCAAAGG‑TTCGTTTCATTCCAAACACGCCGATATCGCTGGCAGAAAGc                                                                                   <  1:1673402/150‑1 (MQ=255)
                                                                 cccAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGcc                                                                                  <  2:631205/151‑1 (MQ=255)
                                                                      cGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCt                                                                             <  2:33348/151‑1 (MQ=255)
                                                                      cGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCt                                                                             >  2:1570030/1‑151 (MQ=255)
                                                                      cGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGCCCCAGATCTTTAATGCAGGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCg                                                                               <  2:855025/149‑1 (MQ=255)
                                                                         aTGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCttt                                                                           <  2:156393/150‑1 (MQ=255)
                                                                           gACCCAGCGCGATCGCGCCCCCTTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGTATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCtttttt                                                                        <  2:377048/151‑1 (MQ=255)
                                                                            aCCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTc                                                                       >  2:809319/1‑151 (MQ=255)
                                                                             ccaaGCGCGATCGCGACACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCa                                                                      <  1:1221359/148‑1 (MQ=255)
                                                                              ccAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTc                                                                       >  2:449323/1‑149 (MQ=255)
                                                                              ccAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTACAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAg                                                                     <  1:442659/151‑1 (MQ=255)
                                                                                      accgcgCCACCGTTGAGTTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCT‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGtt                                                             <  2:1518602/149‑1 (MQ=255)
                                                                                        cgcgCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGttt                                                            <  1:579302/150‑1 (MQ=255)
                                                                                        cgcgCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGttt                                                            <  1:694047/150‑1 (MQ=255)
                                                                                          cgcCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCg                                                          >  2:1246131/1‑150 (MQ=255)
                                                                                          cgcCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGa                                                         <  1:1049171/151‑1 (MQ=255)
                                                                                            ccACCGTTGAGGTTGATTTTCCCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGCCGGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAACCACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTc                                                           <  1:912932/147‑1 (MQ=255)
                                                                                                cGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCgg                                                   >  2:1301703/1‑151 (MQ=255)
                                                                                                 gTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGa                                                  <  2:781105/151‑1 (MQ=255)
                                                                                                   tGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAAc                                                >  1:84589/1‑151 (MQ=255)
                                                                                                       gTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCggg                                            >  2:1010821/1‑151 (MQ=255)
                                                                                                       gTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCggg                                            >  2:965293/1‑151 (MQ=255)
                                                                                                          gATCTTCTCGTCAATCGGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCGGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGATGCCGGAACCGGGcc                                          <  1:883680/150‑1 (MQ=255)
                                                                                                           aTCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGt                                        >  2:1497483/1‑151 (MQ=255)
                                                                                                                 tcGTCAATCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAAccc                                   >  1:1619985/1‑150 (MQ=255)
                                                                                                                      aaTCTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATc                             >  2:1382997/1‑151 (MQ=255)
                                                                                                                         cTGCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGAt                          >  2:63787/1‑151 (MQ=255)
                                                                                                                           gCTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATg                         <  1:934399/150‑1 (MQ=255)
                                                                                                                            cTCAATTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATggg                       >  2:969926/1‑151 (MQ=255)
                                                                                                                               aaTTAGTCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTc                     >  1:826517/1‑150 (MQ=255)
                                                                                                                               aaTTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTc                     >  1:1215426/1‑150 (MQ=255)
                                                                                                                                     tCCCAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCAcaacc               >  1:885631/1‑150 (MQ=255)
                                                                                                                                      cccAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCAcaaccaa             <  2:751495/151‑1 (MQ=255)
                                                                                                                                      cccAGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCAcaaccaa             <  1:1446182/151‑1 (MQ=255)
                                                                                                                                         aGATCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCACAACCAACg           >  1:1079267/1‑150 (MQ=255)
                                                                                                                                            tCTTTAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCACAACCAACGACCg       >  1:1234006/1‑151 (MQ=255)
                                                                                                                                               ttAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCACAACCAACGACCGcc     >  2:80710/1‑150 (MQ=255)
                                                                                                                                               ttAATGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCACAACCAACGACCGc      <  1:809253/149‑1 (MQ=255)
                                                                                                                                                 aaTGCATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCACAACCAACGACCGCCATc  >  2:1190633/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCG‑TTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCACAACCAACGACCGCCATC  >  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: ‑