Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,457,131 T→A 6.8% I677L (ATA→TTA)  fadJ ← 3‑hydroxyacyl‑CoA dehydrogenase FadJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,457,1310TA 6.8% 368.3 / 12.2 118I677L (ATA→TTA) fadJ3‑hydroxyacyl‑CoA dehydrogenase FadJ
Reads supporting (aligned to +/- strand):  ref base T (61/49);  new base A (4/4);  total (65/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

Only 100 of 124 total aligned reads displayed.

TTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTC  >  NC_000913/2456984‑2457279
                                                                                                                                                   |                                                                                                                                                    
ttCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTg                                                                                                                                                    <  1:1028769/150‑1 (MQ=255)
 tCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGc                                                                                                                                                   >  1:574910/1‑150 (MQ=255)
   aGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCaa                                                                                                                                                 <  1:922376/150‑1 (MQ=255)
    gcgcGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAAc                                                                                                                                                >  1:1629205/1‑150 (MQ=255)
     cgcgGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAAcc                                                                                                                                               <  2:144200/150‑1 (MQ=255)
     cgcgGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAAcc                                                                                                                                               <  1:1009707/150‑1 (MQ=255)
      gcgGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCa                                                                                                                                              >  1:214075/1‑150 (MQ=255)
       cgGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCAATAATGCAACGAc                                                                                                                                             <  2:1605065/150‑1 (MQ=255)
         gATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACtt                                                                                                                                           >  1:735008/1‑150 (MQ=255)
         gATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATAATGCGTGGCAAGTCGTTGAATTATTGCACCCAAtt                                                                                                                                           >  2:582527/1‑150 (MQ=255)
           ttCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCg                                                                                                                                         <  1:1045803/150‑1 (MQ=255)
             cATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGcc                                                                                                                                       >  2:1317242/1‑150 (MQ=255)
              atatAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCAATAATGCAACGAc                                                                                                                                             <  1:769911/143‑1 (MQ=255)
              atatAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCAATAATGCAACGAc                                                                                                                                             >  2:769911/1‑143 (MQ=255)
               tataGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCcg                                                                                                                                     >  1:1517133/1‑150 (MQ=255)
                  aGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCgcgc                                                                                                                                  <  1:114060/150‑1 (MQ=255)
                   gCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGcc                                                                                                                                 >  2:967063/1‑150 (MQ=255)
                       tGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCgaga                                                                                                                             <  2:1629205/150‑1 (MQ=255)
                        gACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCgagag                                                                                                                            <  1:368321/150‑1 (MQ=255)
                        gACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCgagag                                                                                                                            <  2:717935/150‑1 (MQ=255)
                         aCCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCgagaga                                                                                                                           <  1:1332465/150‑1 (MQ=255)
                          ccTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCAATAATGCAACGACttgccacg                                                                                                                                     <  1:241680/139‑7 (MQ=255)
                          ccTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCAATAATGCAACGACttgccacg                                                                                                                                     >  2:241680/1‑133 (MQ=255)
                           cttcttATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAAt                                                                                                                         <  1:1096204/150‑1 (MQ=255)
                            ttcttATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGAGTAAAACGGGAACCATACTGCGTGGCAAGTGGTTGCAATAATGCAACGACttgcc                                                                                                                                        <  1:1029479/134‑4 (MQ=255)
                            ttcttATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGAGTAAAACGGGAACCATACTGCGTGGCAAGTGGTTGCAATAATGCAACGACttgcc                                                                                                                                        >  2:1029479/1‑131 (MQ=255)
                             tcttATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCg                                                                                                                       >  1:1347205/1‑150 (MQ=255)
                              cttATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGa                                                                                                                      >  2:803171/1‑150 (MQ=255)
                              cttATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGa                                                                                                                      <  1:276387/150‑1 (MQ=255)
                               ttATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGat                                                                                                                     <  1:68330/150‑1 (MQ=255)
                                   tGCGGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGGGCCCATCTCGAGCAAACGCTCGCAAGGGGTAAAACGGGAGCCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCCCTTCGCCCGCGCCGAGAGAATCGCTTTAg                                                                                                                 >  2:719679/1‑150 (MQ=255)
                                   tGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAg                                                                                                                 >  2:160638/1‑150 (MQ=255)
                                    gCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGc                                                                                                                >  2:1094272/1‑150 (MQ=255)
                                      aGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCgg                                                                                                              <  2:1421878/150‑1 (MQ=255)
                                         tCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCTCTTCGCCCGCGCCGAGAGAATCGATATAGCGGAAc                                                                                                           >  2:1197339/1‑150 (MQ=255)
                                           agttgcagttgTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACgg                                                                                                         >  2:1411228/1‑150 (MQ=255)
                                           agttgcagttgTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACgg                                                                                                         >  1:838539/1‑150 (MQ=255)
                                                      ttttCCAAAAACTTTCCCCCCGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGACCCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGaa                                                                                              >  2:1570047/1‑150 (MQ=255)
                                                      ttttACAAAAACATACCCCACGAGCGCACATCTCGACCAAACGCTAGCAAGCGGTAAAACGGGAAACATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGaa                                                                                              >  2:393202/1‑150 (MQ=255)
                                                         tCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATg                                                                                           >  1:426074/1‑150 (MQ=255)
                                                         tCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATg                                                                                           <  1:1309751/150‑1 (MQ=255)
                                                          ccAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATgg                                                                                          >  1:994211/1‑150 (MQ=255)
                                                          ccAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATgg                                                                                          >  2:821367/1‑150 (MQ=255)
                                                           cAAAAACTTTCCCCACGCGCGCCCATCTCGCCCACACGCTCGCAAGGGGTAAAACGGGACCCATACTGCGTGGCAAGTCGTTGCATTATTGCACCCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCCCCGAGAAATGGc                                                                                         >  2:1061675/1‑150 (MQ=255)
                                                           cAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGACCCATACTGCGTGGCAAGTCGTTGCATTATTGCACCCACTTCGCCCGCGCCGAGAGATTCGATATAGCGGGACGGTCCACCGAGAACTggg                                                                                         >  2:1109627/1‑149 (MQ=255)
                                                           cAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGc                                                                                         <  1:160638/150‑1 (MQ=255)
                                                           cAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGc                                                                                         >  1:1221092/1‑150 (MQ=255)
                                                            aaaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCg                                                                                        <  1:776358/150‑1 (MQ=255)
                                                             aaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCgg                                                                                       <  1:1317242/150‑1 (MQ=255)
                                                             aaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCgg                                                                                       >  1:1193615/1‑150 (MQ=255)
                                                             aaaaCTTTCCCCACGCGCGCCCAACTCGACCAACCGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATAAAGCGGAACGGTCCACCTAGAAATGGCgg                                                                                       >  2:609017/1‑150 (MQ=255)
                                                                aCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGaaa                                                                                    <  1:505779/150‑1 (MQ=255)
                                                                 cTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGATGATCGGTCCACCGAGAAATGGCGGaaaa                                                                                   <  2:479687/150‑1 (MQ=255)
                                                                  tttCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAAc                                                                                  >  1:650961/1‑150 (MQ=255)
                                                                      cccACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAAt                                                                              >  1:954937/1‑150 (MQ=255)
                                                                       ccACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATg                                                                             >  2:92416/1‑150 (MQ=255)
                                                                       ccACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATg                                                                             <  2:214075/150‑1 (MQ=255)
                                                                       ccACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATg                                                                             >  1:200375/1‑150 (MQ=255)
                                                                             gcgcCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAAt                                                                       <  2:1427500/150‑1 (MQ=255)
                                                                                  cATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACcgc                                                                  <  1:1587039/150‑1 (MQ=255)
                                                                                     ctcGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTGGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGcc                                                               <  2:465533/150‑1 (MQ=255)
                                                                                      tcGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCa                                                              >  2:781431/1‑150 (MQ=255)
                                                                                      tcGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCa                                                              >  1:103000/1‑150 (MQ=255)
                                                                                       cGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCaa                                                             >  1:201154/1‑150 (MQ=255)
                                                                                         aCCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAata                                                           <  1:744273/150‑1 (MQ=255)
                                                                                            aaaCGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATcc                                                        <  1:217746/150‑1 (MQ=255)
                                                                                            aaaCGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATcc                                                        >  2:652989/1‑150 (MQ=255)
                                                                                            aaaCGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATcc                                                        >  1:1501916/1‑150 (MQ=255)
                                                                                               cGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCg                                                     <  2:1513919/150‑1 (MQ=255)
                                                                                               cGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCg                                                     <  2:1540916/150‑1 (MQ=255)
                                                                                               cGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCg                                                     <  1:1108789/150‑1 (MQ=255)
                                                                                               cGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGCATGGGGGAAAACCCATGCCAAGTACCGCGCCGAGATCCCCgg                                                    >  2:877891/1‑149 (MQ=255)
                                                                                                   cGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTcac                                                 >  2:930117/1‑150 (MQ=255)
                                                                                                     cAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTcacgc                                               <  1:381255/150‑1 (MQ=255)
                                                                                                         gggTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCacgcacgc                                           >  2:574264/1‑150 (MQ=255)
                                                                                                          ggTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCt                                          >  1:1383861/1‑150 (MQ=255)
                                                                                                             aaaaCGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACg                                       >  2:1180800/1‑150 (MQ=255)
                                                                                                             aaaaCGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTACGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACg                                       >  1:568810/1‑150 (MQ=255)
                                                                                                              aaaCGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCACATACCGCGCCAATATCCCCGTCCCGCACGCTACgg                                      >  2:412824/1‑150 (MQ=255)
                                                                                                              aaaCGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACgg                                      <  2:838539/150‑1 (MQ=255)
                                                                                                               aaCGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGa                                     <  2:1517133/150‑1 (MQ=255)
                                                                                                                aCGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGAt                                    >  1:75749/1‑150 (MQ=255)
                                                                                                                      aCCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTg                              >  1:1380794/1‑150 (MQ=255)
                                                                                                                          tACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCtca                          <  1:860691/150‑1 (MQ=255)
                                                                                                                              gCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCaa                      >  1:1407213/1‑150 (MQ=255)
                                                                                                                              gCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCaa                      >  2:776306/1‑150 (MQ=255)
                                                                                                                               cGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAac                     <  1:1411228/150‑1 (MQ=255)
                                                                                                                               cGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAac                     <  2:426074/150‑1 (MQ=255)
                                                                                                                                gTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGTGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAaca                    <  2:340848/150‑1 (MQ=255)
                                                                                                                                gTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAaca                    <  1:967063/150‑1 (MQ=255)
                                                                                                                                gTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAaca                    <  2:463159/150‑1 (MQ=255)
                                                                                                                                  ggCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAAcaca                  >  1:392701/1‑150 (MQ=255)
                                                                                                                                          gTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTg          >  1:1333505/1‑150 (MQ=255)
                                                                                                                                              cATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCttc      >  1:1327889/1‑150 (MQ=255)
                                                                                                                                               attattGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCttca     >  1:493553/1‑150 (MQ=255)
                                                                                                                                               attattGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCttca     <  2:398358/150‑1 (MQ=255)
                                                                                                                                                  attGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAAcc                                <  2:1206094/120‑1 (MQ=255)
                                                                                                                                                  attGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAAcc                                >  1:1206094/1‑120 (MQ=255)
                                                                                                                                                  attGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCttcattc  >  1:1492633/1‑150 (MQ=255)
                                                                                                                                                  attGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCttcattc  >  1:1516455/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTC  >  NC_000913/2456984‑2457279

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