Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,136,732 A→T 5.5% L21Q (CTG→CAG)  fadE ← medium‑long chain acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,136,7320AT 5.5% 339.4 / 7.9 111L21Q (CTG→CAG) fadEmedium‑long chain acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (55/48);  new base T (3/3);  total (60/51)
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.34e-01

Only 100 of 114 total aligned reads displayed.

GTTTGCGTCGCCACTGCGGGAGCGCAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCTGTGAACGGCATGTTCAGAACAGTGACAGG  >  NC_002947/2136583‑2136881
                                                                                                                                                     |                                                                                                                                                     
gTTTGCGTCGCCACTGCGGGAGCGCAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCa                                                                                                                                                       <  2:1470884/150‑1 (MQ=255)
gTTTGCGTCGCCACTGCGGGAGCGCAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCa                                                                                                                                                       <  2:1470869/150‑1 (MQ=255)
 tttGCGTCGCCACTGCGGGAGCGCAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAg                                                                                                                                                      >  1:410844/1‑150 (MQ=255)
      gTCGCCACTGCGGGAGCGCAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGAACCTGCTGACGACATATGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCTGGCGcc                                                                                                                                                 >  2:62872/1‑150 (MQ=255)
         gCCACTGCGGGAGCGCAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGt                                                                                                                                              <  2:330978/150‑1 (MQ=255)
            aCTGCGGGAGCGCAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGCAGATCGTCAAGGGCGCCTGGCGcccttgca                                                                                                                                           <  1:62872/150‑7 (MQ=255)
               gCGGGAGCGCAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCACCAGCACATAGGCCGCCCCGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGa                                                                                                                                        <  1:1355330/150‑1 (MQ=255)
                 gggAGCGCAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAgg                                                                                                                                      <  1:880907/150‑1 (MQ=255)
                 gggAGCGCAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAgg                                                                                                                                      <  1:2293574/150‑1 (MQ=255)
                 gggAGCGCAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAgg                                                                                                                                      >  2:2293574/1‑150 (MQ=255)
                      cgcAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTAcgc                                                                                                                                 <  2:410844/150‑1 (MQ=255)
                      cgcAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTAcgc                                                                                                                                 <  2:1578510/150‑1 (MQ=255)
                       gcAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTAcgcc                                                                                                                                >  1:577450/1‑150 (MQ=255)
                        cAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTAcgccg                                                                                                                               >  2:1894175/1‑150 (MQ=255)
                         aaCCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTAcgccgc                                                                                                                              <  1:1098357/150‑1 (MQ=255)
                           ccAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTAcgccgcgc                                                                                                                            >  1:1649972/1‑150 (MQ=255)
                           ccAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTAcgccgcgc                                                                                                                            <  1:1894175/150‑1 (MQ=255)
                                     cTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAg                                                                                                                  >  1:583846/1‑150 (MQ=255)
                                      ttCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGc                                                                                                                 >  1:785843/1‑150 (MQ=255)
                                       tCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGca                                                                                                                <  2:2192817/150‑1 (MQ=255)
                                           ggCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGt                                                                                                                                     <  1:1160243/125‑1 (MQ=255)
                                           ggCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGt                                                                                                                                     >  2:1160243/1‑125 (MQ=255)
                                             ccAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTAcgccgcgcc                                                                                                                           >  1:375814/1‑133 (MQ=255)
                                             ccAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTAcgccgcgcc                                                                                                                           <  2:375814/133‑1 (MQ=255)
                                             ccAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCAccacca                                                                                                          >  2:1083459/1‑150 (MQ=255)
                                               aGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCaccaac                                                                                                        <  2:583846/150‑1 (MQ=255)
                                                gCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGGGGTACGCCGCGCCGATCACCAGCACCaccaacc                                                                                                       >  1:2199478/1‑150 (MQ=255)
                                                 caccacCCAGATCAGCGCCAGCAGCCAGCCGGTGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGgcg                                                                                                                                                   <  1:159830/105‑1 (MQ=255)
                                                 caccacCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGgcg                                                                                                                                                   >  2:159830/1‑105 (MQ=255)
                                                    caccCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAgg                                                                                                                                      <  2:1877844/115‑1 (MQ=255)
                                                    caccCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAgg                                                                                                                                      >  1:1877844/1‑115 (MQ=255)
                                                    caccCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCacaa                                                                                                   <  1:567490/150‑1 (MQ=255)
                                                     accCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCacaac                                                                                                  >  1:2054940/1‑150 (MQ=255)
                                                     accCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCacaac                                                                                                  >  1:1775991/1‑150 (MQ=255)
                                                      cccAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCacaacc                                                                                                 >  2:425487/1‑149 (MQ=255)
                                                       ccAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAg                                                                                                >  2:1710419/1‑150 (MQ=255)
                                                       ccAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAg                                                                                                >  1:1153729/1‑150 (MQ=255)
                                                       ccAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAg                                                                                                >  1:1646002/1‑150 (MQ=255)
                                                        cAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGCGAATCTGCAAGAGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGc                                                                                               >  2:192822/1‑150 (MQ=255)
                                                            tCAGCGCCAGCAGCCAGCCGGGGGCCCAGGTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACa                                                                                           >  1:1244490/1‑150 (MQ=255)
                                                               gcgcCGGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGt                                                                                                                                          >  2:336497/1‑100 (MQ=255)
                                                               gcgcCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATaa                                                                                        <  2:1186333/150‑1 (MQ=255)
                                                                   cagcagCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCTGGCGccct                                                                                                                                               >  2:1273147/1‑89 (MQ=255)
                                                                   cagcagCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCTGGCGccct                                                                                                                                               <  1:1273147/91‑3 (MQ=255)
                                                                   cagcagCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCgg                                                                                    >  1:444053/1‑150 (MQ=255)
                                                                          cagccGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTcctcct                                                                             >  2:1477055/1‑150 (MQ=255)
                                                                            gccGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTg                                                                           <  2:2112848/150‑1 (MQ=255)
                                                                                ggggCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAAc                                                                       >  2:554507/1‑150 (MQ=255)
                                                                                ggggCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAAc                                                                       <  1:160727/150‑1 (MQ=255)
                                                                                     aCTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAgg                                                                  <  2:444053/150‑1 (MQ=255)
                                                                                     aCTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAgg                                                                  >  1:654655/1‑150 (MQ=255)
                                                                                     aCTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAgg                                                                  <  2:17077/150‑1 (MQ=255)
                                                                                        cctgAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCt                                                               <  1:108529/149‑1 (MQ=255)
                                                                                         ctgAAACTGCCCATCAAAAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCtt                                                              >  2:629194/1‑150 (MQ=255)
                                                                                         ctgAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAg                                                                                                >  1:1469025/1‑116 (MQ=255)
                                                                                         ctgAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAg                                                                                                <  2:1469025/116‑1 (MQ=255)
                                                                                           taaaaTGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGa                                                            <  1:191537/149‑1 (MQ=255)
                                                                                                  cccATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGcac                                                                                                               >  2:1645911/1‑92 (MQ=255)
                                                                                                  cccATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGcac                                                                                                               <  1:1645911/92‑1 (MQ=255)
                                                                                                  cccATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTc                                                     <  1:425487/150‑1 (MQ=255)
                                                                                                  cccATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTc                                                     <  2:2054940/150‑1 (MQ=255)
                                                                                                    cATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCacaaca                                                                                                 <  1:34564/104‑1 (MQ=255)
                                                                                                     aTCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTAACCCATTCCAg                                                  >  1:492602/1‑150 (MQ=255)
                                                                                                      tCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGTGTGAGGTACGCCGCGCCGATCACCAGCACCACCAGGAAAAACACCTACATAATCGGGCCACCTTGGAACCAAGGGCTTGAGCCATTCCAGa                                                 >  2:1243656/1‑150 (MQ=255)
                                                                                                        aaCAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACa                                               <  1:1493656/150‑1 (MQ=255)
                                                                                                          cAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCtt                                                                            <  1:1417184/119‑1 (MQ=255)
                                                                                                           aGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCAc                                                                                                                     >  1:1981617/1‑77 (MQ=255)
                                                                                                           aGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCAc                                                                                                                     <  2:1981617/77‑1 (MQ=255)
                                                                                                           aGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCt                                            <  1:1710419/150‑1 (MQ=255)
                                                                                                           aGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCt                                            <  1:1083459/150‑1 (MQ=255)
                                                                                                              acaTAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGATCCTCCTTGAAACAAAGGGCTTGAGCCATTCCAGACAGCTTAg                                         >  1:1978397/1‑150 (MQ=255)
                                                                                                               caTAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGa                                        >  2:74194/1‑150 (MQ=255)
                                                                                                                aTAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGAc                                       >  1:795936/1‑150 (MQ=255)
                                                                                                                aTAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGACCGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGAc                                       >  1:1496307/1‑150 (MQ=255)
                                                                                                                 tAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCgg                                                                                                                                               >  1:327301/1‑45 (MQ=255)
                                                                                                                 tAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCgg                                                                                                                                               <  2:327301/45‑1 (MQ=255)
                                                                                                                       gccACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGc                                >  1:2437349/1‑150 (MQ=255)
                                                                                                                        ccACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCt                               <  1:1477055/150‑1 (MQ=255)
                                                                                                                         cACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCtg                              <  2:293425/150‑1 (MQ=255)
                                                                                                                            gATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTAcgccgcgc                                                                                                                            <  2:245960/53‑1 (MQ=255)
                                                                                                                            gATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTAcgccgcgc                                                                                                                            >  1:245960/1‑53 (MQ=255)
                                                                                                                             aTGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCTGTGaa                          <  2:654655/150‑1 (MQ=255)
                                                                                                                              tGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACTGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCTGTGAAc                         <  1:1348017/150‑1 (MQ=255)
                                                                                                                                 cGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCTGTGAACGGc                      <  2:1646002/150‑1 (MQ=255)
                                                                                                                                  gagaATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGCGCCAGGCCAGACAGCTTAGACGCAATGCTGTGAACGGCa                     >  2:2017058/1‑150 (MQ=255)
                                                                                                                                  gagaATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCTGTGAACGGCa                     >  1:439484/1‑150 (MQ=255)
                                                                                                                                  gagaATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCTGTGAACGGCa                     >  2:819355/1‑150 (MQ=255)
                                                                                                                                   agaATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCTGTGAACGGCAt                    >  1:1918691/1‑150 (MQ=255)
                                                                                                                                         tGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCt                               >  2:1790254/1‑133 (MQ=255)
                                                                                                                                         tGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCt                               <  1:1790254/133‑1 (MQ=255)
                                                                                                                                          gccaGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCTGTGAACGGCATGTTCAGa             >  2:387475/4‑150 (MQ=255)
                                                                                                                                           cAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCTGTGAACGGCATGTTCAGaa            <  2:1918691/150‑1 (MQ=255)
                                                                                                                                            accggCGCCTGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACa                                               >  1:778265/4‑114 (MQ=255)
                                                                                                                                            accggCGCCTGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACa                                               <  2:778265/111‑1 (MQ=255)
                                                                                                                                              gggCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAg                                                           <  1:312633/100‑1 (MQ=255)
                                                                                                                                              gggCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAg                                                           >  2:312633/1‑100 (MQ=255)
                                                                                                                                                gcgcCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCTGTGAACGGTATGTTCAGAACAGTg       >  2:153772/1‑150 (MQ=255)
                                                                                                                                                  gcCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCTGTGAACGGCATGTTCAGAACAGTGAc     >  2:1265733/1‑150 (MQ=255)
                                                                                                                                                    cAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCTGTGAACGGCATGTTCAGAACAGTGACAg   >  2:464612/1‑150 (MQ=255)
                                                                                                                                                     aGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCTGTGAACGGCATGTTCAGAACAGTGACAgg  <  1:1265733/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTTGCGTCGCCACTGCGGGAGCGCAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCTGTGAACGGCATGTTCAGAACAGTGACAGG  >  NC_002947/2136583‑2136881

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