Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,050,741 C→T 6.1% E518K (GAG→AAG)  pcrA ← ATP‑dependent DNA helicase PcrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,050,7410CT 6.1% 190.4 / 2.8 65E518K (GAG→AAG) pcrAATP‑dependent DNA helicase PcrA
Reads supporting (aligned to +/- strand):  ref base C (42/19);  new base T (2/2);  total (44/21)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GGAAATTCAAGGCCCTTAGCCGAGTGCATCGTCATTAGTGTTACGCCATTTTCAGTATCTGCCTCATCAATATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGATAGCCAGACTTTTGTAATACTTCATCAACTATTTCATGAATTTCTAAAAA  >  CP000730/2050602‑2050879
                                                                                                                                           |                                                                                                                                          
ggAAATTCAAGGCCCTTAGCCGAGTGCATCGTCATTAGTGTTAAGCCATTTTCAGTATCTGCCTCATCAATATCAGCTACTACCGATAAATCCGTTAAAAAAGAAATTAATGACTGTTATCCTAGTGGGGTGGTCTCCTCa                                                                                                                                           >  2:122032/1‑141 (MQ=255)
 taaaCTCCAGGCCCTTCGGCGAGATCATCGTCATTAGTGTTATGCCATGTTCAGTATCTTCCTCATCAATATCAGCTACTAATGAAAAAAACATAAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCAt                                                                                                                                          <  2:368005/140‑1 (MQ=255)
 gAAATTCAAGGCCCTTAGCCGAGTGCATCGTCATTAGTGTTACGCCATTTTCAGTATCTGCCTCATCAATATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCAt                                                                                                                                          >  1:28792/1‑141 (MQ=255)
  aaaTTCAAGGCCCTTAGCCGAGTGCATCGTCATTAGTGTTACGCCATTTTCAGTATCTGCCTCATCAATATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTTATGGGGTATTTTCCTCATa                                                                                                                                         >  1:183858/1‑141 (MQ=255)
  aaaTTCAAGGCCCTTAGCCGAGTGCATCGTCATTAGTGTTACGCCATTTTCAGTATCTGCCTCATCAATATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATa                                                                                                                                         >  1:287803/1‑141 (MQ=255)
  aaaTTCAAGGCCCTTAGCCGAGTGCATCGTCATTAGTGTTACGCCATTTTCAGTATCTGCCTCATCAATATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATa                                                                                                                                         >  2:547925/1‑141 (MQ=255)
  aaaTTCAAGGCCCTTAGCCGAGTGCATCGTCATTAGTGTTACGCCATTTTCAGTATCTGCCTCATCAATATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATa                                                                                                                                         >  1:132103/1‑141 (MQ=255)
  aaaTTCAAGGCCCTTAGCCGAGTGCATCGTCATTAGTGTTACGCCATTTTCAGTATCTGCCTCATCAATATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATa                                                                                                                                         >  1:571964/1‑141 (MQ=255)
     ttaaaGGCCCTTACCCGAGTGTATCGTCATTCGTGTTCCGCCATTTTCAGTATCTGCCTCATCAATATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTc                                                                                                                                      <  1:315606/138‑1 (MQ=255)
              cTTAGCCGAGTGCATCGTCATTAGTGTTACGCCATTTTCAGTATCTGCCTCATCAATATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTAc                                                                                                                             >  2:561898/1‑141 (MQ=255)
                   ccGAGTGCATCGTCATTAGTGTTACGCCATTTTCAGTATCTGCCTCATCAATATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGCCTCTTCGTGTAATGGGGGATTTTCCTCATAGTCTTTTGGTACTGACa                                                                                                                        >  1:361827/1‑141 (MQ=255)
                       gTGCATCGTCATTAGTGTTACGCCATTTTCAGTATCTGCCTCATCAATATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTCGGTACTGACATaaa                                                                                                                    <  1:363128/141‑1 (MQ=255)
                          cATCGTCATTAGTGTTACGCCATTTTCAGTATCTGCCTCATCAATATCAGCTACTAACGATAAATCCGTTTAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGCCTTTTGGTACTGACATAGATTc                                                                                                                 >  2:197637/1‑141 (MQ=255)
                          cATCGTCATTAGTGTTACGCCATTTTCAGTATCTGCCTCATCAATATCAGCTACTAACAATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCGTAGTCTTTTGGTACTGACATAAATTc                                                                                                                 >  1:424956/1‑141 (MQ=255)
                          cATCGTCATTAGTGATACGCCATTTTCAGTATCTGCCTCATCAAGATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTc                                                                                                                 >  1:379928/1‑141 (MQ=255)
                           aTCGTCATTAGTGTTACGCCATTTTCAGTATCTGCCTCATCAATATCAGCTGCTAACGATAAATCCGTTAAAAAGTTAATTAGTGACTGTTCTTTGATTGTGGTATTTTCCTCATAGTCTTTTGGTACTGCAATAAATTCa                                                                                                                >  1:397566/1‑141 (MQ=255)
                                 aTTAGTGTTACGCCATTTTCAGTATCTGCCTCATCAATATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGAtt                                                                                                          >  1:641842/1‑140 (MQ=255)
                                 aTTAGTGTTACGCCATTTTCAGTATCTGCCTCATCAATATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATg                                                                                                          >  1:559188/1‑141 (MQ=255)
                                           cGCCATTTTCAGTATCTGCCTCATCAATATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTTATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATc                                                                                                >  1:666324/1‑141 (MQ=255)
                                           cGCCATTTTCAGTATCTGCCTCATCAATATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTTATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATc                                                                                                >  1:666540/1‑141 (MQ=255)
                                            gCCATTTTCAGTATCTGCCTCATCAATATCTGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCt                                                                                               <  2:443741/141‑1 (MQ=255)
                                               aTTTTCAGTATCTGCCTCATCAATATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACt                                                                                            >  1:376855/1‑141 (MQ=255)
                                                    cAGTATCTGCCTCATCAATATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCgag                                                                                       >  1:665311/1‑141 (MQ=255)
                                                       tATCTGCCTCATCAATATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGAtt                                                                                    >  1:399373/1‑141 (MQ=255)
                                                       tATCTGCCTCATCAATATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAAAtt                                                                                    >  2:546647/1‑141 (MQ=255)
                                                         tCTGCCTCATCAATATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCt                                                                                  <  2:106816/141‑1 (MQ=255)
                                                             ccTCATCAATATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATg                                                                              >  1:613039/1‑141 (MQ=255)
                                                              cTCATCAATATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGt                                                                             <  2:402533/141‑1 (MQ=255)
                                                                 atcaATATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTAtt                                                                          <  2:347560/141‑1 (MQ=255)
                                                                  tcaATAACAGCTAATAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATATTCCTCATAGTCCTTTGGTAATGTCCTAAATTCCTCGAGTATTTCTAATCTACTTCGCGCTTCTCatttattt                                                                         >  2:90111/1‑141 (MQ=255)
                                                                         cAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGCTTCTAATGTATTTTCACGtt                                                                  >  1:38216/1‑141 (MQ=255)
                                                                             tactaACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAg                                                              <  1:399238/141‑1 (MQ=255)
                                                                              actaACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGc                                                             >  1:539074/1‑141 (MQ=255)
                                                                              actaACGATAAATCCGTTAAAAAGTTAATTAATGACTGCTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGc                                                             >  1:446770/1‑141 (MQ=255)
                                                                                 aaCGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCAtt                                                          >  2:48350/1‑141 (MQ=255)
                                                                                     aTAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTtctc                                                      <  2:339375/141‑1 (MQ=255)
                                                                                       aaaTCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTTATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGa                                                    <  2:666324/141‑1 (MQ=255)
                                                                                       aaaTCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTTATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGa                                                    <  2:666540/141‑1 (MQ=255)
                                                                                        aaTCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGAt                                                   >  2:573577/1‑141 (MQ=255)
                                                                                         aTCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTTTCGATa                                                  >  1:619085/1‑141 (MQ=255)
                                                                                          tCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGATAg                                                 <  2:591701/141‑1 (MQ=255)
                                                                                           ccGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGATAGc                                                <  1:547925/141‑1 (MQ=255)
                                                                                           ccGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATATTACTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGGTTCTAATGTATTTTCACGTTCAAGCATTTCTCCATAGc                                                >  2:676777/1‑141 (MQ=255)
                                                                                            cGTTAAAAAGTTAATTAATGACTGTTCTTCTACTGGGGTATTTTTCTCATAGTCTCTTGGTACTGACATAACTTCATCGATGTTTTCTCATGTTTTTCGAGCTTCTAATGTATTATCACGTTCACCCATTTCTCGATCGcc                                               >  2:310249/1‑141 (MQ=255)
                                                                                             tttAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGATAGCCa                                              <  2:451856/140‑1 (MQ=255)
                                                                                                 aaaaGTTAATTAATGACTGGTCTTCTAATGGGGTATTTTCCTCATAGTGTTTTGGTACTGAAATCAAATCATAGATGTTTTCTGATCTACGTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGATAGCCAGACt                                          >  1:518380/1‑141 (MQ=255)
                                                                                                  aaaGTTCATTAATTACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGAAATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGATAGCCAGACtt                                         <  2:397566/141‑1 (MQ=255)
                                                                                                     gTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGATAGCCAGACttttt                                      >  2:544777/1‑140 (MQ=255)
                                                                                                       taattaatGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGATAGCCAGACTTTTGTa                                    >  2:314284/1‑141 (MQ=255)
                                                                                                          ttaatGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGATAGCCAGACTTTTGTAATa                                 <  1:594283/141‑1 (MQ=255)
                                                                                                              tGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGCTTCTAATGTATTTTCACGTTCAAGCATTTCTCGATGGCCAGGCTTTTGTATTACTtc                             >  2:157911/1‑141 (MQ=255)
                                                                                                                   tttctttTAATGGGGTATTTTCCTCATTGTCTTTTGGTACTGACATAAATTCATCGATGTTCTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGATAGCCAGACTTTTGTAATACTTCATCaa                        <  1:577236/140‑1 (MQ=255)
                                                                                                                   gTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGATAGCCAGACTTTTGTAATACTTCATCaa                        <  1:395696/141‑1 (MQ=255)
                                                                                                                         ctAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGATAGCCAGACTTTTGTAATACTTCATCAACTAttt                  <  2:183858/141‑1 (MQ=255)
                                                                                                                         ctAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGATAGCCAGACTTTTGTAATACTTCATCAACTAttt                  >  2:598823/1‑141 (MQ=255)
                                                                                                                                 gTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGATAGCCAGACTTTTGTAATACTTCATCAACTATTTCATGatct          >  2:645188/1‑138 (MQ=255)
                                                                                                                                   aTTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGATAGCCAGACTTTTGTAATACTTCATCAACTATTTCATGAATTTc        <  2:38216/141‑1 (MQ=255)
                                                                                                                                     tttCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGATAGCCAGACTTTTGTAATACTTCATCAACTATTTCATGAATTTCTa      >  1:183786/1‑141 (MQ=255)
                                                                                                                                      ttCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGATAGCCAGACTTTTGTAATACTTCATCAACTATTTCATGAATTTCTaa     >  2:387053/1‑141 (MQ=255)
                                                                                                                                      ttCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGATAGCCAGACTTTTGTAATACTTCATCAACTATTTCATGAATTTCTaa     >  1:69100/1‑141 (MQ=255)
                                                                                                                                       tactctTAGTCTTTTGGTATTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGATAGCCAGACTTTTGTAATACTTCGTCAACTATTTCATGAATTTCTaaa    <  1:442773/139‑1 (MQ=255)
                                                                                                                                       tCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGATAGCCAGACTTTTGTAATACTTCATCAACTATTTCATGAATTTCTtaa    >  1:641413/1‑141 (MQ=255)
                                                                                                                                       tCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCAAATGTATTTTCACGTTCAAGCATTTCTCGATAGCCAGACTTTTGTAATACTTCATCAACTATTTCATGAATTTCTaaa    >  2:250664/1‑141 (MQ=255)
                                                                                                                                       tCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACATCCAGACCCTAATGTATTTTCACGTTCAAGCATTTCTCGATAGCCAGACTTTTGTAATACTTCATCAACTATTTCATGAATTTCTaaa    <  1:680014/141‑1 (MQ=255)
                                                                                                                                         cTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGATAGCCAGACTTTTGTAATACTTCATCAACTATTTCATGAATTTCTaaaaa  >  1:433590/1‑141 (MQ=255)
                                                                                                                                         cTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTATAATGGATTTTCACGTTCAAGCATTTCTCGATAGCCAGACTTTTGTAATACTTCATCAACTATTTCATGAATTTGTaaaaa  >  1:95686/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                          
GGAAATTCAAGGCCCTTAGCCGAGTGCATCGTCATTAGTGTTACGCCATTTTCAGTATCTGCCTCATCAATATCAGCTACTAACGATAAATCCGTTAAAAAGTTAATTAATGACTGTTCTTCTAATGGGGTATTTTCCTCATAGTCTTTTGGTACTGACATAAATTCATCGATGTTTTCTAATCTACTTCGAGATTCTAATGTATTTTCACGTTCAAGCATTTCTCGATAGCCAGACTTTTGTAATACTTCATCAACTATTTCATGAATTTCTAAAAA  >  CP000730/2050602‑2050879

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