Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,233,544 T→C 6.4% E174G (GAA→GGA)  USA300HOU_2105 ← HemK family methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,233,5440TC 6.4% 171.1 / 2.1 64E174G (GAA→GGA) USA300HOU_2105HemK family methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (34/24);  new base C (2/2);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.67e-01

AGCATATCCATGGTTATCTGCAAACAATGCCTGATGTGGTTCGAATCTCGTAACCGTTGGAGACATCGTAACCATATCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCAAGTGAAATATCAGTAGCAATAACATTTAAATCCGGCTTTTCACATTTC  >  CP000730/2233414‑2233684
                                                                                                                                  |                                                                                                                                              
aGCATATCCATGGTTATCTGCAAACAATGCCTGATGTGGTTCGAATCTCGTAACCGTTGGAGACATCGTAACCATATCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCAttaattaa                                                                                                                                      <  2:410632/141‑1 (MQ=255)
      tCCATGGTTATCTGCAAACAATGCCTGATGTGGTTCGAATCTCGTAACCGTTGGAGACATCGTAACCATATCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCTTTCAACTTGATAGCTTCATTAATTAAAGACtt                                                                                                                                >  2:571093/1‑141 (MQ=255)
        cATGGTTATCTGCAAACAATGCCTGATGTGGTTCGAATCTCGTAACCGTTGGAGACATCGTAACCATATCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTAAACTTGATACCTTCATTAATTAAGGGCTTTa                                                                                                                              >  1:13040/1‑141 (MQ=255)
           ggTTTTCTGCAAACAATGCCTGATGTGGTTCGAATCTCGTAACCGTTGGAGACATCGTAACCATATCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATg                                                                                                                           <  2:452598/141‑1 (MQ=255)
            gTTATCTGCAAACAATGCCTGATGTGGTTCGAATCTCGTAACCGTTGGAGACATCGTAACCATATCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGc                                                                                                                          >  1:132274/1‑141 (MQ=255)
                  tGCAAACAATGCCTGATGTGGTTCGAATCTCGTAACCGTTGGAGACATCGTAACCATATCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATcccc                                                                                                                    >  1:166672/1‑141 (MQ=255)
                   gCAAATAATGCCTGATGTGGTTCGAATCTCGTAACCGTTGGAGACATCGTAACCATATCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCt                                                                                                                   <  2:598100/141‑1 (MQ=255)
                    cAAACAATGCCTGATGTGGTTCGAATCTCGTAACCGTTGGAGACATCGTAACCATATCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTCAACCTGATACCTTCATTAATTAAGGACTTTAATGCATCCCCTg                                                                                                                  >  2:197654/1‑141 (MQ=255)
                      aaCAATGCCTGATGTGGTTCGAATCTCGTAACCGTTGGAGACATCGTAACCATATCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGtt                                                                                                                >  1:736392/1‑141 (MQ=255)
                            gCCTGATGTGGTTCGAATCTCGCAACCTCTGGAGACATCGTAACCATATCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAAt                                                                                                         <  1:478925/141‑1 (MQ=255)
                                    tgGTTCGAATCTCGTAACCGTTGGAGACATCGTAACCATATCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTCATGCATCCCCTGTTAAAAATTGTATTtg                                                                                                <  2:543971/141‑1 (MQ=255)
                                    tgGTCCGAATCCCGTAACCGTTGGAGACATCGTAACCATATCGTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTtg                                                                                                <  1:652316/141‑1 (MQ=255)
                                       ttCGAATCTCGTAACCGTTGGAGACATCGTAACCATATCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTtgtgt                                                                                             >  1:702621/1‑140 (MQ=255)
                                       ttCGAATCTCGTAACCGATGGAGACATCGTAACCATATCTTTTTCATCTATATATGGTAGGTTAGATATCAAGCCGTTCAAATTGATACATTCATTAATTATCGGCTTTAATACATACCCTGTTAAAAATCGTAtgtaggg                                                                                              >  2:458571/1‑135 (MQ=255)
                                             tctcGTAACCGTTGGAGACATCGTAACCATATCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATg                                                                                       <  2:722971/141‑1 (MQ=255)
                                              ctcGTAACCGTTGGAGACATCGTAACCATATCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGc                                                                                      <  1:87974/141‑1 (MQ=255)
                                                 gTAACCGTTGGAGACATCGTAACCATATCTTTTTCATCTCTATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTtc                                                                                   >  2:671820/1‑141 (MQ=255)
                                                 gTAACCGTTGGAGACATCGTAACCATATCTTTTTCATCTATATATAGTGGATTAGATATCAAGCCGTTCACCTTGATACCTTCATTACTTAAGGACTTTAATGAATCCCCCGTTAAAAAGTGTATTTGtgattagtgcgcc                                                                                   >  2:156089/1‑133 (MQ=255)
                                                  tAACCGGTGGAGACATCGTAACCATATCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTtct                                                                                  <  1:300052/141‑1 (MQ=255)
                                                   aaCCGTTGGAGACATCGTAACCATATCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTtctc                                                                                 >  2:19121/1‑141 (MQ=255)
                                                    aCCGTTGGAGACATCGTAACCATATCTTTTTTTTCTATTTATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTTTATTTGTGATTGATGCTTCTCa                                                                                >  2:242535/1‑141 (MQ=255)
                                                    aCCGTTGGAGACATCGTAACCATATCTTTTTCAGCTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCa                                                                                <  2:258713/141‑1 (MQ=255)
                                                     ccGTTGGAGACATCGTAACCATATCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTCCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAg                                                                               >  2:9864/1‑141 (MQ=255)
                                                     ccGTTGGAGACATCGTAACCATATCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTCCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAg                                                                               >  2:9859/1‑141 (MQ=255)
                                                                 tCGTAACCATATCTTTTTCATTTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTCCGAg                                                                   <  2:7749/141‑1 (MQ=255)
                                                                 tCGTAACATTATCTTTTTCATCTATGTAGGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAg                                                                   <  2:119098/141‑1 (MQ=255)
                                                                     aaCCATATCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCat                                                               >  2:211155/1‑141 (MQ=255)
                                                                       ccATATCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCatat                                                             >  2:603730/1‑141 (MQ=255)
                                                                       ccATATCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCatat                                                             >  1:278503/1‑141 (MQ=255)
                                                                         atatCTTTTTCATCTATTTATGGTGGATTAGTTAGCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCTCCTGTTAAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATAtt                                                            <  2:721442/141‑1 (MQ=255)
                                                                         atatCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTc                                                           >  2:275415/1‑141 (MQ=255)
                                                                         atatCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTc                                                           <  2:713831/141‑1 (MQ=255)
                                                                              tttttCATCTATATATGGTGGATTAGATATCAAGCCGCTCAACCTCACACTGTCAGTAATTATGGGCTTTAATGCATCCCCTGTTAAAAATGGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGc                                                      <  2:349783/141‑1 (MQ=255)
                                                                               ttttCATCTATATATGGTGGATTAGAGATCAAGCCGTTCAACTTGATACCTGCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCt                                                     <  1:290550/141‑1 (MQ=255)
                                                                                     tCTATATATGGTGGATTAGTTATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCACATTCATTGCTTCAAGt                                               <  1:381444/141‑1 (MQ=255)
                                                                                         tatatGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTATTGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTTCTTCAAGTGaaa                                           >  2:515254/1‑141 (MQ=255)
                                                                                             tggtggATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCAAGTGAAATATc                                       >  1:455633/1‑141 (MQ=255)
                                                                                             tggtggATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCAAGTGAAATATc                                       >  2:134314/1‑141 (MQ=255)
                                                                                                   tttAGGTATCAAGCCGTTCAACCTGATCCCTTCATTAATGTAGGGCTTTCAGGCATCCTCTGTTAAAAATTGTATTTGTGATTGGTGCTTCTCAGCATTATTACGATCCATATTCATTGCTTCAAGTGAAATATCAGTAGc                                 <  2:267696/140‑1 (MQ=255)
                                                                                                   aTTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCGCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGATCCATATTTATTGCTTCAACTGAAATATCAGTAGc                                 >  1:465436/1‑141 (MQ=255)
                                                                                                   aTTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCAAGTGAAATATCAGTAGc                                 <  1:11224/141‑1 (MQ=255)
                                                                                                   aTTAGATATCAAGCCGGTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCAAGTGAAATATCAGTAGc                                 <  1:11156/141‑1 (MQ=255)
                                                                                                     tagccatCAAGCCGTTCAACTTGATACCTACATTAATTAAGTGCTTTAATTCATCCCCTTTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCGTTGCTTCAAGTGAAATATCAGTAGCaa                               >  2:640481/6‑141 (MQ=255)
                                                                                                       gATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCAAGTGAAATATCAGTAGCAATa                             >  2:313143/1‑141 (MQ=255)
                                                                                                        atatCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCAAGTGGAATATCAGTAGCAATaa                            >  2:730577/1‑141 (MQ=255)
                                                                                                        atatCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCAAGTGAAATATCAGTAGCAATaa                            >  1:308041/1‑141 (MQ=255)
                                                                                                         tatCAAGCCGTTCAACTTGATACCTCCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCAAGTGAAATATCAGTAGCAATAAc                           <  1:9859/141‑1 (MQ=255)
                                                                                                         tatCAAGCCGTTCAACTTGATACCTCCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCAAGTGAAATATCAGTAGCAATAAc                           <  1:9864/141‑1 (MQ=255)
                                                                                                           tCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCAAGTGAAATATCAGTAGCAAAAACAt                         >  2:32055/1‑141 (MQ=255)
                                                                                                                ccGTTCAACTTGATACGTTCATTAATTAAGGGCTTTAAAGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATAGCTTCAAGTGAAATATCAGTAGCAATAACATTTaaa                    >  1:427625/1‑141 (MQ=255)
                                                                                                                ccGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCCTTATTACGACCCATATTCATTGCTCCAACTGAAATAACAGTAGCAATAACATTTaac                    >  2:678587/1‑140 (MQ=255)
                                                                                                                ccGTTCAACTTGATACATTCATTAATTAAGGGCTTTAAGGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATAGCTTCAAGTGAAATATCAGTAGCAATAACATTTaaa                    >  1:427644/1‑141 (MQ=255)
                                                                                                                    tCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCAAGTGAAATATCAGTAGCAATAACATTTAAATCCg                >  1:649745/1‑141 (MQ=255)
                                                                                                                    tCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCAAGTGAAATATCAGTAGCAATAACATTTAAATCCg                >  2:661888/1‑141 (MQ=255)
                                                                                                                        cTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCACGTAAAATATCAGTAGCAATAACGTTTAAATCCGGcct            >  1:437965/1‑139 (MQ=255)
                                                                                                                        cTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCAAGTGAAATATCAGTAGCAATAACATTTAAATCCGGCtt            >  1:738139/1‑141 (MQ=255)
                                                                                                                        cTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATATGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCAAGTGAAATATCAGTAGCAGTAACATTTAAATCCGGCtt            >  1:511391/1‑141 (MQ=255)
                                                                                                                        cTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATATGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCAAGTGAAATATCAGTAGCAATAACATTTAAATCCGGCtt            >  1:511416/1‑141 (MQ=255)
                                                                                                                         ttGATACCTTTCTTAATTAAGGGCTTTAATGCATCCAATCTTCAACATTGTCTCGTATGATGTTCGCCTGCCGCATTAATAAGTGTATTATTCATTGCTTCAAGTGAAATATCAGTAGCAATAACATTTAAATCCGGCttt           <  2:540673/141‑1 (MQ=255)
                                                                                                                         ttGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCAAGTGAAATATCAGTAGCAATAACATTTAAATCCGGCttt           >  1:321085/1‑141 (MQ=255)
                                                                                                                         ttGATACCTTCATTAATGAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCAAGTGAAATATCAGTAGCAATAACATTTAAATCCGGCttt           <  2:34583/141‑1 (MQ=255)
                                                                                                                          tGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCAAGTGAAATATCAGTAGCAATAACATTTAAATCCGGCtttt          <  1:19121/141‑1 (MQ=255)
                                                                                                                          atacaCCTTCATTAATTAAGTGCTTGAATGCAAACAAAAATAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCAAGTGAAATATCAGTAGCAATAACATTTAAATCCGGCtttt          <  2:619150/137‑1 (MQ=255)
                                                                                                                            aTAACTTCATCAATTGAGTGCTTTAATGCATCCCCCGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCAAGTGAAATATCAGTAGCAAGAACATTTAAATCCGGCTTTTca        <  2:574193/141‑1 (MQ=255)
                                                                                                                                 ttCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCAAGTGAAATATCAGTAGCAATAACATTTAAATCCGGCTTTTCACCttt   >  2:457215/1‑141 (MQ=255)
                                                                                                                                  tCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTCTTTGTGATTGATGCTTCTCAGCATTATTACGAGCCACCTTCATTGCTTCAAGTGAAATATCAGTAGCAATAACATTTAAATCCGGCTTTTCACATTTc  <  2:162501/141‑1 (MQ=255)
                                                                                                                                  tCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCAAGTGAAATATCAGTAGCAATAACATTTAAATCCAGCTATTCACATTTc  <  2:520017/141‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                              
AGCATATCCATGGTTATCTGCAAACAATGCCTGATGTGGTTCGAATCTCGTAACCGTTGGAGACATCGTAACCATATCTTTTTCATCTATATATGGTGGATTAGATATCAAGCCGTTCAACTTGATACCTTCATTAATTAAGGGCTTTAATGCATCCCCTGTTAAAAATTGTATTTGTGATTGATGCTTCTCAGCATTATTACGAGCCATATTCATTGCTTCAAGTGAAATATCAGTAGCAATAACATTTAAATCCGGCTTTTCACATTTC  >  CP000730/2233414‑2233684

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