Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000731 25,663 Δ1 bp 100% coding (316/333 nt) USA300HOU_pUSA300HOUMR0031 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073125,6630A.99.2% 6143.6 / 24.0 1241coding (316/333 nt)USA300HOU_pUSA300HOUMR0031hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/8);  new base . (672/559);  total (674/567)
Fisher's exact test for biased strand distribution p-value = 1.87e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

Only 100 of 1416 total aligned reads displayed.

ATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAA  >  CP000731/25573‑25753
                                                                                          |                                                                                          
aTTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATa                                                                                            >  1:348646/1‑91 (MQ=255)
 ttGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaa                                                                                           <  2:2368295/91‑1 (MQ=255)
  tGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaa                                                                                          >  1:3908681/1‑91 (MQ=255)
   gAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaa                                                                                         >  2:5522100/1‑91 (MQ=255)
   gAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaa                                                                                         <  2:3240615/91‑1 (MQ=255)
     aGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaa                                                                                       <  1:5516579/91‑1 (MQ=255)
     aGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaa                                                                                       >  1:4954419/1‑91 (MQ=255)
     aGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaa                                                                                       >  2:4177495/1‑91 (MQ=255)
     aGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaa                                                                                       >  2:2713823/1‑91 (MQ=255)
      ggATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaat                                                                                      >  1:1489042/1‑90 (MQ=255)
          cTTGGTTTTCTCATTATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATAc                                                                                 >  2:3986031/1‑91 (MQ=255)
          cTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATAc                                                                                 >  2:2263198/1‑91 (MQ=255)
          cTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACGATACTACAAAT‑AAAAAATATAc                                                                                 <  2:4298690/91‑1 (MQ=255)
           ttGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACa                                                                                <  2:3870261/91‑1 (MQ=255)
           ttGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACa                                                                                <  1:3365856/91‑1 (MQ=255)
           ttGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACa                                                                                >  2:3377624/1‑91 (MQ=255)
             ggTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATc                                                                              >  1:5179214/1‑91 (MQ=255)
             ggTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATc                                                                              >  2:1067374/1‑91 (MQ=255)
               ttttCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTa                                                                            >  1:4271568/1‑91 (MQ=255)
               ttttCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTa                                                                            >  1:3775947/1‑91 (MQ=255)
                 ttCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTaaa                                                                          >  1:2000543/1‑91 (MQ=255)
                   ctcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAgg                                                                        <  2:1730199/91‑1 (MQ=255)
                   ctcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAgg                                                                        >  2:1462788/1‑91 (MQ=255)
                    tcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGt                                                                       <  2:581487/91‑1 (MQ=255)
                     cAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGtt                                                                      >  1:3087164/1‑91 (MQ=255)
                      aaTATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTa                                                                     >  1:1111563/1‑91 (MQ=255)
                        tatCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTaaa                                                                   >  2:80928/1‑91 (MQ=255)
                          tCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAAtt                                                                 <  1:2560161/91‑1 (MQ=255)
                          tCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAAtt                                                                 <  1:1697500/91‑1 (MQ=255)
                          tCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAAtt                                                                 <  1:1575090/91‑1 (MQ=255)
                            cAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTaa                                                               >  1:3011594/1‑91 (MQ=255)
                            cAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTaa                                                               >  1:1892902/1‑91 (MQ=255)
                            cAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTACACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTaa                                                               >  1:4676659/1‑91 (MQ=255)
                             aGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAAt                                                              >  2:4551474/1‑91 (MQ=255)
                              gAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATg                                                             <  1:3201936/91‑1 (MQ=255)
                              gAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATg                                                             >  2:3665576/1‑91 (MQ=255)
                               aaaTTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGa                                                            <  1:2271319/91‑1 (MQ=255)
                                 aTTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGaaa                                                          <  1:2205459/91‑1 (MQ=255)
                                 aTTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGaaa                                                          <  1:5284912/91‑1 (MQ=255)
                                  ttttCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGaaaa                                                         <  1:2107469/91‑1 (MQ=255)
                                  ttttCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGaaaa                                                         <  2:5604288/91‑1 (MQ=255)
                                    ttCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAAt                                                       <  2:1590550/91‑1 (MQ=255)
                                    ttCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAAt                                                       <  2:3718919/91‑1 (MQ=255)
                                    ttCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAAt                                                       <  2:5262838/91‑1 (MQ=255)
                                     tCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATa                                                      >  2:1137675/1‑91 (MQ=255)
                                      cTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAg                                                     >  2:3572861/1‑91 (MQ=255)
                                       tttATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGa                                                    >  1:277506/1‑91 (MQ=255)
                                       tttATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATAAAAAATAGa                                                    <  2:737414/91‑1 (MQ=255)
                                        ttATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGaa                                                   <  2:1380282/91‑1 (MQ=255)
                                            aaTTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTag                                               <  2:4740589/91‑1 (MQ=255)
                                             aTTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTaga                                              <  1:1174683/91‑1 (MQ=255)
                                                gTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGaa                                           >  2:4937143/1‑91 (MQ=255)
                                                gTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGaa                                           <  1:5068429/91‑1 (MQ=255)
                                                gTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGaa                                           <  1:2440144/91‑1 (MQ=255)
                                                 tAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAat                                          >  2:4645911/1‑91 (MQ=255)
                                                 tAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAat                                          >  2:1940775/1‑91 (MQ=255)
                                                  aGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAata                                         >  2:5132379/1‑91 (MQ=255)
                                                     aCTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAg                                      >  1:4574358/1‑91 (MQ=255)
                                                     aCTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAg                                      >  2:3061424/1‑91 (MQ=255)
                                                      cTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGa                                     <  1:954830/91‑1 (MQ=255)
                                                       tATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGaa                                    <  1:1388189/91‑1 (MQ=255)
                                                       tATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGaa                                    >  1:3991876/1‑91 (MQ=255)
                                                           tCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACa                                >  2:3456090/1‑91 (MQ=255)
                                                             ttATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaa                              >  2:892094/1‑91 (MQ=255)
                                                             ttATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaa                              <  2:3164156/91‑1 (MQ=255)
                                                             ttATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaa                              >  1:4382732/1‑91 (MQ=255)
                                                             ttATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaa                              <  2:2494061/91‑1 (MQ=255)
                                                              tATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaaa                             >  2:2043464/1‑91 (MQ=255)
                                                              tATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaaa                             <  1:5639813/91‑1 (MQ=255)
                                                               aTGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaaaa                            >  2:1252861/1‑91 (MQ=255)
                                                               aTGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaaaa                            <  2:4441401/91‑1 (MQ=255)
                                                               aTGCAGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaaaa                            <  1:2330873/91‑1 (MQ=255)
                                                                tGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAg                           >  1:808729/1‑91 (MQ=255)
                                                                tGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAg                           <  1:4093745/91‑1 (MQ=255)
                                                                tGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAg                           >  1:372099/1‑91 (MQ=255)
                                                                 gCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAgg                          <  2:2946044/91‑1 (MQ=255)
                                                                 gCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAgg                          <  1:478204/91‑1 (MQ=255)
                                                                  cTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGa                         >  1:1298507/1‑91 (MQ=255)
                                                                    gTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAAt                       <  1:1804042/91‑1 (MQ=255)
                                                                       tAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTtc                    <  1:686119/91‑1 (MQ=255)
                                                                       tAATTACTATAATACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTtc                    <  2:670843/91‑1 (MQ=255)
                                                                         aTTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAAAGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTtctc                  <  1:2700067/91‑1 (MQ=255)
                                                                         aTTACTATAATACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAAAAAAAAAAAAAAAAAGGAATTtctc                  <  2:1641011/91‑1 (MQ=255)
                                                                         aTTAATATACTAAAAAT‑AAAAAATATACATCTAAAGATTAAATTAATGAAAAATAGAATTAAAGAAAAAAAAAAAAAAAAAGGAATTtctc                  <  2:2812265/91‑1 (MQ=255)
                                                                          ttACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCa                 >  2:1828322/1‑91 (MQ=255)
                                                                          ttACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCa                 <  1:4188266/91‑1 (MQ=255)
                                                                           tactatactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCaa                >  2:5032400/1‑91 (MQ=255)
                                                                            actatactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAg               <  1:5404423/91‑1 (MQ=255)
                                                                            actatactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAg               <  2:4823008/91‑1 (MQ=255)
                                                                            actatactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAg               <  2:408490/91‑1 (MQ=255)
                                                                            actatactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAg               <  1:361495/91‑1 (MQ=255)
                                                                            actatacgaCAAAT‑AATAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATAGAGAAAACAAAAAGGAATTTCTCAAg               <  2:5379956/91‑1 (MQ=255)
                                                                            aatatactaCAAAT‑AAAAAATATACATCTAAAAGTTAAATTAATGAAAAATAGAATTAGAGAAAATAGAAAACAAAAAGGAATTTCTCAAg               <  2:2586641/89‑1 (MQ=255)
                                                                              tatactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGaa             >  2:215351/1‑91 (MQ=255)
                                                                               atactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGaaa            >  1:4504674/1‑91 (MQ=255)
                                                                                tactacAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAAc           <  1:2934182/91‑1 (MQ=255)
                                                                                 actacAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACt          >  1:101913/1‑91 (MQ=255)
                                                                                   aacAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACttt        <  2:3935870/90‑1 (MQ=255)
                                                                                    acAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTa       <  1:61801/91‑1 (MQ=255)
                                                                                      aaaaaaaaaaATATACATCTAAAAGTTAAATTAATGAAAAATAGATTTAGAGAATTAAAAAAAAAAAAAGGAATTTCTCAAGAAACTTTaa      <  2:5441135/91‑1 (MQ=255)
                                                                                          |                                                                                          
ATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAA  >  CP000731/25573‑25753

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