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.3% 6939.2 / 31.5 1399coding (316/333 nt)USA300HOU_pUSA300HOUMR0031hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/10);  new base . (748/637);  total (751/648)
Fisher's exact test for biased strand distribution p-value = 4.44e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1602 total aligned reads displayed.

ATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAA  >  CP000731/25573‑25753
                                                                                          |                                                                                          
aTTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATa                                                                                            <  1:96481/91‑1 (MQ=255)
aTTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATa                                                                                            >  1:158477/1‑91 (MQ=255)
aTTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATa                                                                                            <  2:1082342/91‑1 (MQ=255)
aTTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATa                                                                                            <  2:2502009/91‑1 (MQ=255)
 ttGAAGGATCTTTGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaa                                                                                           <  2:1629625/91‑1 (MQ=255)
 ttGAAGGAACTTTGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTTTGCTGTTTAATTACTATAATACAAATaa                                                                                           <  2:6367945/91‑1 (MQ=255)
  tGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaa                                                                                          >  1:595891/1‑91 (MQ=255)
  tGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaa                                                                                          >  1:3798489/1‑91 (MQ=255)
   gAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaa                                                                                         >  1:5234905/1‑91 (MQ=255)
     aGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaa                                                                                       <  1:3553902/91‑1 (MQ=255)
     aGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaa                                                                                       <  1:6405805/91‑1 (MQ=255)
       gATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaata                                                                                     >  2:6182562/1‑89 (MQ=255)
        aTCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAGTTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaatat                                                                                    <  1:3687658/91‑4 (MQ=255)
         tCTTTGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaatata                                                                                   <  2:5135850/91‑5 (MQ=255)
         tCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATGACTATACTACAAATaaaaaatata                                                                                   >  1:6643288/1‑87 (MQ=255)
          cTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATAc                                                                                 >  2:2587029/1‑91 (MQ=255)
          cTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATAc                                                                                 <  1:1860358/91‑1 (MQ=255)
          cTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATAc                                                                                 >  2:3303733/1‑91 (MQ=255)
          cTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATAc                                                                                 >  1:3438833/1‑91 (MQ=255)
            tGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACAt                                                                               <  1:5140454/91‑1 (MQ=255)
            tGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACAt                                                                               >  2:1197451/1‑91 (MQ=255)
            tGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACAt                                                                               >  1:6312961/1‑91 (MQ=255)
            tGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACAt                                                                               <  2:706722/91‑1 (MQ=255)
             ggTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATc                                                                              >  1:2389336/1‑91 (MQ=255)
                tttCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTaa                                                                           >  2:82911/1‑91 (MQ=255)
                 ttCTCAATATCCAAAAATTTTCTTTATAATGGTGGCATTATTTCTTTGGCTGTTTAATTATTAAATTAAAAAAAAAAAATATACTTCTaaa                                                                          <  2:1956133/91‑1 (MQ=255)
                   ttcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAgg                                                                        <  1:6770899/90‑1 (MQ=255)
                   ctcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAgg                                                                        <  2:2714136/91‑1 (MQ=255)
                   ctcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAgg                                                                        <  1:880875/91‑1 (MQ=255)
                   ctcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAgg                                                                        <  2:5482105/91‑1 (MQ=255)
                    tcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGt                                                                       >  1:168910/1‑91 (MQ=255)
                     cAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGtt                                                                      >  1:4249167/1‑91 (MQ=255)
                      aaTATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTa                                                                     <  2:3518089/91‑1 (MQ=255)
                           ccAGAAATTTTCTTTATAATTGGAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTa                                                                >  2:5340409/1‑91 (MQ=255)
                            cAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTaa                                                               <  1:3335840/91‑1 (MQ=255)
                            cAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTaa                                                               >  1:6137527/1‑91 (MQ=255)
                             aGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAAt                                                              >  1:276667/1‑91 (MQ=255)
                              gAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATg                                                             <  1:2205997/91‑1 (MQ=255)
                               aaaTTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGa                                                            >  2:5612816/1‑91 (MQ=255)
                               aaaTTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGa                                                            >  2:4160967/1‑91 (MQ=255)
                                aaTTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGaa                                                           <  1:5521274/91‑1 (MQ=255)
                                  ttttCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGaaaa                                                         <  2:3823451/91‑1 (MQ=255)
                                  ttttCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGaaaa                                                         <  2:4284335/91‑1 (MQ=255)
                                     tCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATa                                                      >  1:1426618/1‑91 (MQ=255)
                                      cTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAg                                                     >  2:5791149/1‑91 (MQ=255)
                                      cTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAg                                                     >  2:3946345/1‑91 (MQ=255)
                                       tttATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGa                                                    <  2:1440584/91‑1 (MQ=255)
                                        ttATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGaa                                                   >  2:4653694/1‑91 (MQ=255)
                                         ttgaaTTGTAGCACTAAGTCGTATGCTGTTTAATTACTATACTAAAAAAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAAt                                                  <  2:1724289/88‑1 (MQ=255)
                                          ataATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAAtt                                                 <  1:2490501/91‑1 (MQ=255)
                                              ttGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTagag                                             >  1:1006491/1‑91 (MQ=255)
                                              ttGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTagag                                             <  1:1108083/91‑1 (MQ=255)
                                              ttGTAGCACTATTTCTTATGCTGTTTAATTACTAAACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTagag                                             <  2:4874775/91‑1 (MQ=255)
                                               tGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAgaga                                            <  1:3156708/91‑1 (MQ=255)
                                                gTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGa                                                    >  2:6622472/1‑82 (MQ=255)
                                                gTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGaa                                           <  1:152876/91‑1 (MQ=255)
                                                gTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGaa                                           <  2:5927686/91‑1 (MQ=255)
                                                gTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGaa                                           <  2:2076918/91‑1 (MQ=255)
                                                gTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGaa                                           <  1:5364978/91‑1 (MQ=255)
                                                gTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGaa                                           >  2:3617484/1‑91 (MQ=255)
                                                gTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGaa                                           <  1:4107876/91‑1 (MQ=255)
                                                 tAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAat                                          <  1:4938672/91‑1 (MQ=255)
                                                  aGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAata                                         >  2:3434618/1‑91 (MQ=255)
                                                   gCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAatat                                        <  1:5041912/91‑1 (MQ=255)
                                                    cACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAAtata                                       <  1:4381122/91‑1 (MQ=255)
                                                    cACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAAtata                                       >  2:2415018/1‑91 (MQ=255)
                                                    cACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAAtata                                       >  2:1527634/1‑91 (MQ=255)
                                                    aaCTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAAtata                                       <  2:478916/90‑1 (MQ=255)
                                                     aCTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTCATGAAAAATAGAATTAGAGCATATAg                                      >  2:4412895/1‑91 (MQ=255)
                                                     aCTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAg                                      >  2:4995278/1‑91 (MQ=255)
                                                     aCTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAg                                      >  1:2911955/1‑91 (MQ=255)
                                                      cTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGa                                     >  2:6078416/1‑91 (MQ=255)
                                                      cTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGa                                     >  1:2653889/1‑91 (MQ=255)
                                                       tATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGaa                                    >  1:4675107/1‑91 (MQ=255)
                                                         tttCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGaaaa                                  >  2:5984218/1‑91 (MQ=255)
                                                            cTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaa                               <  1:6929993/91‑1 (MQ=255)
                                                             ttATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaa                              >  2:5206965/1‑91 (MQ=255)
                                                             ttATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaa                              >  2:2314816/1‑91 (MQ=255)
                                                             ttATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaa                              <  1:205973/91‑1 (MQ=255)
                                                              tATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaaa                             <  2:4022748/91‑1 (MQ=255)
                                                                tGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAg                           >  1:5615630/1‑91 (MQ=255)
                                                                tGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAg                           >  2:3055408/1‑91 (MQ=255)
                                                                  cTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGa                         >  2:2164508/1‑91 (MQ=255)
                                                                    gTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATCGAATTAGAGAATATAGAAAACAAAAAGGAAt                       >  1:6520633/1‑91 (MQ=255)
                                                                          ttACTATACTCCAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAACAGGAATTTCTCa                 >  2:4518247/1‑91 (MQ=255)
                                                                          ttACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCa                 <  2:3151789/91‑1 (MQ=255)
                                                                          ttACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCa                 <  1:4811494/91‑1 (MQ=255)
                                                                          ttACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCa                 <  2:6481720/91‑1 (MQ=255)
                                                                          ttACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCa                 <  1:3925408/91‑1 (MQ=255)
                                                                          ttACTATACTACAAATAAAAAATATACAACTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCa                 <  1:1971675/91‑1 (MQ=255)
                                                                           tactatactaCAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCaa                >  1:4564239/1‑91 (MQ=255)
                                                                            actatactaCAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAg               <  1:755087/91‑1 (MQ=255)
                                                                            actatactaCAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAAAAAACAAAAAGGAATTTCTCAAg               <  2:6840842/91‑1 (MQ=255)
                                                                             ctatactaCAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGa              <  2:1801313/91‑1 (MQ=255)
                                                                                 actacAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACt          <  1:5904238/91‑1 (MQ=255)
                                                                                 actacAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACt          >  1:1299846/1‑91 (MQ=255)
                                                                                 actacAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACt          >  1:3030647/1‑91 (MQ=255)
                                                                                 actacAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACt          >  2:6743189/1‑91 (MQ=255)
                                                                                   tacAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACttt        <  2:2875632/91‑1 (MQ=255)
                                                                                    acAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTa       <  1:575085/91‑1 (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: