Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 609,782 A→C 5.3% L356R (CTG→CGG)  PP_0525 ← B12 family TonB‑dependent receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947609,7820AC 5.3% 220.0 / 4.6 76L356R (CTG→CGG) PP_0525B12 family TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base A (39/33);  new base C (2/2);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

GATTCTTGTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTCCACTGCTCG  >  NC_002947/609633‑609931
                                                                                                                                                     |                                                                                                                                                     
gATTCTTGTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAAcagca                                                                                                                                                       <  2:770895/150‑1 (MQ=255)
gATTCTTGTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAAcagca                                                                                                                                                       <  2:1368703/150‑1 (MQ=255)
  ttCTTGTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAgg                                                                                                                                                     <  1:748541/150‑1 (MQ=255)
   tCTTGTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGc                                                                                                                                                    >  1:719298/1‑150 (MQ=255)
       gTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCtgtt                                                                                                                                                <  2:800783/150‑1 (MQ=255)
           tGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTgtgt                                                                                                                                            <  1:985237/150‑1 (MQ=255)
                   gCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCCCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAg                                                                                                                                    >  2:1159550/1‑150 (MQ=255)
                     cAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAggg                                                                                                                                  >  2:653699/1‑150 (MQ=255)
                         cTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAg                                                                                                                              <  2:1484530/150‑1 (MQ=255)
                         cTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAg                                                                                                                              >  1:1484530/1‑150 (MQ=255)
                             gagtgaaaGTATTCGTCATGAAAGCGTTGCTGGACGAAAAGCGCATGAGTCAAGCGTCTGTATTCCGTAAAATCCTCACGCACATGGACACGGTCTTCTTACCAGTCGCTCCCCAACAGCAGGCTGTTATTTTAACCCAGGGCTAGgttg                                                                                                                          >  2:34700/6‑150 (MQ=255)
                               gtAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGgttgtt                                                                                                                        >  2:1587345/1‑150 (MQ=255)
                                   aaGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGAGGCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAAGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGc                                                                                                                    >  2:1574043/1‑150 (MQ=255)
                                         tCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGtt                                                                                                              <  1:1529695/150‑1 (MQ=255)
                                           gTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGc                                                                                                            >  2:1454159/1‑150 (MQ=255)
                                           gTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGc                                                                                                            <  1:167718/150‑1 (MQ=255)
                                            tCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCCTGCtgttgggt                                                                                                                                            >  2:389782/1‑114 (MQ=255)
                                            tCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCCTGCTGTTgggt                                                                                                                                            <  1:389782/117‑1 (MQ=255)
                                             cATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCCCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCtt                                                                                                          <  1:1159550/150‑1 (MQ=255)
                                                 aaaGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCt                                                                                                                      <  2:1307827/134‑1 (MQ=255)
                                                 aaaGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCt                                                                                                                      >  1:1307827/1‑134 (MQ=255)
                                                     cGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGc                                                                                                                                 <  1:1751565/119‑1 (MQ=255)
                                                     cGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGc                                                                                                                                 >  2:1751565/1‑119 (MQ=255)
                                                      gATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCggt                                                                                                 >  1:988533/1‑150 (MQ=255)
                                                        tGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCCTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCggtag                                                                                               <  1:215828/150‑1 (MQ=255)
                                                            ggACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTc                                                                                                      <  2:1360416/139‑1 (MQ=255)
                                                            ggACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTc                                                                                                      >  1:1360416/1‑139 (MQ=255)
                                                             gACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTcgcg                                                                                                   <  2:1097633/141‑1 (MQ=255)
                                                             gACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTcgcg                                                                                                   >  1:1097633/1‑141 (MQ=255)
                                                                 aaCAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGaaa                                                                                      >  1:1564031/1‑150 (MQ=255)
                                                                  aCAGCGCATGATTCCAGCTGATGTTTTCCGTGAACCCCGCACTCGAATGGACACGGTCTCCGTAACTGAGTTTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAg                                                                                     <  1:34700/150‑1 (MQ=255)
                                                                    aGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGt                                                                                   <  2:1564031/150‑1 (MQ=255)
                                                                       gcATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTgcg                                                                                >  2:78752/1‑150 (MQ=255)
                                                                             ttCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAg                                                                                     >  1:1017797/1‑139 (MQ=255)
                                                                             ttCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAg                                                                                     <  2:1017797/139‑1 (MQ=255)
                                                                             ttCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACATCCACCCGAGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAg                                                                          >  1:19992/1‑150 (MQ=255)
                                                                              tCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGt                                                                         <  1:1587345/150‑1 (MQ=255)
                                                                               ccAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTg                                                                        <  1:116788/150‑1 (MQ=255)
                                                                                      cTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACGAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGc                                                                 >  2:377855/1‑150 (MQ=255)
                                                                                        gTCTTCCGTGAAGTCCTCACTCGCATGGACATGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCtt                                                               <  2:322928/150‑1 (MQ=255)
                                                                                        gTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCtt                                                               >  1:625288/1‑150 (MQ=255)
                                                                                        gTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCtt                                                               <  2:625288/150‑1 (MQ=255)
                                                                                        gTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCtt                                                               <  2:719298/150‑1 (MQ=255)
                                                                                           ttCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGtc                                                            >  2:522023/1‑150 (MQ=255)
                                                                                           ttCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGtc                                                            <  1:78752/150‑1 (MQ=255)
                                                                                            tCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGtct                                                           >  1:834595/1‑150 (MQ=255)
                                                                                              cGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAgg                                                                                                                             >  2:633229/1‑82 (MQ=255)
                                                                                              cGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAgg                                                                                                                             <  1:633229/82‑1 (MQ=255)
                                                                                              cGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTAt                                                                                           <  1:1327988/116‑1 (MQ=255)
                                                                                              cGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTAt                                                                                           >  2:1327988/1‑116 (MQ=255)
                                                                                                       cTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGa                                                                   <  2:1243401/131‑1 (MQ=255)
                                                                                                       cTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGa                                                                   >  1:1243401/1‑131 (MQ=255)
                                                                                                        tCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCAGTTTTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTGTTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATcg                                               >  1:49505/1‑150 (MQ=255)
                                                                                                            tCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGt                                           >  1:327040/1‑150 (MQ=255)
                                                                                                              gCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGttt                                         >  1:30868/1‑150 (MQ=255)
                                                                                                              gCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGttt                                         >  2:643383/1‑150 (MQ=255)
                                                                                                              gCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGttt                                         >  2:453691/1‑150 (MQ=255)
                                                                                                                   gACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGGTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCt                                    >  1:1571112/1‑150 (MQ=255)
                                                                                                                      acGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATcgcg                                             >  1:474392/1‑138 (MQ=255)
                                                                                                                      acGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATcgcg                                             <  2:474392/138‑1 (MQ=255)
                                                                                                                       cGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCagttt                                <  1:67698/150‑3 (MQ=255)
                                                                                                                        ggTGTTCTTACCAGTCGCTGCCCAAGAGCAGGCTGTTGTGTTCACCCCGGGCTATGTTGTTCTGCCAAGTGGCTTGGTCGCGGTAGGTATTGAAAGGTGCGGTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTGGCtgtgtg                               <  2:39465/150‑3 (MQ=255)
                                                                                                                        ggTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTggg                               >  1:1141029/1‑150 (MQ=255)
                                                                                                                        ggTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTggg                               >  1:1875182/1‑150 (MQ=255)
                                                                                                                            ttCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAg                           >  1:1336360/1‑150 (MQ=255)
                                                                                                                               gTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGt                                                                         <  2:1518071/101‑1 (MQ=255)
                                                                                                                               gTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGt                                                                         >  1:1518071/1‑101 (MQ=255)
                                                                                                                                tACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTACTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGGTTTCGCTGTGGGCCAGCTCa                       >  1:825080/1‑150 (MQ=255)
                                                                                                                                tACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGATAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGTTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGGCAGCTCa                       >  2:128566/1‑150 (MQ=255)
                                                                                                                                  ccAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCaaa                     <  2:834595/150‑1 (MQ=255)
                                                                                                                                   cAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATcgc                                              <  1:1681447/124‑1 (MQ=255)
                                                                                                                                   cAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATcgc                                              >  2:1681447/1‑124 (MQ=255)
                                                                                                                                   cAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAAc                    >  1:1840722/1‑150 (MQ=255)
                                                                                                                                   cAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAAc                    <  1:643383/150‑1 (MQ=255)
                                                                                                                                      tCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGa                 <  2:327040/150‑1 (MQ=255)
                                                                                                                                            cacAACAGCCTGCTGTTGGGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCtgt                                  <  1:1694602/125‑1 (MQ=255)
                                                                                                                                            cacAACAGCCTGCTGTTGGGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCtgt                                  >  2:1694602/3‑127 (MQ=255)
                                                                                                                                               aaCAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTCCAc        <  1:1454159/150‑1 (MQ=255)
                                                                                                                                                     agGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTCCACTGCTCg  <  1:653699/150‑1 (MQ=255)
                                                                                                                                                     agGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTCCACTGCTCg  >  2:1091800/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATTCTTGTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTCCACTGCTCG  >  NC_002947/609633‑609931

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