Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,521,608 G→A 6.6% intergenic (+548/+18) pbp2 → / ← USA300HOU_1389 penicillin‑binding protein 2/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,521,6080GA 6.6% 217.0 / 2.8 71intergenic (+548/+18)pbp2/USA300HOU_1389penicillin‑binding protein 2/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (37/29);  new base A (2/3);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

AGATATCTTAAAGCATATCTCAAGATATCCATTATGATGTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTAAAGGTCTTCATTTGTCGATAGCCCATAAAATGATATTTACGTTCAATCATAT  >  CP000730/1521470‑1521741
                                                                                                                                          |                                                                                                                                     
aGATATCTTAAAGCATATCTCAAGATATCCATTATGATGTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGa                                                                                                                                     >  2:793246/1‑141 (MQ=255)
 gATATCTTAAAGCATATCTCAAGATATCCATTATGATGTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGaa                                                                                                                                    <  1:434990/141‑1 (MQ=255)
  atatCTTAAAGCATATCTCAAGATATCCATTATGATGTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGaaa                                                                                                                                   >  1:958905/1‑141 (MQ=255)
           aGCATATCTCAAGATATCCATTATGATGTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATc                                                                                                                          >  2:270887/1‑141 (MQ=255)
            gCATATCTCAAGATATCCATTATGATGTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCt                                                                                                                         >  1:416126/1‑141 (MQ=255)
            gCATATCTCAAGATATCCATTATGATGTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCt                                                                                                                         <  1:69471/141‑1 (MQ=255)
             cATATCTCAAGATATCCATTATGATGTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCtt                                                                                                                        >  1:537491/1‑141 (MQ=255)
                atCTCAAGATATCCATTATGATGTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTc                                                                                                                     <  1:343124/141‑1 (MQ=255)
                 gctcAAGATATCCATTATGATGTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCt                                                                                                                    >  2:428475/2‑141 (MQ=255)
                      aGATATCCATTATGATGTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATTGATTTTCAACATTAGGAAATGATCAATCTTGTCTATAcc                                                                                                               <  1:1062977/141‑1 (MQ=255)
                      aGATATCCATTATGATGTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATAcc                                                                                                               >  2:440628/1‑141 (MQ=255)
                      aGATATCCATTATGATGTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATAcc                                                                                                               <  1:203260/141‑1 (MQ=255)
                      aGATATCCATTATGATGTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATAcc                                                                                                               <  2:842906/141‑1 (MQ=255)
                         tatCCATTATGATGTAGATGGCCAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCCTGTCTATACCTGc                                                                                                            <  2:848614/141‑1 (MQ=255)
                         tatCCATTATGATGTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGc                                                                                                            >  1:1056377/1‑141 (MQ=255)
                         tatCCATTATGATGTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGc                                                                                                            >  2:307415/1‑141 (MQ=255)
                          atCCATTATGATGTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCt                                                                                                           <  1:41472/141‑1 (MQ=255)
                            ccATTATGATGTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTa                                                                                                         <  1:802144/141‑1 (MQ=255)
                            ccATTATGATGTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATAGTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTa                                                                                                         >  1:556016/1‑141 (MQ=255)
                             cATTATGATGTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTAt                                                                                                        <  2:282791/141‑1 (MQ=255)
                               ttATGATGTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTc                                                                                                      <  2:722250/141‑1 (MQ=255)
                                tATGATGTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATATTGTCTATACCTGCTTATTCt                                                                                                     >  2:541031/1‑141 (MQ=255)
                                      gTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAAcc                                                                                               >  2:311337/1‑141 (MQ=255)
                                      gTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAAATACAAATATGTGTGCATGAATTAGTTTTACATATTTTCAACATTAAGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAAcc                                                                                               <  1:1067006/141‑1 (MQ=255)
                                        aGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTATCTATACCTGCTTATTCTTTAACCtt                                                                                             >  2:858556/1‑141 (MQ=255)
                                         gATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTg                                                                                            >  1:166476/1‑141 (MQ=255)
                                            ggACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGttt                                                                                         >  1:809277/1‑141 (MQ=255)
                                                   ttttATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACt                                                                                  >  1:609550/1‑141 (MQ=255)
                                                    tttATATTAAGAAGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTa                                                                                 <  2:1026399/141‑1 (MQ=255)
                                                      tataTTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGc                                                                               >  2:770886/1‑141 (MQ=255)
                                                       atatTAAGACGTTGTAATGCCACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGca                                                                              <  2:1050074/141‑1 (MQ=255)
                                                         atTAAGACGTTGTAATGACACGACACCTAACATACATACAAATATGTGTGCATGAATTAGTTTTAACTATTTTCAACCCTAAGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGcatt                                                                            <  1:103362/141‑1 (MQ=255)
                                                         atTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGcatt                                                                            <  1:465493/141‑1 (MQ=255)
                                                         atTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGcatt                                                                            <  1:359501/141‑1 (MQ=255)
                                                            aaGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGcattcat                                                                         >  2:355220/1‑141 (MQ=255)
                                                            aaGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGcattcat                                                                         >  1:926597/1‑141 (MQ=255)
                                                             aGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGcattcatt                                                                        >  2:743735/1‑141 (MQ=255)
                                                              gACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCAttt                                                                       <  2:798070/141‑1 (MQ=255)
                                                               aCGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACATTGTTTGATAACTAGCCTTCAtttt                                                                      >  2:316796/1‑141 (MQ=255)
                                                                 gTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTaa                                                                    >  1:964495/1‑141 (MQ=255)
                                                                 gTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTaa                                                                    >  1:155505/1‑141 (MQ=255)
                                                                   tGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATc                                                                  <  1:442828/141‑1 (MQ=255)
                                                                   tGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATc                                                                  >  1:1494/1‑141 (MQ=255)
                                                                   tGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATc                                                                  <  2:992765/141‑1 (MQ=255)
                                                                           cacgacacCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGt                                                          >  1:1001108/1‑141 (MQ=255)
                                                                              gacacCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGtttt                                                       >  2:1041179/1‑141 (MQ=255)
                                                                                cacCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTaa                                                     <  2:461174/141‑1 (MQ=255)
                                                                                cacCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTaa                                                     >  2:723011/1‑141 (MQ=255)
                                                                                 acCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGACCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTaaa                                                    <  1:60032/141‑1 (MQ=255)
                                                                                     aaCATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTAAAGGTc                                                <  1:963467/141‑1 (MQ=255)
                                                                                     aaCATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTAAAGGTc                                                >  1:247018/1‑141 (MQ=255)
                                                                                     aaCATAAATACAAATATGTGTGCATGAATTAGTTTTACATATTTTCAACATTAAGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTAAAGGTc                                                <  1:29074/141‑1 (MQ=255)
                                                                                       cATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTAAAGGTCtt                                              >  2:834149/1‑141 (MQ=255)
                                                                                          atttATAATTATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTAAAGGTCTTCAt                                           >  2:349036/3‑141 (MQ=255)
                                                                                          aaTTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTAAAGGTCTTCAt                                           <  1:938287/141‑1 (MQ=255)
                                                                                            ttATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTAAAGGTCTTCAttt                                         <  2:656203/141‑1 (MQ=255)
                                                                                              ataAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTAAAGGTCTTCATTTGt                                       <  1:59270/141‑1 (MQ=255)
                                                                                              ataAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTAAAGGTCTTCATTTGt                                       >  2:93586/1‑141 (MQ=255)
                                                                                                 aaTATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTAAAGGTCTTCATTTGTCGa                                    >  1:948905/1‑141 (MQ=255)
                                                                                                 aaTATGTGTGCATGAATTAGTTTTACATATTTTCAACATTAAGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTAAAGGTCTTCATTTGTCGa                                    >  1:198461/1‑141 (MQ=255)
                                                                                                   tatGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTAAAGGTCTTCATTTGTCGATa                                  >  1:996181/1‑141 (MQ=255)
                                                                                                    atGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTAAAGGTCTTCATTTGTCGATAg                                 <  2:764778/141‑1 (MQ=255)
                                                                                                     tgtgtgCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTAAAGGTCTTCATTTGTCGATAGc                                <  2:996181/141‑1 (MQ=255)
                                                                                                                aTTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTAAAGGTCTTCATTTGTCGATAGCCCATAAAATGa                     <  1:963440/141‑1 (MQ=255)
                                                                                                                   aGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGGTTTAAAGGTCTTCATTTGTCGATAGCCCATAAAATGatat                  <  2:935666/141‑1 (MQ=255)
                                                                                                                              ttttCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTAAAGGTCTTCATTTGTCGATAGCCCATAAAATGATATTTACGTTCaat       >  1:57667/1‑141 (MQ=255)
                                                                                                                              ttttCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTAAAGGTCTTCATTTGTCGATAGCCCATAAAATGATATTTACGTTCaat       >  2:316631/1‑141 (MQ=255)
                                                                                                                                ttCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTAAAGGTCTTCATTTGTCGATAGCCCATAAAATGATATTTACGTTCaatca     >  1:188119/1‑141 (MQ=255)
                                                                                                                                ttCAACATTAAGAAATGATCAATCTTGTCTATACCTGCTTATTCTTCAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTAAAGGTCTTCATTTGTCGATAGCCCATAAAATGATATTTACGTTCaatca     >  2:723285/1‑141 (MQ=255)
                                                                                                                                  cAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTAAAGGTCTTCATTTGTCGATAGCCCATAAAATGATATTTACGTTCaatcata   <  1:131613/141‑1 (MQ=255)
                                                                                                                                   aaCATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTAAAGGTCTTCATTTGTCGATAGCCCATAAAATGATATTTACGTTCaatcatat  >  2:1058221/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                     
AGATATCTTAAAGCATATCTCAAGATATCCATTATGATGTAGATGGACAGTTTTTATATTAAGACGTTGTAATGACACGACACCTAACATAATTATAAATATGTGTGCATGAATTAGTTTTATATATTTTCAACATTAGGAAATGATCAATCTTGTCTATACCTGCTTATTCTTTAACCTTGTTTGATAACTAGCATTCATTTTAATCAATTCAGTTTTAAAGGTCTTCATTTGTCGATAGCCCATAAAATGATATTTACGTTCAATCATAT  >  CP000730/1521470‑1521741

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