Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 191,912 Δ1 bp 7.2% coding (1606/2157 nt) paxB ← toxin secretion ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947191,9120C. 7.2% 214.9 / 11.5 69coding (1606/2157 nt)paxBtoxin secretion ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (27/37);  new base . (3/2);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01

GTAGCGGGCGCCGCTGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGACC  >  NC_002947/191763‑192054
                                                                                                                                                     |                                                                                                                                               
gTAGCGGGCGCCGCTGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTc                                                                                                                                                 <  2:1418190/150‑1 (MQ=255)
    cGGGCGCCGCTGACCAGGTTGTCACGCACGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTAc                                                                                                                                             >  2:713277/1‑150 (MQ=255)
       gcgcCGCTGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAgg                                                                                                                                         <  1:97726/150‑1 (MQ=255)
         gccgcTGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGcc                                                                                                                                       <  1:852010/150‑1 (MQ=255)
              tGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCcgacga                                                                                                                                  <  1:888531/150‑1 (MQ=255)
              tGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCcgacga                                                                                                                                  >  2:761738/1‑150 (MQ=255)
              tGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCcgacga                                                                                                                                  <  1:926427/150‑1 (MQ=255)
               gACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGAt                                                                                                                                 <  1:801281/150‑1 (MQ=255)
                   aGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGc                                                                                                                             >  2:118917/1‑150 (MQ=255)
                         tCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGcc                                                                                                                       >  2:1164459/1‑150 (MQ=255)
                               aaGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCg                                                                                                                                                >  1:928976/1‑120 (MQ=255)
                               aaGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCg                                                                                                                                                <  2:928976/120‑1 (MQ=255)
                                    gCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCtt                                                                                                            <  1:905632/150‑1 (MQ=255)
                                     ccGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCATCAGCCAGCCACCATCGGCCGCTTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTg                                                                                                           <  1:175288/150‑1 (MQ=255)
                                      cGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCACCGGCCTGTACGAGGCcgatgg                                                                                                                                  >  1:377449/1‑123 (MQ=255)
                                      cGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCACCGGCCTGTACGAGGCcgatgg                                                                                                                                  <  2:377449/126‑4 (MQ=255)
                                        gCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCGCCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCg                                                                                                        >  1:1008592/1‑150 (MQ=255)
                                          cagcagCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGa                                                                                                      <  1:594327/150‑1 (MQ=255)
                                               gcTGGATCTCCTGGGGCACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGc                                                                                                 >  2:1128998/1‑150 (MQ=255)
                                                  ggATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGc                                                                                              <  2:953895/150‑1 (MQ=255)
                                                  ggATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGc                                                                                              >  1:148495/1‑150 (MQ=255)
                                                     tatCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGc                                                                                           <  2:1026200/150‑1 (MQ=255)
                                                       tCCTGGGGGACGTAGCCGATGTTGTGGCGCTGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCAGAGCCGCTGCGGCCg                                                                                         >  1:1328450/1‑150 (MQ=255)
                                                       tCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGCACAGCCCGAGCATGAGTTTGTCCAGCGTGCCGGTTCCGGAGCCGCTGAGGACg                                                                                         >  2:727570/1‑150 (MQ=255)
                                                       tCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATAGCTCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCg                                                                                         >  2:1267674/1‑150 (MQ=255)
                                                         cTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCgat                                                                                       >  1:765916/1‑150 (MQ=255)
                                                          tGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCgatg                                                                                      >  2:1319032/1‑150 (MQ=255)
                                                            ggggACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCgatgat                                                                                    >  2:1393564/1‑150 (MQ=255)
                                                                  gTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCga                                                                                        >  1:6097/1‑140 (MQ=255)
                                                                  gTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCga                                                                                        <  2:6097/140‑1 (MQ=255)
                                                                     gCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCtt                                                                           >  2:445929/1‑150 (MQ=255)
                                                                             tgtgGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTtctc                                                                        >  2:694114/1‑145 (MQ=255)
                                                                             tgtgGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTtctc                                                                        <  1:694114/145‑1 (MQ=255)
                                                                             tgtgGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAg                                                                   <  2:289163/150‑1 (MQ=255)
                                                                              gtgGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGGTCAATCTGCTGCGGCCGATGATGCCGACCTTCTCGCCAgg                                                                  >  1:784985/1‑150 (MQ=255)
                                                                               tgGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAggg                                                                 <  2:616322/150‑1 (MQ=255)
                                                                                   gcAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGa                                                             <  1:762311/150‑1 (MQ=255)
                                                                                   gcAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGa                                                             <  1:761738/150‑1 (MQ=255)
                                                                                            tGACATCGATCTGGCGACTGTCGACGACACAACCCCGCAGCGAGCCACCCTCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGt                                                    <  2:784985/150‑1 (MQ=255)
                                                                                             gACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCGCCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGtt                                                   <  2:1008592/150‑1 (MQ=255)
                                                                                             gACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGtt                                                   <  1:1036231/150‑1 (MQ=255)
                                                                                             gACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGtt                                                   <  1:255553/150‑1 (MQ=255)
                                                                                                aTCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCa                                                                    >  2:1221744/1‑130 (MQ=255)
                                                                                                aTCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCa                                                                    <  1:1221744/130‑1 (MQ=255)
                                                                                                 tCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGata                                               <  1:1319032/150‑1 (MQ=255)
                                                                                                      cTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATAttctt                                          <  1:1055937/150‑1 (MQ=255)
                                                                                                       tGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATAttcttc                                         <  2:814395/150‑1 (MQ=255)
                                                                                                           gAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAggg                                     >  1:685796/1‑150 (MQ=255)
                                                                                                                  gACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGt                                                    <  1:276971/128‑1 (MQ=255)
                                                                                                                  gACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGt                                                    >  2:276971/1‑128 (MQ=255)
                                                                                                                   aCGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTc                             <  1:1164459/150‑1 (MQ=255)
                                                                                                                      ccATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTctgc                          <  1:445929/150‑1 (MQ=255)
                                                                                                                       cATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTctgct                         >  1:1268339/1‑150 (MQ=255)
                                                                                                                        aTCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTctgctg                        >  1:137084/1‑150 (MQ=255)
                                                                                                                         tCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTctgc                          <  2:1044655/147‑1 (MQ=255)
                                                                                                                         tCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTctgc                          >  1:1044655/1‑147 (MQ=255)
                                                                                                                            aCCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTc                    <  2:148495/150‑1 (MQ=255)
                                                                                                                            aCCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCACTGCGGCCGATGATGCCGACCTTCTCGGCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTc                    <  1:713277/150‑1 (MQ=255)
                                                                                                                             ccAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCAGATGATGCCGACCTTCTCGCCAGGGCGGCTGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCg                   <  1:1393564/150‑1 (MQ=255)
                                                                                                                               agcagcGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGa                 <  2:685796/150‑1 (MQ=255)
                                                                                                                               agcagcGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGa                 >  1:64424/1‑150 (MQ=255)
                                                                                                                               agcagcGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTATGCTGGttgtga                 <  2:765916/150‑5 (MQ=255)
                                                                                                                                       gCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGAAATACGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGaa         <  1:727570/150‑1 (MQ=255)
                                                                                                                                          atcgtCGGCCTGTACGAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTc      >  1:570800/5‑150 (MQ=255)
                                                                                                                                             aTCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGAc   >  1:718580/1‑150 (MQ=255)
                                                                                                                                              tCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGTTGGGGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGAcc  <  1:155189/150‑1 (MQ=255)
                                                                                                                                              tCGGCCTGTACGAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGa    >  2:237712/1‑148 (MQ=255)
                                                                                                                                              tCGGCCTGTACGAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGa    <  1:237712/148‑1 (MQ=255)
                                                                                                                                                    tCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATaa              >  1:713236/1‑132 (MQ=255)
                                                                                                                                                    tCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATaa              <  2:713236/132‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GTAGCGGGCGCCGCTGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGACC  >  NC_002947/191763‑192054

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: