Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,125,862 C→T 7.3% S700S (AGC→AGT)  katG → catalase‑peroxidase HPI, heme b‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,125,8620CT 7.3% 147.6 / 2.2 55S700S (AGC→AGT) katGcatalase‑peroxidase HPI, heme b‑containing
Reads supporting (aligned to +/- strand):  ref base C (27/23);  new base T (2/2);  total (29/26)
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.42e-01

ACGAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGT‑AATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGA  >  CP009273/4125714‑4126010
                                                                                                                                                    |                                                                                                                                                     
aCGAGGGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGc                                                                                                                                                       <  2:667535/149‑1 (MQ=255)
  gAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGACGAGAACACCAAACCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGa                                                                                                                                                     <  2:673920/149‑1 (MQ=255)
       gAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGccc                                                                                                                                                >  2:57643/1‑149 (MQ=255)
                      aTCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGGTCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTa                                                                                                                                 >  2:132219/1‑149 (MQ=255)
                         gAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAGTTTACGGCCAGCCGTGCGGGTCGGGTGTTTGGTTCTAACTCCGTCCGGCGTGCGGTGGCGGAAGGTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAg                                                                                                                              >  2:7776/1‑149 (MQ=255)
                            agagCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGTGATGCCCACGAGAAGTTTGTTAAAGACt                                                                                                                           >  2:88533/1‑149 (MQ=255)
                            agagCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGTGATGCCCACGAGAAGTTTGTTAAAGACt                                                                                                                           >  2:839759/1‑149 (MQ=255)
                               gCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTGTGTTAAAGACTTCg                                                                                                                        >  1:306163/1‑149 (MQ=255)
                                cTGTTCGAAGGCCGTGACCGTGAGACCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGAAAGTTTACGCCAGTAGCGATGCCCACGAGGAGTTTGTTAAAGACATcag                                                                                                                       >  1:596389/1‑147 (MQ=255)
                                cTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGt                                                                                                                       >  2:433372/1‑149 (MQ=255)
                                 tGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTg                                                                                                                      >  2:758871/1‑149 (MQ=255)
                                     cGAAGGCCGTGACCGTGAAAGCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTggcgg                                                                                                                  <  1:141895/149‑1 (MQ=255)
                                     cGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTggcgg                                                                                                                  >  2:521963/1‑149 (MQ=255)
                                       aaGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTggcggcg                                                                                                                >  2:290523/1‑148 (MQ=255)
                                         ggCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATg                                                                                                              <  1:684583/149‑1 (MQ=255)
                                            cgtgaccgtgaAACCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCCGTAGCGATGCCCACGAGAAGTTTGTTAAAGGCTTTGTGGCGGCATGGGt                                                                                                           >  2:426851/1‑149 (MQ=255)
                                            cgtgaccgtgaAACCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGt                                                                                                           >  2:426833/1‑149 (MQ=255)
                                              tgaccgtgaAACCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGGGa                                                                                                         >  1:644795/1‑149 (MQ=255)
                                                      aaaCCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGtgatg                                                                                                 >  2:617410/1‑149 (MQ=255)
                                                       aaCCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGtgatga                                                                                                >  2:662145/1‑149 (MQ=255)
                                                         ccGGCGAAGTGAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGCCTTCGTGGCGGCATGGGTGAAAGTGATGAAc                                                                                              >  2:225122/1‑149 (MQ=255)
                                                           ggCGAAGTGAAATTGTCGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCt                                                                                            <  2:585516/149‑1 (MQ=255)
                                                            gCGAAGTGAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAGGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGGCTGCGTGGCGGCATGGGGGAAAGTGAGAAAcccc                                                                                           >  1:138856/1‑147 (MQ=255)
                                                            gCGAAGTGAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTc                                                                                           <  2:794688/149‑1 (MQ=255)
                                                            gCGAAGTGAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTc                                                                                           <  1:390257/149‑1 (MQ=255)
                                                              gaggTGAAATTTACGGCCGGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGAGGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGa                                                                                         <  1:760823/146‑1 (MQ=255)
                                                                   gAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGGTAAAGACTTCGTGGCGGCATGGGGGGAAGTGAGGAACCTCGACCGtt                                                                                    >  2:275646/1‑149 (MQ=255)
                                                                    aaaTTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGttt                                                                                   >  2:55987/1‑149 (MQ=255)
                                                                    aaaGTTACGGCCAGCCGTGCGGATCTGGTGGTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGttt                                                                                   <  2:415153/149‑1 (MQ=255)
                                                                     aaTTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGAAGAACCTCGACCGTTTc                                                                                  <  2:549438/149‑1 (MQ=255)
                                                                             gCCAGCCGGGCGGAGCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTCAAGATTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACctgct                                                                          <  2:170878/149‑1 (MQ=255)
                                                                                   cgggcgcATCTGGTGTTTGGTTCTAACTCCGTCCGGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCCGT‑AATc                                                                   <  1:532590/146‑1 (MQ=255)
                                                                                        ggATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGT‑AATCTGAcc                                                              <  1:521963/149‑1 (MQ=255)
                                                                                           tCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGT‑AATCTGACCCCg                                                           <  1:103279/149‑1 (MQ=255)
                                                                                           tCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGT‑AATCTGACCCCg                                                           <  1:758871/149‑1 (MQ=255)
                                                                                            cTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGT‑AATCTGACCCCGt                                                          <  2:740288/149‑1 (MQ=255)
                                                                                            cTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGT‑AATCTGACCCCGt                                                          <  2:518921/149‑1 (MQ=255)
                                                                                                 gtTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGTGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGATCTCCTTT‑AAGCTGACCACGTTCAGc                                                     <  1:839759/149‑1 (MQ=255)
                                                                                                 gtTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGTGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGT‑AATCTGACCCCGTTCAGc                                                     <  1:88533/149‑1 (MQ=255)
                                                                                                  tttGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGT‑CATCTGACCCCGTTCAGCg                                                    >  2:15107/1‑149 (MQ=255)
                                                                                                     ggTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGT‑AATCTGACCCCGTTCAGCGgct                                                 >  2:662783/1‑149 (MQ=255)
                                                                                                      gttctggagccGTCCTGCGTGGGGTGGCGGAAGTTTACGCCAGTAGAGCTGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGT‑AATCTGACCCCGTTCAGCGgctg                                                <  2:235251/140‑1 (MQ=255)
                                                                                                       ttCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGT‑AATCTGACCCCGTTCAGCGgctgc                                               >  2:471080/1‑149 (MQ=255)
                                                                                                           aaCTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGT‑AATCTGACCCCGTTCAGCGGCtgcgtgc                                           >  1:854318/1‑149 (MQ=255)
                                                                                                             cTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACNAGAAGTTTGTTAAAGACTTCGCGGCGGCGTGGCTGGAAGTGATGAACCTCGACCGCTCCGACCTCCTGT‑AACCTGACCCCCATCAGCGAGcgcttgcta                                         >  1:450047/1‑148 (MQ=255)
                                                                                                              tCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGGTATACCACGACCGTTTCGACCTGCTGT‑AATCTGACCCCGTTCAGCGGCTGCTTGCTgg                                        >  2:202377/1‑149 (MQ=255)
                                                                                                                      gcgtgcgGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGT‑AATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCt                                <  2:543996/149‑1 (MQ=255)
                                                                                                                       cgtgcgGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGT‑AATCTGACCCCGTTCTGCGGCTGCTTGCTGGCAGTCGCTg                               <  2:439064/149‑1 (MQ=255)
                                                                                                                         tgcgGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGGACCTCGACCGTTTCGACCTGCTGT‑AATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTgca                             >  1:562501/1‑147 (MQ=255)
                                                                                                                             gTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGT‑AATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGtt                         >  2:419459/1‑149 (MQ=255)
                                                                                                                                       tttAGGCCAGGTGCGATGCCCACGAGAAGTCTGTTAAAGTCTTCGTGGCAGCATGGGTGACAGTGCTTTTCCCGGGCCCTTTTGACC‑CTTGTGGATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGc               <  2:421025/149‑1 (MQ=255)
                                                                                                                                       tttACGGCAGTAGCGATGCGCACGAGAAGTTTGTGAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGT‑AATCTGACCCCGTTCAGCGGCCGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGc               <  1:701604/149‑1 (MQ=255)
                                                                                                                                       tttACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGT‑AATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGc               <  2:273943/149‑1 (MQ=255)
                                                                                                                                              cAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGT‑AATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCGGCGTATAGt        >  2:647949/1‑149 (MQ=255)
                                                                                                                                                  aGCGATGTCCACGATAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGT‑AATCTGACTCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGc    <  1:415090/149‑1 (MQ=255)
                                                                                                                                                  aGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGT‑AATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGc    <  1:471080/149‑1 (MQ=255)
                                                                                                                                                    cGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGCAAGTGCTGAACCTCGCCCGTTTCGACCTGCTGT‑AATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGCACGTTCCTTGCCAGCGTCTCGTCGGcga  >  1:383610/1‑149 (MQ=255)
                                                                                                                                                    cGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGT‑AATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGcga  >  2:442966/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACGAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGATCTGGTGTTTGGTTCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGT‑AATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGA  >  CP009273/4125714‑4126010

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