Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000731 11,100 C→A 5.9% intergenic (+105/‑159) blaI → / → USA300HOU_pUSA300HOUMR0014 beta‑lactamase regulator BlaI/recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073111,1000CA 5.9% 1617.9 / 63.7 454intergenic (+105/‑159)blaI/USA300HOU_pUSA300HOUMR0014beta‑lactamase regulator BlaI/recombinase
Reads supporting (aligned to +/- strand):  ref base C (222/205);  new base A (5/22);  total (227/227)
Fisher's exact test for biased strand distribution p-value = 1.09e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 463 total aligned reads displayed.

AAACATAATACATATGTAAACCATACTTGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAATTATGTGTTCAAAAACCCAACAAAAAA  >  CP000731/11011‑11190
                                                                                         |                                                                                          
aaaCATAATACATATGTAAACCATACTTGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTcc                                                                                           <  1:760642/91‑1 (MQ=255)
  gtcttttgtgaacatagattttttgttgggtttttgaacacataattctatgtaatacttacttatttagacaaaaaacaAGAGTGTCCAg                                                                                         >  1:97767‑M1/81‑91 (MQ=255)
  aCATAATACATATGTAAACCATACTTGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAg                                                                                         >  1:1021086/1‑91 (MQ=255)
      aataCATATGTAAACCATACTTGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGaaaa                                                                                     >  1:715352/1‑91 (MQ=255)
       ataCATATGTAAACCATACTTGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAAt                                                                                    >  1:2161094/1‑91 (MQ=255)
        taCATATGTAAACCATACTTGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATa                                                                                   >  1:1637184/1‑91 (MQ=255)
         aCATATGTAAACCATACTTGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATaa                                                                                  >  1:562375/1‑91 (MQ=255)
         aCATATGTAAACCATACTTGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATaa                                                                                  >  1:1236705/1‑91 (MQ=255)
         aCATATGTAAACCATACTTGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATaa                                                                                  <  1:2529976/91‑1 (MQ=255)
         aCATATGTAAACCATACTTGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATaa                                                                                  <  1:1405679/91‑1 (MQ=255)
          cATATGTAAACCATACTTGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAAc                                                                                 >  1:3570999/1‑91 (MQ=255)
            tatGTAAACCATACTTGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCa                                                                               <  1:515608/91‑1 (MQ=255)
            tatGTAAACCATACTTGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCa                                                                               <  1:3115521/91‑1 (MQ=255)
             atGTAAACCATACTTGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCAt                                                                              <  1:1143129/91‑1 (MQ=255)
             acatagattttttgttgggtttttgaacacataattctatgtaatacttacttatttagacaaaaaacaAGAGTGTACAGAGAATAACCAt                                                                              >  1:2089016‑M1/70‑91 (MQ=255)
              tGTAAACCATACTTGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCAtt                                                                             <  1:2557981/91‑1 (MQ=255)
              catagattttttgttgggtttttgaacacataattctatgtaatacttacttatttagacaaaaaacaAGAGTGTCCAGAAAATAACCAtt                                                                             >  1:3074186‑M1/69‑91 (MQ=255)
              catagattttttgttgggtttttgaacacataattctatgtaatacttacttatttagacaaaaaacaAGAGTGTCCAGAAAATAACCAtt                                                                             >  1:2900701‑M1/69‑91 (MQ=255)
              catagattttttgttgggtttttgaacacataattctatgtaatacttacttatttagacaaaaaacaAGAGTGTCCAGAAAATAACCAtt                                                                             >  1:1309685‑M1/69‑91 (MQ=255)
              catagattttttgttgggtttttgaacacataattctatgtaatacttacttatttagacaaaaaacaAGAGTGTCCAGAAAATAACCAtt                                                                             >  1:4060234‑M1/69‑91 (MQ=255)
               gTAAACCATACTTGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCAttt                                                                            <  1:1745642/91‑1 (MQ=255)
               atagattttttgttgggtttttgaacacataattctatgtaatacttacttatttagacaaaaaacaAGAGTGTCCAGAAAATAACCAttt                                                                            >  1:321638‑M1/68‑91 (MQ=255)
                tAAACCATACTTGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCAtttt                                                                           <  1:1005147/91‑1 (MQ=255)
                  gattttttgttgggtttttgaacacataattctatgtaatacttacttatttagacaaaaaacaAGAGTGTCCAGAAAATAACCATTTTCt                                                                         >  1:2587493‑M1/65‑91 (MQ=255)
                   attttttgttgggtttttgaacacataattctatgtaatacttacttatttagacaaaaaacaAGAGTGTCCAGAAAATAACCATTTTCTg                                                                        >  1:2379167‑M1/64‑91 (MQ=255)
                   aCCATACTTGACTCTTATTAAATTAGTTTTGGTAATGACATTTAGAATATATTCTTATAGATGAGAGTGTCCAGAAAATAACCATTTTCTg                                                                        <  1:3310344/91‑1 (MQ=255)
                    ccATACTTGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGt                                                                       >  1:1475567/1‑91 (MQ=255)
                     tttttgttgggtttttgaacacataattctatgtaatacttacttatttagacaaaaaacaAGAGTGTCCAGAAAATAACCATTTTCTGTa                                                                      >  1:842969‑M1/62‑91 (MQ=255)
                      aTACTTGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTac                                                                     >  1:2619329/1‑91 (MQ=255)
                      aTACTTGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTac                                                                     >  1:3884014/1‑91 (MQ=255)
                       tACTTGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTaca                                                                    >  1:1981418/1‑91 (MQ=255)
                       tACTTGACTCTTATTAAATTAGTTTTGTGAATGACATTTAAAATATATTATTAAAGATGAGAGTGTACAGAAAATGGTTATTTTCTGGaca                                                                    <  1:274183/91‑1 (MQ=25)
                         cTTGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACAct                                                                  >  1:3656256/1‑91 (MQ=255)
                          ttGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACActc                                                                 <  1:2727160/91‑1 (MQ=255)
                           tGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACActct                                                                <  1:2196312/91‑1 (MQ=255)
                            gACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTACAGAAAATGGTTATTTTCTGGACActctc                                                               <  1:2870768/91‑2 (MQ=21)
                             aCTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTg                                                              <  1:2685340/91‑1 (MQ=255)
                             aCTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTACAGAAAATGGTTATTTTCTGGACActctca                                                              <  1:1689181/91‑3 (MQ=21)
                             aCTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTACAGAAAATGGTTATTTTCTGGACActctca                                                              <  1:1927070/91‑3 (MQ=21)
                              ctctTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGt                                                             >  1:2856387/1‑91 (MQ=255)
                              ctctTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGt                                                             >  1:2804334/1‑91 (MQ=255)
                                  tattaAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGttttt                                                         >  1:1779945/1‑91 (MQ=255)
                                      aaaTTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTc                                                     <  1:3087467/91‑1 (MQ=255)
                                      aaaTTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTc                                                     <  1:3607707/91‑1 (MQ=255)
                                        aTTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTa                                                   <  1:2461063/91‑1 (MQ=255)
                                        aTTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTa                                                   >  1:3710838/1‑91 (MQ=255)
                                         ttAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTaa                                                  >  1:1871994/1‑91 (MQ=255)
                                          tAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTaaa                                                 >  1:639607/1‑91 (MQ=255)
                                           aGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAAt                                                >  1:800790/1‑91 (MQ=255)
                                             ttttGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAAtaa                                              <  1:812120/91‑1 (MQ=255)
                                             ttttGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAAtaa                                              <  1:1279661/91‑1 (MQ=255)
                                             ttttGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAAtaa                                              <  1:1765533/91‑1 (MQ=255)
                                               ttGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAAtaagt                                            <  1:3930893/91‑1 (MQ=255)
                                               ttGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAAtaagt                                            <  1:22990/91‑1 (MQ=255)
                                                tGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAAtaagta                                           >  1:380297/1‑91 (MQ=255)
                                                tGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAAtaagta                                           >  1:1090014/1‑91 (MQ=255)
                                                tGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAAtaagta                                           >  1:3161496/1‑91 (MQ=255)
                                                 gTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAAtaagtaa                                          >  1:915313/1‑91 (MQ=255)
                                                  ttAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAAtaagtaag                                         >  1:3831075/1‑91 (MQ=255)
                                                     aTGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTAt                                      <  1:487004/91‑1 (MQ=255)
                                                      tGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTAtt                                     >  1:1906878/1‑91 (MQ=255)
                                                       gACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTa                                    <  1:2258783/91‑1 (MQ=255)
                                                        aCATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTGTTAc                                   >  1:3461530/1‑91 (MQ=255)
                                                         tacttacttatttagacaaaaaacaAGAGTGTACAGAAAATGGTTATTTTCTGGACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACa                                  <  1:3659643‑M1/66‑1 (MQ=255)
                                                         cATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACa                                  >  1:695756/1‑91 (MQ=255)
                                                          aTTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACAt                                 >  1:86972/1‑91 (MQ=255)
                                                           tttAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATa                                <  1:2764552/91‑1 (MQ=255)
                                                            ttAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAg                               <  1:1539149/91‑1 (MQ=255)
                                                            ttAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAg                               <  1:3046386/91‑1 (MQ=255)
                                                             tAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGa                              <  1:3901914/91‑1 (MQ=255)
                                                              aaaaTATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGaa                             >  1:3040442/1‑91 (MQ=255)
                                                              aaaaTATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGaa                             >  1:375744/1‑91 (MQ=255)
                                                               aaaTATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAAt                            <  1:470976/91‑1 (MQ=255)
                                                               aaaTATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAAt                            >  1:2819232/1‑91 (MQ=255)
                                                               aaaTATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAAt                            >  1:3260252/1‑91 (MQ=255)
                                                                aaTATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAAtt                           <  1:3361141/91‑1 (MQ=255)
                                                                aaTATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAAtt                           >  1:2420017/1‑91 (MQ=255)
                                                                 atatatTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAATTa                          <  1:2445432/91‑1 (MQ=255)
                                                                    tattatTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAATTAtgt                       >  1:3410329/1‑91 (MQ=255)
                                                                    tattatTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAATTAtgt                       <  1:1445136/91‑1 (MQ=255)
                                                                    tattatTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAATTAtgt                       >  1:1493109/1‑91 (MQ=255)
                                                                      ttattaAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAATTATgtgt                     <  1:2016991/91‑1 (MQ=255)
                                                                        attaAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAATTATGTGTTc                   <  1:3507478/91‑1 (MQ=255)
                                                                         caaaaaacaAGAGTGTACAGAAAATGGTTATTTTCTGGACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAATTATGTGTTCa                  <  1:2991669‑M1/82‑1 (MQ=255)
                                                                              gATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAATTATGTGTTCAAAAAc             <  1:2328305/91‑1 (MQ=255)
                                                                               aTGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAATTATGTGTTCAAAAAcc            <  1:2343751/91‑1 (MQ=255)
                                                                                tGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAATTATGTGTTCAAAAAccc           <  1:1116962/91‑1 (MQ=255)
                                                                                caagagTGTACAGAAAAGGGTTATTTACTGGACACTCTTGTTTTTCGTCGAAAAAAGAAAGTATTACATAGAAGTATGTGTTCAAAAAccc           >  1:223632/3‑91 (MQ=12)
                                                                                 gagaGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAATTATGTGTTCAAAAACCca          >  1:733642/1‑91 (MQ=255)
                                                                                  AGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAATTATGTGTTCAAAAACCcaa         <  1:126847‑M1/91‑1 (MQ=255)
                                                                                  AGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAATTATGTGTTCAAAAACCcaa         >  1:1579055‑M1/1‑91 (MQ=255)
                                                                                   GAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAATTATGTGTTCAAAAACCcaac        >  1:403346‑M1/1‑91 (MQ=255)
                                                                                   GAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAATTATGTGTTCAAAAACCcaac        >  1:3141575‑M1/1‑91 (MQ=255)
                                                                                    AGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAATTATGTGTTCAAAAACCcaaca       >  1:1204667‑M1/1‑91 (MQ=255)
                                                                                     GTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAGGTAAGTATTACATAGAATTATGTGTTCAAAAACCcaacaa      <  1:1677409‑M1/91‑1 (MQ=255)
                                                                                     GTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAATTATGTGTTCAAAAACCcaacaa      >  1:414005‑M1/1‑91 (MQ=255)
                                                                                     GTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAATTATGTGTTCAAAAACCcaacaa      <  1:1177510‑M1/91‑1 (MQ=255)
                                                                                     GTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAATTATGTGTTCAAAAACCcaacaa      <  1:2787250‑M1/91‑1 (MQ=255)
                                                                                       GTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAATTATGTGTTCGAAAACCCAACaaaa    >  1:2925137‑M1/1‑91 (MQ=255)
                                                                                        TCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAATTATGTGTTCAAAAACCCAACaaaaa   >  1:4028956‑M1/1‑91 (MQ=255)
                                                                                         |                                                                                          
AAACATAATACATATGTAAACCATACTTGACTCTTATTAAATTAGTTTTGTTAATGACATTTAAAATATATTATTAAAGATGAGAGTGTCCAGAAAATAACCATTTTCTGTACACTCTTGTTTTTTGTCTAAATAAGTAAGTATTACATAGAATTATGTGTTCAAAAACCCAACAAAAAA  >  CP000731/11011‑11190

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