Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 587,957 G→T 58.6% A621S (GCG→TCG)  rpoB → DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730587,9570GT58.6% 37.6 / 66.5 67A621S (GCG→TCG) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (13/15);  new base T (29/10);  total (42/25)
Fisher's exact test for biased strand distribution p-value = 2.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ACAAGCAAACTCTAAATTAGATGAAAATGGTCGTTTCATGGATGATGAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATTTGTTGGTACAGGTATGGAACACGT  >  CP000730/587818‑588087
                                                                                                                                           |                                                                                                                                   
acaAGCAAACACTAAATTAGATGAAAATTGTCGTTTAATGGATGATGAAGTTGTACGTCGTTTCCGTGGTAACAACACAGTTAAGGCTAAAGTAAAAAAGGATTATAGGGATGTATGGCTGAAGCATGTTGTGTCagcagc                                                                                                                                    <  2:349850/141‑1 (MQ=255)
    caataCTCTAAATTAGATGGAAATGGTCGTTTCATGGATGATGAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACa                                                                                                                                <  2:354299/137‑1 (MQ=255)
      aaaatcaaaatcAGATGAAAATCGTAGATTAAAAGATAATGAAGTAATATGTAGATCAAGTGGAAACAATAAAGTTAAAAAAAAAAAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGc                                                                                                                              <  2:712496/141‑1 (MQ=255)
      aaaCTCTAAATTAGATGAAAATGGTCGTTTCATGGATGATGAAGTTGTATGTCGTTTCCGTGGTAGCAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGATGTTTCAGCAGCGACAGc                                                                                                                              >  2:131750/1‑141 (MQ=255)
      aaaCTCTAAATTAGATGAAAATGGTCGTTTCATGGATGATGAAGTTATATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGAATCGACAGc                                                                                                                              >  2:503416/1‑141 (MQ=255)
       atagctAAATTAGATGAAAATGGTCGTTTCATGGATGATGAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCa                                                                                                                             <  1:190550/137‑1 (MQ=255)
       aaCTCTAAATTAGATGAAAATGGTCGTTTCATGGATGATGAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCa                                                                                                                             >  2:485922/1‑141 (MQ=255)
        aCTCTAAATTAGATGAAAATGGACGTTTCATGGATGATGAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGAAAGCAt                                                                                                                            >  2:739920/1‑141 (MQ=255)
         ctctAAATTAGATGAAAATGGTCGTTTCATGGATGATGAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATg                                                                                                                           <  1:201907/141‑1 (MQ=255)
           ctAAATTAGATGAAACTGGTCGTTTCATGGATGATGAAGTTGTATGTCGTTTCCGTGGTCACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGGTGTTTCAGCATCGACAGCATGTa                                                                                                                         >  1:391183/1‑141 (MQ=255)
             aaaTTAGATGAAAATGGTCGTTTCATGGAGGATGAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTAtt                                                                                                                       <  2:595269/141‑1 (MQ=255)
                       aaaaTGGTCGTTTCATGGATGATGAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGCCAGCATGTATTCCAGTCTTAg                                                                                                            >  2:615768/1‑141 (MQ=255)
                       aaaaTGGTCGTTTCATGGATGATGAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCAACAGCATGTATTACATTCTTAg                                                                                                            >  2:17051/1‑141 (MQ=255)
                        aaaTGGTCGTTTCATGGATGATGAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGa                                                                                                           >  1:760400/1‑141 (MQ=255)
                            ggTCGTTTCATGGATGATGAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGTAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAat                                                                                                       >  2:405434/1‑141 (MQ=255)
                            ggTCGTTTCATGGATGATGAAGTTGTATGTCGTGTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAat                                                                                                       >  2:612712/1‑141 (MQ=255)
                             gTCGTTTCATGGATGATGAAGTTGTATGTCGTTTCCGTGGTAACAATCCAGTTATGGCTAACGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTATTAGAAAatg                                                                                                      >  1:86523/1‑141 (MQ=255)
                              tCGTTTCATGGATGATGAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAatga                                                                                                     >  1:401601/1‑141 (MQ=255)
                                gTTTCATGGATGATGAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAatgatg                                                                                                   <  1:387853/141‑1 (MQ=255)
                                  ttCATGGATGATGAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGAc                                                                                                 >  1:107398/1‑141 (MQ=255)
                                      tGGATGATGAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGAATATAGGGATGTATCGCCGAAGCAAGTTGTTTCAGAATCGACAGCATGTATTCCAATCTTAGACAAagatgaacgaa                                                                                             >  2:584395/1‑136 (MQ=255)
                                       ggATGATGAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCaaa                                                                                            <  1:730479/141‑1 (MQ=255)
                                          tgatgaAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCg                                                                                         >  2:127331/1‑141 (MQ=255)
                                              gaAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAAGCATGACTCAAACCGAGCa                                                                                     >  2:528378/1‑141 (MQ=255)
                                               aaGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCATGCAt                                                                                    >  2:525658/1‑141 (MQ=255)
                                                 gTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTg                                                                                  >  2:734846/1‑141 (MQ=255)
                                                 gTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTg                                                                                  <  1:210895/141‑1 (MQ=255)
                                                     tATGTCGTTTCCGTGGCAACAATACAGTTATGGCTAAACAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATgg                                                                              <  1:9049/141‑1 (MQ=255)
                                                      aTGTCGTGTCCGTGGTAAAAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATggg                                                                             <  2:33574/141‑1 (MQ=255)
                                                          cGTTTCCGTTGTAAAAATACAGTTATGGCTAAAGAAAAAATGTATATCATATATGTATCTCCTTAGCAAGATGTATTAGCAGAGAGAGCATTTATTCCATGCTTAGAATATGATGACTCAAACCGTGTATTGATGGGTGCg                                                                         <  2:590086/141‑1 (MQ=255)
                                                          cGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCg                                                                         <  1:43555/141‑1 (MQ=255)
                                                              tCCGTGGTAACAATGCAGTTATGGCTAAAGAAAAAATGGATTATATGTATGTATCGCCGAACCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACa                                                                     <  2:38179/141‑1 (MQ=255)
                                                               ccGTGGTAACAATACAGTTAGGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGCTGTCTCAAACCGTGCATTGATGGGTGCGAACAt                                                                    <  2:134057/141‑1 (MQ=255)
                                                                cGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATg                                                                   >  2:572308/1‑141 (MQ=255)
                                                                 gTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCCACCCGTGCATTGATGGGTGCGAACATga                                                                  <  2:661275/141‑2 (MQ=255)
                                                                    gTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGCCACAGCCTGTATTACATTTGTTGAAAATGATGACTCAAGCCGTGCATTGATGGGTGCGAACATGCAAc                                                               <  2:166997/141‑1 (MQ=255)
                                                                     tAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGGGAACATGAAACg                                                              >  2:367930/1‑141 (MQ=255)
                                                                     tAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACg                                                              <  1:63814/141‑1 (MQ=255)
                                                                     tAACAATACAGTGATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCGGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGCTGGGTGCGAACATTCAAcc                                                              >  2:618292/1‑140 (MQ=255)
                                                                           tACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGc                                                        >  1:606356/1‑141 (MQ=255)
                                                                             cAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAg                                                      <  1:670926/141‑1 (MQ=255)
                                                                                   tGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCCACGGCCAGCAGTACCtt                                                >  1:232358/1‑141 (MQ=255)
                                                                                    ggCTAAAGAAAAAATGAATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGGGCGAACATGCAACGCCAAGCCGTGCCttt                                               >  1:469371/1‑141 (MQ=255)
                                                                                      cTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGGACATGCCACGCCAAGCAGTGCCTTtga                                             >  2:703364/1‑141 (MQ=255)
                                                                                           gAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAAt                                        >  1:423968/1‑141 (MQ=255)
                                                                                            aaaaaaTGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATGCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATc                                       >  1:777117/1‑141 (MQ=255)
                                                                                              aaaaTGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCa                                     >  2:773037/1‑141 (MQ=255)
                                                                                              aaaaTGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGCGCCTTTTATGAATCCa                                     >  2:165054/1‑141 (MQ=255)
                                                                                                   ggATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGc                                >  1:438230/1‑141 (MQ=255)
                                                                                                   ggATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCCTTAATGAATCCAGAAGc                                >  2:258208/1‑141 (MQ=255)
                                                                                                   ggATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGc                                >  1:687456/1‑141 (MQ=255)
                                                                                                          atGGATGTATCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGAAGTGCCTGTGATGAATCCAGAAGCACCAttt                         <  2:760400/141‑1 (MQ=255)
                                                                                                           tggatggaTCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGAGCGAACATGCAACGTCAAACAGTGCCGTTGATGAATCCAGAAGAACCATttg                        >  2:496014/1‑141 (MQ=255)
                                                                                                            ggATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCCACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATgtgt                       >  2:472681/1‑141 (MQ=255)
                                                                                                            ggATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTAAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATttgt                       >  2:482802/1‑141 (MQ=255)
                                                                                                             gATGTATCGCGGAAGCCAGTTGTTTCAGCAGCGACAGCATGTATTCTTTTTCTTTTGAAAGATGAATAAAACCGTTCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATttgtt                      <  1:149357/141‑1 (MQ=255)
                                                                                                                  aTCGCCGAAGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATTTGTTGGTAc                 >  1:484113/1‑141 (MQ=255)
                                                                                                                   tCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATTTGTTGGTACa                <  1:43585/141‑1 (MQ=255)
                                                                                                                         aaGCAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGCACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAAAAGCACCACTTGTTGGCTCCGGTAtc          >  2:614232/1‑140 (MQ=255)
                                                                                                                          aggaaGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAAGGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATTTGTTGGTACAGGTATgg         <  2:184512/138‑1 (MQ=255)
                                                                                                                            cAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATTTGTTGGTACAGGTATGGaa       >  1:69371/1‑141 (MQ=255)
                                                                                                                            cAAGTTGTTTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATTTGTTGGTACAGGTATGGaa       <  1:91974/141‑1 (MQ=255)
                                                                                                                            cAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATTTGTTGGTACAGGTATGGaa       <  1:384851/141‑1 (MQ=255)
                                                                                                                            cAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACCCAAACCGTGAATTGATGGTTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATTTGTTGGTACAGGTATGGaa       <  1:775648/141‑1 (MQ=255)
                                                                                                                               gttgttTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATTTGTTGGTACAGGTATGGCacac    >  1:591785/1‑141 (MQ=255)
                                                                                                                               gttgttTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATTTGTTGGTACAGGGATGGAaccc    >  2:110065/1‑139 (MQ=255)
                                                                                                                               gttgttTCAGCATCGACAGCATGTATTCCATTCTTAGAAAATCATGCGCCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATTTGTTGGTAGAGGTATGGAacac    >  1:615382/1‑141 (MQ=255)
                                                                                                                                 tgttTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATTTGTTGGTAAAGGTATGGAACACGt  >  2:414741/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                   
ACAAGCAAACTCTAAATTAGATGAAAATGGTCGTTTCATGGATGATGAAGTTGTATGTCGTTTCCGTGGTAACAATACAGTTATGGCTAAAGAAAAAATGGATTATATGGATGTATCGCCGAAGCAAGTTGTTTCAGCAGCGACAGCATGTATTCCATTCTTAGAAAATGATGACTCAAACCGTGCATTGATGGGTGCGAACATGCAACGTCAAGCAGTGCCTTTGATGAATCCAGAAGCACCATTTGTTGGTACAGGTATGGAACACGT  >  CP000730/587818‑588087

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: