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.98.4% 5331.4 / 62.5 1090coding (316/333 nt)USA300HOU_pUSA300HOUMR0031hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/18);  new base . (576/495);  total (576/514)
Fisher's exact test for biased strand distribution p-value = 1.11e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1279 total aligned reads displayed.

ATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAA  >  CP000731/25573‑25753
                                                                                          |                                                                                          
aTTGAAGGATCTTGGTTTTCTCAATATCCATAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATa                                                                                            <  2:3949558/91‑1 (MQ=255)
aTTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATa                                                                                            >  2:1516854/1‑91 (MQ=255)
  tGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaa                                                                                          >  1:630297/1‑91 (MQ=255)
  tGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaa                                                                                          >  2:1345779/1‑91 (MQ=255)
   gAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaa                                                                                         >  2:82880/1‑91 (MQ=255)
   gAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaa                                                                                         <  2:2045966/91‑1 (MQ=255)
   gAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaa                                                                                         <  2:5103338/91‑1 (MQ=255)
    aaGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaa                                                                                        >  2:4191693/1‑91 (MQ=255)
     aGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaa                                                                                       >  1:2755621/1‑91 (MQ=255)
     aGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaa                                                                                       >  1:196532/1‑91 (MQ=255)
     aGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaa                                                                                       >  1:783011/1‑91 (MQ=255)
       gATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaata                                                                                     >  1:3607802/1‑89 (MQ=255)
        aTCTTGGTTTTCTCAATATCCAGAACTTTTCTTTATAATTGTAGCACTACTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaatat                                                                                    <  2:5394600/91‑4 (MQ=255)
        aTCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCAATATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaatat                                                                                    <  2:5657206/91‑4 (MQ=255)
         tCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaatata                                                                                   <  1:5207008/91‑5 (MQ=255)
         tCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaatata                                                                                   <  2:1579526/91‑5 (MQ=255)
          cTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATAc                                                                                 <  1:5073973/91‑1 (MQ=255)
             ggTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATc                                                                              >  2:303972/1‑91 (MQ=255)
             ggTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATc                                                                              >  2:145767/1‑91 (MQ=255)
             ggTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATc                                                                              >  1:4440640/1‑91 (MQ=255)
               ttttCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTa                                                                            >  2:215048/1‑91 (MQ=255)
               ttttCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTa                                                                            <  2:3566494/91‑1 (MQ=255)
               ttttCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTa                                                                            <  2:4419353/91‑1 (MQ=255)
                tttCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTaa                                                                           >  2:4297818/1‑91 (MQ=255)
                 ttCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTaaa                                                                          >  2:1813510/1‑91 (MQ=255)
                  tctcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAg                                                                         <  1:3305329/91‑1 (MQ=255)
                  tctcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAg                                                                         <  2:4675182/91‑1 (MQ=255)
                   ctcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAgg                                                                        <  2:2280534/91‑1 (MQ=255)
                    tcAATATCCATAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTAGACTACAAATAAAAAATATACATCTAAAGGt                                                                       <  2:4618424/91‑1 (MQ=255)
                    tcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGt                                                                       <  2:2433091/91‑1 (MQ=255)
                    tcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGt                                                                       <  2:906007/91‑1 (MQ=255)
                    tcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGt                                                                       >  1:1129417/1‑91 (MQ=255)
                    tcAATATCAAGAAATTTTATTAATAATCGTAGAACTATTTCTTATGCTGTATAAATAATATACTACAAATAAAAAATATACATCTACAGGt                                                                       <  2:4828950/91‑1 (MQ=255)
                     cAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGtt                                                                      >  1:1238406/1‑91 (MQ=255)
                       atatCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTaa                                                                    <  2:1051056/91‑1 (MQ=255)
                        tatCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTaaa                                                                   <  1:1709970/91‑1 (MQ=255)
                        tatCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTaaa                                                                   >  1:92205/1‑91 (MQ=255)
                            cAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGTTTAAATTaa                                                               >  1:5308024/1‑91 (MQ=255)
                            cAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTaa                                                               <  1:2528504/91‑1 (MQ=255)
                            cAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTaa                                                               <  1:4968003/91‑1 (MQ=255)
                            cAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTaa                                                               <  2:4521673/91‑1 (MQ=255)
                             aGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAAt                                                              <  1:267839/91‑1 (MQ=255)
                             aGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAAt                                                              <  1:1801096/91‑1 (MQ=255)
                              gAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATg                                                             <  1:4340859/91‑1 (MQ=255)
                              gAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATg                                                             >  2:4759670/1‑91 (MQ=255)
                              aaaaTTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATg                                                             <  1:4651729/90‑1 (MQ=255)
                               aaaTTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGa                                                            <  2:3745007/91‑1 (MQ=255)
                               aaaTTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGa                                                            >  1:1003772/1‑91 (MQ=255)
                                aaTTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGaa                                                           >  2:1873222/1‑91 (MQ=255)
                                    ttCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAAt                                                       <  1:5544536/91‑1 (MQ=255)
                                     tCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATa                                                      >  1:3985869/1‑91 (MQ=255)
                                      cTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAg                                                     >  1:4748318/1‑91 (MQ=255)
                                      cTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAg                                                     <  1:861605/91‑1 (MQ=255)
                                      cTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAg                                                     >  1:4532727/1‑91 (MQ=255)
                                      cTTTATAATTGTAGCACTATTTCTTATGATGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAg                                                     >  2:5566647/1‑91 (MQ=255)
                                       tttATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGa                                                    <  1:2070578/91‑1 (MQ=255)
                                        ttATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGaa                                                   >  1:3802623/1‑91 (MQ=255)
                                           taATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTa                                                <  1:389296/91‑1 (MQ=255)
                                            aaTTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTag                                               <  2:704344/91‑1 (MQ=255)
                                                gTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGaa                                           <  1:2192925/91‑1 (MQ=255)
                                                gTAGCACTATTTCTTATGCTGTTTAATTACTAAACTAAAAATAAAAAATATACATCTAAAGGTAAAAAAAAAGAAAAATAGAATTAGAGaa                                           <  2:2668783/91‑1 (MQ=255)
                                                gTAGCACTATTTCTTATGCTGTATAATTACTATACTACAAATAAAAAATATACATCTCAAGGTTAAATTAATGACAAATAGAATAAGAGaa                                           >  2:3266288/1‑91 (MQ=255)
                                                  aGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAata                                         >  2:1680215/1‑91 (MQ=255)
                                                   gCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAatat                                        >  2:521168/1‑91 (MQ=255)
                                                   gCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAatat                                        >  1:699865/1‑91 (MQ=255)
                                                     aCTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAg                                      >  2:593878/1‑91 (MQ=255)
                                                     aCTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAg                                      >  2:3642545/1‑91 (MQ=255)
                                                     aCTATTTCTTATGCTAATAAATAACTATACTACAAAAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAAtataa                                      <  2:385408/91‑2 (MQ=255)
                                                      cTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAGATTAATGAAAAATAGAATTAGAGAATATAGa                                     <  1:5688685/91‑1 (MQ=255)
                                                        aTTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGaaa                                   <  1:1597178/91‑1 (MQ=255)
                                                           tCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACa                                <  1:2848990/91‑1 (MQ=255)
                                                           tCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACa                                >  2:3362435/1‑91 (MQ=255)
                                                             ttATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaa                              >  2:2363483/1‑91 (MQ=255)
                                                               aTGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATaaaaaaaaaaaaaa                            <  2:2787895/91‑1 (MQ=255)
                                                               aTGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaaaa                            >  2:4995658/1‑91 (MQ=255)
                                                                tGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAg                           <  1:235132/91‑1 (MQ=255)
                                                                tGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAg                           <  1:5367313/91‑1 (MQ=255)
                                                                 gCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAgg                          <  2:5493127/91‑1 (MQ=255)
                                                                 gCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAgg                          <  1:1449329/91‑1 (MQ=255)
                                                                 gCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAgg                          <  2:1259476/91‑1 (MQ=255)
                                                                  cTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGa                         <  2:1142852/91‑1 (MQ=255)
                                                                  cTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGa                         <  1:3215467/91‑1 (MQ=255)
                                                                  cTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAAAAAAAAAAAAAAAAAGGa                         <  2:2549125/91‑1 (MQ=255)
                                                                    gctgaaTTAATATACAACAAATAAAAGATATTGAGCTGAAGGTTAAATTAATGAAATATAGAAATAGAGAATATAGAAAACAAAAAGGAAt                       <  2:3452843/87‑1 (MQ=255)
                                                                       tAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTtc                    >  1:5413309/1‑91 (MQ=255)
                                                                         aTTACTAAAATAAAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTtctc                  <  2:2899378/91‑1 (MQ=255)
                                                                          ttACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCa                 <  2:3129191/91‑1 (MQ=255)
                                                                           tactatactaCAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCaa                <  1:2273119/91‑1 (MQ=255)
                                                                           tactatactaCAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCaa                >  1:3005866/1‑91 (MQ=255)
                                                                           tactatactaCAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCaa                <  2:1946068/91‑1 (MQ=255)
                                                                           tactatactaCAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCaa                >  1:3085687/1‑91 (MQ=255)
                                                                            actatactaCAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAg               >  1:1366259/1‑91 (MQ=255)
                                                                              tatactaCAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGaa             >  1:4126961/1‑91 (MQ=255)
                                                                              tatactaCAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGaa             >  2:403692/1‑91 (MQ=255)
                                                                                 actacAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACt          >  1:4263596/1‑91 (MQ=255)
                                                                                 actacAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACt          <  2:4100915/91‑1 (MQ=255)
                                                                                  ctacAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACtt         >  2:5333255/1‑91 (MQ=255)
                                                                                   tacAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACttt        <  1:2431403/91‑1 (MQ=255)
                                                                                     cAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAATCTTTaa      >  1:4865537/1‑91 (MQ=255)
                                                                                     cAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTaa      >  2:5772643/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: