Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,726,562 (T)5→4 coding (971/2574 nt) clpB ← ATP‑dependent chaperone ClpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,726,5580T.98.9% 414.8 / ‑1.8 94coding (975/2574 nt)clpBATP‑dependent chaperone ClpB
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (46/47);  total (47/47)
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.92e-01

Only 100 of 107 total aligned reads displayed.

TCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGT  >  NZ_CP009273/2726469‑2726643
                                                                                         |                                                                                     
tCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCt                                                                                       >  2:2782372/1‑90 (MQ=255)
 cTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCtt                                                                                      <  2:1829467/90‑1 (MQ=255)
 cTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCtt                                                                                      >  2:2874891/1‑90 (MQ=255)
      aGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCttttcaa                                                                                 <  2:1051116/90‑4 (MQ=255)
      aGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCttttcaa                                                                                 <  2:1458085/90‑4 (MQ=255)
      aGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCttttcaa                                                                                 <  2:2052145/90‑4 (MQ=255)
      aGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCttttcaa                                                                                 <  2:2799586/90‑4 (MQ=255)
       ggCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCttttcaat                                                                                >  1:2982293/1‑86 (MQ=255)
       ggCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCttttcaat                                                                                >  2:586399/1‑86 (MQ=255)
         ccACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGt                                                                             <  2:487619/90‑1 (MQ=255)
           aCGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTAc                                                                           >  2:505353/1‑90 (MQ=255)
           aCGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTAc                                                                           >  1:1409984/1‑90 (MQ=255)
           aCGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTAc                                                                           >  1:1010263/1‑90 (MQ=255)
           aCGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTAc                                                                           >  1:1217567/1‑90 (MQ=255)
           aCGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTAc                                                                           >  2:2506168/1‑90 (MQ=255)
           aCGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTAc                                                                           >  2:2605653/1‑90 (MQ=255)
           aCGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTAc                                                                           >  1:643042/1‑90 (MQ=255)
               aGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGc                                                                       <  1:505353/90‑1 (MQ=255)
                 aaTCGCACTGGCATCTTCAACAGAAGGCTCGACAACAAACCCTTTCTGGAAACGACCTTCCAGCGCAGCATCTTTTCAATGTACTGGCGa                                                                     >  1:2255940/1‑90 (MQ=255)
                  aTCGCAATGGTATCTTCAGCAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCgct                                                                    >  1:2063369/1‑88 (MQ=255)
                    cGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGatat                                                                  <  1:649210/90‑1 (MQ=255)
                    cGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGatat                                                                  <  2:160672/90‑1 (MQ=255)
                     gCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATAtt                                                                 >  1:40283/1‑90 (MQ=255)
                          ggTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCttttcaa                                                                                 >  2:2713792/1‑67 (MQ=255)
                          ggTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCttttcaa                                                                                 <  1:2713792/70‑4 (MQ=255)
                          ggTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTAATGTACTGGCGATATTCGTCa                                                            >  1:1211044/1‑90 (MQ=255)
                          ggTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCa                                                            >  1:2025543/1‑90 (MQ=255)
                          ggTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCa                                                            >  2:1346253/1‑90 (MQ=255)
                          ggTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCa                                                            >  2:3012472/1‑90 (MQ=255)
                            tATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAg                                                          >  2:1438125/1‑90 (MQ=255)
                             aTCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGc                                                         <  1:1154381/90‑1 (MQ=255)
                             aTCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGc                                                         <  2:1503693/90‑1 (MQ=255)
                             aTCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGc                                                         >  2:168508/1‑90 (MQ=255)
                             aTCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGc                                                         >  2:1930867/1‑90 (MQ=255)
                             aTCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGc                                                         >  2:2425388/1‑90 (MQ=255)
                              tCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTAc                                                                           >  1:49332/1‑71 (MQ=255)
                              tCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTAc                                                                           <  2:49332/71‑1 (MQ=255)
                                ttCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGcgtc                                                      >  1:735909/1‑90 (MQ=255)
                                 tCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGcgtcg                                                     <  1:2874891/90‑1 (MQ=255)
                                  cAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCttttcaa                                                                                 <  2:1923156/62‑4 (MQ=255)
                                  cAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGcgtcgt                                                    >  2:1570840/1‑90 (MQ=255)
                                  cAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGcgtcgt                                                    >  1:1220274/1‑90 (MQ=255)
                                   aaCAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTg                                                   <  1:301039/90‑1 (MQ=255)
                                   aaCAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTg                                                   <  1:3012472/90‑1 (MQ=255)
                                   aaCAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTg                                                   <  1:375101/90‑1 (MQ=255)
                                   aaCAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTg                                                   <  1:1346253/90‑1 (MQ=255)
                                   aaCAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTACCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTg                                                   <  2:78260/90‑1 (MQ=255)
                                    aCAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTgg                                                  >  1:1805572/1‑90 (MQ=255)
                                    aCAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTgg                                                  >  1:2397771/1‑90 (MQ=255)
                                    aCAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTgg                                                  >  1:2039511/1‑90 (MQ=255)
                                     cAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGc                                                 <  1:1930867/90‑1 (MQ=255)
                                     cAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGc                                                 <  1:1570840/90‑1 (MQ=255)
                                     cAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGc                                                 <  1:1941162/90‑1 (MQ=255)
                                     cAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGc                                                 <  1:2425388/90‑1 (MQ=255)
                                        aaGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCAcc                                              <  2:332646/90‑1 (MQ=255)
                                        aaGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCAcc                                              <  1:2115279/90‑1 (MQ=255)
                                        aaGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCAcc                                              <  2:531917/90‑1 (MQ=255)
                                        aaGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCAcc                                              <  2:1190445/90‑1 (MQ=255)
                                         aGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCt                                             >  1:1411891/1‑90 (MQ=255)
                                         aGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCt                                             >  2:2689397/1‑90 (MQ=255)
                                         aGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCt                                             >  2:77981/1‑90 (MQ=255)
                                           gCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTAc                                           <  2:1411891/90‑1 (MQ=255)
                                             tCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGc                                         <  2:619962/90‑1 (MQ=255)
                                             tCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGc                                         <  1:1713874/90‑1 (MQ=255)
                                             tCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGc                                         <  2:2717433/90‑1 (MQ=255)
                                             tCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGc                                         <  1:168508/90‑1 (MQ=255)
                                             tCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGc                                         <  2:105561/90‑1 (MQ=255)
                                             tCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGc                                         <  2:324171/90‑1 (MQ=255)
                                             tCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGc                                         <  1:2177082/90‑1 (MQ=255)
                                             tCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGc                                         <  2:1820492/90‑1 (MQ=255)
                                              cGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGc                                                         >  2:2997128/1‑73 (MQ=255)
                                              cGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGc                                                         <  1:2997128/73‑1 (MQ=255)
                                              cGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGc                                         <  1:357796/89‑1 (MQ=255)
                                              cGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGc                                         >  2:357796/1‑89 (MQ=255)
                                              cGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCa                                        >  2:2667418/1‑90 (MQ=255)
                                               ggCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAg                                       <  1:2305812/90‑1 (MQ=255)
                                                     aaacaCTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAAt                                 >  2:2224791/1‑90 (MQ=255)
                                                       acacTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTc                               <  1:661547/90‑1 (MQ=255)
                                                           tttCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCAcc                                              <  1:2110024/71‑1 (MQ=255)
                                                           tttCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCAcc                                              >  2:2110024/1‑71 (MQ=255)
                                                           tttCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTcacca                           >  1:2942068/1‑90 (MQ=255)
                                                            ttCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTcaccac                          <  1:1378747/90‑1 (MQ=255)
                                                            ttCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTcaccac                          <  2:2838226/90‑1 (MQ=255)
                                                            ttCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTcaccac                          >  1:1344282/1‑90 (MQ=255)
                                                            ttCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTcaccac                          <  1:270809/90‑1 (MQ=255)
                                                            ttCTGGAAACGACGTTCCAGCGCAGCATCTTTCCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTcaccac                          >  2:2199701/1‑90 (MQ=255)
                                                              cTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGc                                         >  2:2469297/1‑73 (MQ=255)
                                                              cTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGc                                         <  1:2469297/73‑1 (MQ=255)
                                                              cTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCAcgc                        >  2:2428903/1‑90 (MQ=255)
                                                                     cgacgTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAgcg                 <  2:2374749/90‑1 (MQ=255)
                                                                        cgTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGcgtcg                                                     >  2:2146175/1‑51 (MQ=255)
                                                                        cgTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGcgtcg                                                     <  1:2146175/51‑1 (MQ=255)
                                                                             cAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTc         <  2:2544289/90‑1 (MQ=255)
                                                                               gcgcAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAg       >  2:436301/1‑90 (MQ=255)
                                                                                 gcagcaTCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCa     >  2:2006842/1‑90 (MQ=255)
                                                                                   agcaTCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATg   <  1:2689397/90‑1 (MQ=255)
                                                                                    gcaTCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGc      >  1:1229203/1‑86 (MQ=255)
                                                                                    gcaTCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGc      <  2:1229203/86‑1 (MQ=255)
                                                                                    gcaTCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGt  >  2:2534408/1‑90 (MQ=255)
                                                                                    gcaTCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGt  >  2:2472280/1‑90 (MQ=255)
                                                                                         |                                                                                     
TCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGT  >  NZ_CP009273/2726469‑2726643

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