Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,696,558 T→C 5.9% intergenic (+110/‑130) USA300HOU_2551 → / → copA hypothetical protein/P‑ATPase superfamily P‑type ATPase copper (Cu2+) transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,696,5580TC 5.9% 208.5 / 3.5 68intergenic (+110/‑130)USA300HOU_2551/copAhypothetical protein/P‑ATPase superfamily P‑type ATPase copper (Cu2+) transporter
Reads supporting (aligned to +/- strand):  ref base T (34/30);  new base C (2/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

GAATTTAACCCATCTGATTTCTTTTAAGTCACCTAAGAATTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATACCTATAGGGGGTACATTAGACGTGTAACAATGAATCACAGGGAGGCAATAATGTGGCT  >  CP000730/2696422‑2696693
                                                                                                                                        |                                                                                                                                       
gAATTTAACCCATCTGATTTCTTTTAAGTCACCTAAGAATTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTtgtg                                                                                                                                     <  2:13725/141‑1 (MQ=255)
  aTTTAACCCATCTGATTTCTTTTAAGTCACCTAAGAATTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACTACTTCTTTTATATAGATGCTAAGTAGACAGATTGTTGTGCa                                                                                                                                   >  2:467557/1‑141 (MQ=255)
  aTTTAACCCATCTGATTTCTTTTAAGTCACCTAAGAATTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCa                                                                                                                                   >  2:271986/1‑141 (MQ=255)
   tttAACCCATCTGATTTCTTTTAAGTCACCTAAGAATTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCaa                                                                                                                                  >  2:38232/1‑141 (MQ=255)
     tAACCCATCTGATTTCTTTTAAGTCACCTAAGAATTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATg                                                                                                                                >  1:689884/1‑141 (MQ=255)
     tAACCCATCTGATTTCTTTTAAGTCACCTAAGAATTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGCCGTGCaaac                                                                                                                                <  2:412422/141‑3 (MQ=255)
          cATCTGATTTCTTTTAAGTCACCTAAGAATTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGc                                                                                                                           >  1:387261/1‑141 (MQ=255)
              tGATTTCTTTTAAGTCACCTAAGAATTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGCCGTGCaaacattgcacaa                                                                                                                       <  1:479014/141‑12 (MQ=255)
               gATTTCTTTTAAGTCACCTAAGAATTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGc                                                                                                                      <  1:223956/141‑1 (MQ=255)
                 tttCTTTTAAGTCACCTAAGAATTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCaa                                                                                                                    <  2:388649/141‑1 (MQ=255)
                 tttCCTTTAAGTCACCTAAGAATTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCaa                                                                                                                    <  2:388425/141‑1 (MQ=255)
                   tCTTTTAAGTCACCTAAGAATTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAAtc                                                                                                                  <  1:590540/141‑1 (MQ=255)
                       ttAAGTCACCTAAGAATTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCCAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAAtctctc                                                                                                              >  2:640047/1‑141 (MQ=255)
                        tAAGTCACCTAAGAATTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATctctct                                                                                                             >  1:562624/1‑141 (MQ=255)
                         aaGTCACCTAAGAATTGCAAATCCAGAAGTCATTTAAGTTTTACTTTTCATTCATACATCCTTTATTATTAATTATGACTTTTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCtctt                                                                                                            >  1:244854/1‑141 (MQ=255)
                          aGTCACCTAAGAATTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCtcttt                                                                                                           >  2:477124/1‑141 (MQ=255)
                            tCACCTAAGAATTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCtcttttt                                                                                                         >  1:1479/1‑141 (MQ=255)
                             cACCTAAGAATTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCtctttttc                                                                                                        >  2:638917/1‑141 (MQ=255)
                              aCCTAAGACTTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGCTGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCtctttttct                                                                                                       >  2:336036/1‑141 (MQ=255)
                              aCCTAAGAATTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCtctttttct                                                                                                       >  1:304097/1‑141 (MQ=255)
                                cTAAGAATTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCtctttttcttt                                                                                                     <  1:89381/141‑1 (MQ=255)
                                cTAAGAATTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCtctttttcttt                                                                                                     <  2:583072/141‑1 (MQ=255)
                                  aaGAATTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTctttttcttttt                                                                                                   <  1:50922/141‑1 (MQ=255)
                                      aTTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATACTCATTACACCTTCCATTATATAGTTGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTaaaa                                                                                               <  1:494105/141‑1 (MQ=255)
                                         gCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTg                                                                                            >  1:443585/1‑141 (MQ=255)
                                                cAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGgtaaaag                                                                                     >  2:475989/1‑141 (MQ=255)
                                                  gAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGgtaaaagta                                                                                   >  1:417230/1‑141 (MQ=255)
                                                              gTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGAt                                                                       <  1:44740/141‑1 (MQ=255)
                                                                 ttACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGa                                                                    >  2:531383/1‑141 (MQ=255)
                                                                  tACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGAc                                                                   >  1:91431/1‑141 (MQ=255)
                                                                   aCCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACt                                                                  <  2:226497/141‑1 (MQ=255)
                                                                     cTTTCATTCATACATCCTTTAATATTCATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTta                                                                >  1:668668/1‑141 (MQ=255)
                                                                     cTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTta                                                                >  1:668719/1‑141 (MQ=255)
                                                                     cTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTta                                                                >  2:519780/1‑141 (MQ=255)
                                                                       ttcattcaTACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTtata                                                              >  1:324509/1‑141 (MQ=255)
                                                                        tcattcatACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTtatat                                                             <  1:373081/141‑1 (MQ=255)
                                                                        tcattcatACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTATTACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAGAACGCAACGATTGACTtatat                                                             <  2:22596/141‑1 (MQ=255)
                                                                          attcatACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATAc                                                           <  1:159105/141‑1 (MQ=255)
                                                                          attcatACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATAc                                                           <  1:313498/141‑1 (MQ=255)
                                                                           ttcatACATCCTTTACTATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATAcc                                                          <  2:244854/141‑1 (MQ=255)
                                                                               tACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATACCtata                                                      >  1:363039/1‑141 (MQ=255)
                                                                                aCATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATACCTATAg                                                     >  2:416257/1‑141 (MQ=255)
                                                                                 cATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATACCTATAgg                                                    <  1:128515/141‑1 (MQ=255)
                                                                                  aTCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATACCTATAggg                                                   <  2:1479/141‑1 (MQ=255)
                                                                                    ccTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATACCTATAggggg                                                 <  2:354475/141‑1 (MQ=255)
                                                                                        tAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCtctt                                                                                                            >  2:414411/1‑78 (MQ=255)
                                                                                        tAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCtctt                                                                                                            <  1:414411/78‑1 (MQ=255)
                                                                                           tatTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATACCTATAGGGGGTACATTa                                          >  1:632561/1‑141 (MQ=255)
                                                                                            attaattaCGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATACCTATAGGGGGTACATTAg                                         <  1:454231/141‑1 (MQ=255)
                                                                                                 ttaCGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATACCTATAGGGGGTACATTAGACgtg                                    <  1:687321/141‑1 (MQ=255)
                                                                                                   aCGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATACCTATAGGGGGTACATTAGACGTGTa                                  >  1:192998/1‑141 (MQ=255)
                                                                                                    ctaCTTCTTTTATATAGATGCTAAGTAGAGATATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATACCTATAGGGAGTACATTAGACGTGTaa                                 <  2:337172/139‑1 (MQ=255)
                                                                                                       cttcttTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATACCTATAGGGGGTACATTAGACGTGTAACaa                              >  1:404776/1‑141 (MQ=255)
                                                                                                          cttTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATACCTATAGGGGGTACATTAGACGTGTAACAATGa                           >  2:182179/1‑141 (MQ=255)
                                                                                                              tatataGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTGTTCTTTATAAAATTGTGAGAAGTAAAACGCAACGATTGACTGATATACCTATAGGGGGTACATTAGACGTGTAACAATGAATca                       >  1:420799/1‑141 (MQ=255)
                                                                                                                tataGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAAtctctc                                                                                                              <  2:9880/52‑1 (MQ=255)
                                                                                                                tataGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAAtctctc                                                                                                              >  1:9880/1‑52 (MQ=255)
                                                                                                                tataGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATACCTATAGGGGGTACATTAGACGTGTAACAATGAATcaca                     >  1:298287/1‑141 (MQ=255)
                                                                                                                      tGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATACCTATAGGGGGTACATTAGACGTGTAACAATGAATCACAGGGAgg               <  1:477124/141‑1 (MQ=255)
                                                                                                                       gCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATACCTATAGGGGGTACATTAGACGTGTAACAATGAATCACAGGGAGGc              <  1:338864/141‑1 (MQ=255)
                                                                                                                        gcgttaagagagATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATACCTATAGGGGGTACATTAGACGTGTAACAATGAATCACAGGGAGGCa             <  2:623354/135‑1 (MQ=255)
                                                                                                                        cTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATACCTATAGGGGGTACATTAGACGTGTAACAATGAATCACAGGGAGGCa             <  1:296399/141‑1 (MQ=255)
                                                                                                                        agaaaaagagagATTGCCGTGCAAACATTGCACAACAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATACCTATAGGGGGTACATTAGACGTGTAACAATGAATCACAGGGAGGCa             >  1:436069/7‑141 (MQ=25)
                                                                                                                         tAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATACCTATAGGGGGTACATTAGACGTGTAACAATGAATCACAGGGAGGCaa            <  1:38232/141‑1 (MQ=255)
                                                                                                                            gTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATACCTATAGGGGGTACATTAGACGTGTAACAATGAATCACAGGGAGGCaataa         >  1:626101/1‑141 (MQ=255)
                                                                                                                              agagagATTGCCGTGCAAACATTGCACAACAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATACCTATAGGGGGTACATTAGACGTGTAACAATGAATCACAGGGAGGCAATAAtg       >  2:50632/1‑141 (MQ=37)
                                                                                                                                agagaTTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATACCTATAGGGGGTACATTAGACGTGTAACAATGAATCACAGGGAGGCAATAAtgtg     >  2:597363/1‑141 (MQ=255)
                                                                                                                                   gaTTGTTGTGCAATCCCTGCACGGAAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTATAACGCACCGATTGACTTATATACCTATAGGGGGTACATTAGACGTGTATCCATGAATCACAGCGAGGCAATAATGTGCCt  <  2:167732/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
GAATTTAACCCATCTGATTTCTTTTAAGTCACCTAAGAATTGCAAATCCAGAAGTCATTTAAGTTTTACCTTTCATTCATACATCCTTTAATATTAATTACGACTTCTTTTATATAGATGCTAAGTAGAGAGATTGTTGTGCAATGTTTGCACGGCAATCTCTCTTTTTCTTTTTAAAATTGGTAAAAGTAAAACGCAACGATTGACTTATATACCTATAGGGGGTACATTAGACGTGTAACAATGAATCACAGGGAGGCAATAATGTGGCT  >  CP000730/2696422‑2696693

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