Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000731 25,669 (A)7→6 coding (322/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% 4213.5 / 24.6 919coding (316/333 nt)USA300HOU_pUSA300HOUMR0031hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/5);  new base . (514/393);  total (518/401)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1061 total aligned reads displayed.

ATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAA  >  CP000731/25573‑25753
                                                                                          |                                                                                          
aTTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATa                                                                                            >  1:1825054/1‑91 (MQ=255)
 ttGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaa                                                                                           >  2:35568/1‑91 (MQ=255)
 ttGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaa                                                                                           <  2:143339/91‑1 (MQ=255)
   gAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaa                                                                                         >  1:3538635/1‑91 (MQ=255)
   gAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaa                                                                                         >  1:872386/1‑91 (MQ=255)
   gAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaa                                                                                         >  1:2985844/1‑91 (MQ=255)
     aGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaa                                                                                       >  1:1916250/1‑91 (MQ=255)
     aGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaa                                                                                       >  2:1941078/1‑91 (MQ=255)
     aGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaa                                                                                       <  1:17534/91‑1 (MQ=255)
       gATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaata                                                                                     <  1:1001485/91‑3 (MQ=255)
       gATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaata                                                                                     <  2:1620714/91‑3 (MQ=255)
       gATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaata                                                                                     >  1:1312696/1‑89 (MQ=255)
        aTCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaatat                                                                                    <  2:2639911/91‑4 (MQ=255)
          cTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATAc                                                                                 >  2:176499/1‑91 (MQ=255)
          cTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATAc                                                                                 >  1:311383/1‑91 (MQ=255)
          cTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATAc                                                                                 >  2:2903399/1‑91 (MQ=255)
           ttGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACa                                                                                <  1:1228102/91‑1 (MQ=255)
             ggTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATc                                                                              >  1:607315/1‑91 (MQ=255)
             ggTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATc                                                                              >  2:2328951/1‑91 (MQ=255)
             ggTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATc                                                                              >  2:584474/1‑91 (MQ=255)
             ggTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTACTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATc                                                                              >  1:2883560/1‑91 (MQ=255)
             ggTATTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAGATACATc                                                                              >  2:3214124/1‑91 (MQ=255)
                tttCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTaa                                                                           <  2:77470/91‑1 (MQ=255)
                tttCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTaa                                                                           >  2:2268489/1‑91 (MQ=255)
                   ctcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAgg                                                                        <  2:3038022/91‑1 (MQ=255)
                   ctcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAgg                                                                        >  2:377601/1‑91 (MQ=255)
                   ctcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAgg                                                                        >  2:2599691/1‑91 (MQ=255)
                   ctcAATATCCAAAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAgg                                                                        >  1:918392/1‑91 (MQ=255)
                    tcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGt                                                                       <  1:2352999/91‑1 (MQ=255)
                    tcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGt                                                                       >  1:417419/1‑91 (MQ=255)
                    tcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGt                                                                       <  1:362890/91‑1 (MQ=255)
                     cAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGtt                                                                      >  2:1554318/1‑91 (MQ=255)
                     cAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGtt                                                                      >  2:3423309/1‑91 (MQ=255)
                     cAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGtt                                                                      >  2:2485720/1‑91 (MQ=255)
                          tCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAAtt                                                                 >  1:2605081/1‑91 (MQ=255)
                            cAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTaa                                                               >  2:1971438/1‑91 (MQ=255)
                            cAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTaa                                                               >  2:1165918/1‑91 (MQ=255)
                              gAAATTTTCTTTATAATTGTAGCACTATTTCTTCTGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATg                                                             >  1:664929/1‑91 (MQ=255)
                              gAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATg                                                             <  1:1482284/91‑1 (MQ=255)
                              gAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATg                                                             >  2:1134073/1‑91 (MQ=255)
                               aaaTTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGa                                                            >  1:3420208/1‑91 (MQ=255)
                               aaaTTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGa                                                            <  2:2974103/91‑1 (MQ=255)
                                    ttCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAAt                                                       <  2:927861/91‑1 (MQ=255)
                                    ttCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTAAAAAAAAAAAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAAt                                                       <  2:2718173/91‑1 (MQ=255)
                                     tCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATa                                                      >  1:273840/1‑91 (MQ=255)
                                     tCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATa                                                      >  1:2075734/1‑91 (MQ=255)
                                       tttATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGa                                                    >  2:1337646/1‑91 (MQ=255)
                                        ttATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGaa                                                   >  1:2165309/1‑91 (MQ=255)
                                         tataATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAAt                                                  <  1:1696201/91‑1 (MQ=255)
                                          ataATTGTAGAACTATTTCTTATGCTGTTTAATTACTAAAAAAAAAAAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAAtt                                                 <  2:1683760/91‑1 (MQ=255)
                                           taATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTa                                                <  1:201596/91‑1 (MQ=255)
                                            aaTTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTag                                               >  1:772474/1‑91 (MQ=255)
                                            aaTTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTag                                               <  2:1480964/91‑1 (MQ=255)
                                              ttGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTagag                                             <  2:2522109/91‑1 (MQ=255)
                                               tGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAgaga                                            >  1:1141606/1‑91 (MQ=255)
                                               tGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAgaga                                            <  2:2161425/91‑1 (MQ=255)
                                                gTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGaa                                           >  2:1844645/1‑91 (MQ=255)
                                                gTAGCACTATTTCTTATGCTGTTTAATTACTATAATACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGaa                                           >  1:162467/1‑91 (MQ=255)
                                                  aGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAata                                         >  2:3101549/1‑91 (MQ=255)
                                                   gCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAatat                                        >  1:2854728/1‑91 (MQ=255)
                                                    cACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAAtata                                       >  2:2398993/1‑91 (MQ=255)
                                                     aCTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAg                                      >  2:529588/1‑91 (MQ=255)
                                                     aCTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAg                                      >  1:2811646/1‑91 (MQ=255)
                                                     aCTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAg                                      <  2:1218523/91‑1 (MQ=255)
                                                       tATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGaa                                    >  1:1430594/1‑91 (MQ=255)
                                                          ttCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAAc                                 >  2:680223/1‑91 (MQ=255)
                                                           tCTTATGCTGTTTAATTACAAAAATAAAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACa                                <  2:349106/91‑1 (MQ=255)
                                                             ttATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaa                              >  1:1369161/1‑91 (MQ=255)
                                                             ttATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaa                              >  1:1580257/1‑91 (MQ=255)
                                                              tATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaaa                             >  2:1738649/1‑91 (MQ=255)
                                                               aTGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaaaa                            <  2:844950/91‑1 (MQ=255)
                                                               aTGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaaaa                            >  1:3329669/1‑91 (MQ=255)
                                                                tGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAg                           >  2:1297920/1‑91 (MQ=255)
                                                                tGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAg                           >  2:418937/1‑91 (MQ=255)
                                                                tGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAg                           <  1:3167313/91‑1 (MQ=255)
                                                                 gCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAgg                          >  1:1869700/1‑91 (MQ=255)
                                                                  cTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGa                         <  1:729230/91‑1 (MQ=255)
                                                                      ttAATTACTGTACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAAttt                     <  1:3481719/91‑1 (MQ=255)
                                                                      ttAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAAttt                     >  2:2200067/1‑91 (MQ=255)
                                                                       tAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAACTAGAATTAGAGAATATAGAAAACAAAAAGGAATTtc                    >  2:2861284/1‑91 (MQ=255)
                                                                       tAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTtc                    >  1:516005/1‑91 (MQ=255)
                                                                       tAATTACTAGACTAGAAATAAAAAATATACATCTAAAGGTTAAATTAAGGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTtc                    <  1:3272745/91‑1 (MQ=255)
                                                                        aaTTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTtct                   >  1:2683188/1‑91 (MQ=255)
                                                                         aTTACTATAATACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAAAAAAAAAAAAAAAAAGGAATTtctc                  <  2:3312951/91‑1 (MQ=255)
                                                                          ttACTATACTACAAATGAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCa                 <  2:2789420/91‑1 (MQ=255)
                                                                           tactatactaCAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCaa                >  1:3063143/1‑91 (MQ=255)
                                                                            cctacactacaAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATATAATTAAAGAATATTAAAAACAAAAAGGAATTTCTCAAg               >  2:3382931/2‑91 (MQ=255)
                                                                            actatactaCAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAg               <  1:1094740/91‑1 (MQ=255)
                                                                            actatactaCAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAg               >  2:200853/1‑91 (MQ=255)
                                                                            actatactaCAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAg               >  2:1383296/1‑91 (MQ=255)
                                                                            actatactaCAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAg               <  1:966402/91‑1 (MQ=255)
                                                                              tatactaCGAATCAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATGAGCGAATATAGAAAACAAAAAGGAATTTCTCAAGaa             >  2:466444/1‑91 (MQ=255)
                                                                               atactaCAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGaaa            <  1:1954012/91‑1 (MQ=255)
                                                                                tactacAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAAc           >  1:248701/1‑91 (MQ=255)
                                                                                 actacAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACt          <  1:1244453/91‑1 (MQ=255)
                                                                                 actacAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACt          <  2:1018862/91‑1 (MQ=255)
                                                                                  ctacAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACtt         <  1:25227/91‑1 (MQ=255)
                                                                                   tacAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACttt        >  1:2216691/1‑91 (MQ=255)
                                                                                   tacAAATAAAAAATATACATCTAAAGATTAAATTAATGAAAAATAGAATTAGAGAAAAAAGAAAAAAAAAAGGAATTTCTCAAGAAACttt        <  2:1795650/91‑1 (MQ=255)
                                                                                     cAAATAAAAAATATACATCTAAAGGTCAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTaa      >  2:2070423/1‑91 (MQ=255)
                                                                                          |                                                                                          
ATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAA  >  CP000731/25573‑25753

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