Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,872,089 G→T 91.9% S464Y (TCT→TAT)  gyrB ← DNA gyrase, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,872,0890GT91.9% 121.7 / 7.7 50S464Y (TCT→TAT) gyrBDNA gyrase, subunit B
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (21/25);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01

GCTGAGCGATGTAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATC  >  CP009273/3871886‑3872351
                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                      
gcTGAGCGATGTAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGCAACGACCGATACCACAGCCAAGCGCGGTGATAAGCTTCGCCACTTCCTCATAAGAGAACATCATAACGAAGCGCGCTGTCTCGACGTTGCGGAATTTACCCTTCACGGGTNAGATCGCCTg                                                                                                                                                                                                  >  2:88591/1‑274 (MQ=255)
                                               agATAGAGGAATGTCTACAGAAGCGTATGAATGTGCGATCAGTTGTCGTCCGCATCGGTCGTGATGATGATGCTGTGATAACGCAGTTTGTNCGGGTTGTACTCTTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCGAGCTTGGGAGAGCAGCCCGCGGAGTcccc                                                                                                                                                  <  1:569380/274‑1 (MQ=255)
                                                   agaagaAGGTTAACAGCAGCGTGCGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCtt                                                                                                                                                <  1:499497/273‑1 (MQ=255)
                                                        aagGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGACCCTGCTTCGCAGAGCCGCCCGCGGAATCCCCTCCCACCAg                                                                                                                                         >  1:35804/1‑275 (MQ=255)
                                                                      cGTACGAATGTGCGAGCCGTCGACGTCGGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGAtttt                                                                                                                                                                                                                          >  1:479842/1‑180 (MQ=255)
                                                                      cGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGc                                                                                                                                                                            >  1:761891/1‑226 (MQ=255)
                                                                          cGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTAccc                                                                                                                                                                                                                      <  2:479839/180‑1 (MQ=255)
                                                                          cGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCg                                                                                                                                                                        <  2:761885/226‑1 (MQ=255)
                                                                           gAATGTGCGGGCGGTGGAGTTCCGCATCGGGCAGGATGATGATGGTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGAGACCACAGCCAAGCGCGGTGAGAAGCGTCGCCACTGCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCGTCAGCGGCAGAATCGCCTGGTTCTTGCGGTGACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAgcgc                                                                                                                       <  1:211722/274‑1 (MQ=255)
                                                                                tgCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTAcc                                                                                                                                                                                                                       >  1:361809/1‑173 (MQ=255)
                                                                                tgCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCg                                                                                                                                                                        >  2:287254/1‑220 (MQ=255)
                                                                                tgCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGGTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCg                                                                                                                                                                                         >  1:480429/1‑203 (MQ=255)
                                                                                  cGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCagag                                                                                                                                                                   >  2:790704/1‑223 (MQ=255)
                                                                                    aGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTc                                                                                                                                                                                                                   <  2:361806/173‑1 (MQ=255)
                                                                                    aGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCaga                                                                                                                                                                    <  1:287257/220‑1 (MQ=255)
                                                                                    aGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTa                                                                                                                                                                                     <  2:480426/203‑1 (MQ=255)
                                                                                      ccGTCGACGTCCGCATCGGTCATGATGCTGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAgccgc                                                                                                                                                               <  1:790710/223‑1 (MQ=255)
                                                                                            gcgaccGCATTGGGCATGGTGTTGATGCGGTGATAACTCAGTTGGTCCGGGTTGTACTCGTAACGACCGATACCACCGCGAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTGGAGGATTTTACCCTTCAGCGGCAGGCTCGCCTGGTTCTTGCGGTTACGCCCCTGCTGCGCAGAGCCGCCCGCGGAGTCCCCTCCCACCAGGGACAGTTCGGAAAGCGCGGGATCGAGTTCCTGGCa                                                                                                     <  2:35804/271‑1 (MQ=255)
                                                                                                           atgatgatgatgCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTAc                                                                                                                                                                                    >  1:771655/1‑181 (MQ=255)
                                                                                                               tgatgatgCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGccc                                                                                                                                                                                <  2:771649/181‑1 (MQ=255)
                                                                                                                       cTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTc                                                                                                                                                                                              >  2:104843/1‑159 (MQ=255)
                                                                                                                        tgtgATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTGCAGCGGCAGAATCGCCTGGTTCTTGAGGTTACGCCCCTGCTTCGCAGAGCAGCCCGCGGAGTCCCCGTccacca                                                                                                                                          >  1:153317/1‑210 (MQ=255)
                                                                                                                        tggatgaaCNGAGTTTGTCCGGGGTGTACGCGTCACGACCTATCCCACAGACAGGCGCGGTGATAAGCTTCGCCACTTCCTGATAAGAGATCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCCGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCt                                                                         <  2:333615/269‑1 (MQ=255)
                                                                                                                          tgATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGataag                                                                                                                                                                                                                                                                     >  1:720231/1‑85 (MQ=255)
                                                                                                                          tgATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCATTCTGCCAGTTTGCCCGGCAGGCCCGCTa                                                                        >  1:333618/1‑274 (MQ=255)
                                                                                                                           gATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGc                                                                                                                                                                                          <  1:104843/159‑1 (MQ=255)
                                                                                                                              aaCGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAagagag                                                                                                                                                                                                                                                                 <  2:720225/85‑1 (MQ=255)
                                                                                                                                  ccagttGTCCGGGTTGTCATCCTCCCGACCGATACCACGGCCAGGCGCGGGTATAAGCTTCGCCACTTCCTGATCAGAGAGCATCTTATCGAAGCGCGCTTGCTCGACGTTGAGTATTTTTCCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAgcg                                                               <  2:18797/271‑1 (MQ=255)
                                                                                                                                             ggTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTACCCTTCCCCCAGGTACAGTTCGGAAAGCGCGGGATCGCGTTCCCGCAAGTCTGCCATTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTAcca                                                    >  1:420772/1‑273 (MQ=255)
                                                                                                                                                 gTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGTATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTACGTAGAGCGCACCTTTAcggctg                                                 >  1:18797/1‑272 (MQ=255)
                                                                                                                                                       gTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGa                                                                                                                                                                                                                                       >  1:154539/1‑86 (MQ=255)
                                                                                                                                                           cgcccGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAgcg                                                               <  1:168440/247‑1 (MQ=255)
                                                                                                                                                           cGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGtt                                                                                                                                                                                                                                   <  2:154538/86‑1 (MQ=255)
                                                                                                                                                            gACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAgccgcc                                                                                                                                                              >  1:800943/1‑154 (MQ=255)
                                                                                                                                                                gATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCcgcg                                                                                                                                                          <  2:800937/154‑1 (MQ=255)
                                                                                                                                                                        aGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCttt                                                                                                                                                                                                                                            >  1:419677/1‑64 (MQ=255)
                                                                                                                                                                          ccAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCa                                                                                                                                                                                                                  >  2:239744/1‑88 (MQ=255)
                                                                                                                                                                            aaGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCg                                                                                                                                                                                                                                        <  2:419674/64‑1 (MQ=255)
                                                                                                                                                                              gcgcGTTGATAAGCGTCGCCACTTCCTGTTAAGAGCGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTca                   <  1:88591/275‑1 (MQ=255)
                                                                                                                                                                              gcgcGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCgg                                                                                                                                                                                                              <  1:239747/88‑1 (MQ=255)
                                                                                                                                                                                   tggttAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCATCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCACTTcacgggc              <  1:790978/271‑1 (MQ=255)
                                                                                                                                                                                             tCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGtt                                                                                                                                                                                                                                   <  2:552857/52‑1 (MQ=255)
                                                                                                                                                                                              cGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAgccgcc                                                                                                                                                              >  1:222368/1‑120 (MQ=255)
                                                                                                                                                                                              cGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCt                                                                                                                                                 >  1:559424/1‑133 (MQ=255)
                                                                                                                                                                                              cGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCg                                                                                    >  1:155851/1‑194 (MQ=255)
                                                                                                                                                                                                ccACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCCACTTTACGGCGGGTCATTTCCCGCGCGCGATGCGCCGCCTGACGGGCACGGgcagcacc  >  1:488546/1‑272 (MQ=255)
                                                                                                                                                                                                  aCTTCCTGATAAGAGAGCNTCTTATCGNAGCGCGCTTTCTAGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCCGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAgcg                                                               <  1:548962/211‑1 (MQ=255)
                                                                                                                                                                                                  aCTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCcgcg                                                                                                                                                          <  2:222365/120‑1 (MQ=255)
                                                                                                                                                                                                  aCTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTcca                                                                                                                                             <  2:559421/133‑1 (MQ=255)
                                                                                                                                                                                                  aCTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAAGCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAg                                                                                <  2:155850/194‑1 (MQ=255)
                                                                                                                                                                                                        tGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTccacc                                                                                                                                           >  1:432938/1‑129 (MQ=255)
                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                      
GCTGAGCGATGTAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATC  >  CP009273/3871886‑3872351

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