Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,421,606 T→C 6.1% intergenic (+42/+98) pepF1_2 → / ← phoU Oligoendopeptidase F, plasmid/Phosphate‑specific transport system accessory protein PhoU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,421,6060TC 6.1% 738.3 / 23.3 229intergenic (+42/+98)pepF1_2/phoUOligoendopeptidase F, plasmid/Phosphate‑specific transport system accessory protein PhoU
Reads supporting (aligned to +/- strand):  ref base T (115/100);  new base C (7/7);  total (122/107)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 235 total aligned reads displayed.

CGATGAAGTTGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAA  >  USA300TCH1516_ALE/1421516‑1421696
                                                                                          |                                                                                           
 gATGAAGTTGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGAtttt                                                                                            >  1:983951/1‑91 (MQ=255)
  aTGAAGTTGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTc                                                                                           >  1:2854131/1‑91 (MQ=255)
       gTTGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACAtt                                                                                     >  1:237135/1‑91 (MQ=255)
       gTTGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACAtt                                                                                     <  2:3768687/91‑1 (MQ=255)
         tGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTAc                                                                                   <  1:3144761/91‑1 (MQ=255)
          gAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACt                                                                                  >  2:3282792/1‑91 (MQ=255)
          gAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACt                                                                                  >  1:3733308/1‑91 (MQ=255)
             aaaTTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAAt                                                                               >  1:759935/1‑91 (MQ=255)
                ttAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtgt                                                                            <  2:1060382/91‑1 (MQ=255)
                 tAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtgtg                                                                           <  2:3330748/91‑1 (MQ=255)
                  aaCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTga                                                                          <  2:1303777/91‑1 (MQ=255)
                  aaCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTga                                                                          >  1:3016727/1‑91 (MQ=255)
                  aaCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTga                                                                          >  1:3016502/1‑91 (MQ=255)
                    cAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTgaga                                                                        >  1:1056877/1‑91 (MQ=255)
                      gATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTAc                                                                                   <  2:628239/78‑1 (MQ=255)
                      gATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTAc                                                                                   >  1:628239/1‑78 (MQ=255)
                      gATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAAt                                                                      >  1:1261902/1‑91 (MQ=255)
                       aTGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATc                                                                     <  1:3799758/91‑1 (MQ=255)
                        tGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCt                                                                    <  1:2689913/91‑1 (MQ=255)
                         gAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCtt                                                                   >  2:2502837/1‑91 (MQ=255)
                         gAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTCTCACAATTAGTAATGTGAAAATCtt                                                                   <  2:3791712/91‑1 (MQ=21)
                           aaTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTgg                                                                 >  2:1162201/1‑91 (MQ=255)
                           aaTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTgg                                                                 <  2:2021413/91‑1 (MQ=255)
                            aTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGt                                                                >  1:1585678/1‑91 (MQ=255)
                            aTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGt                                                                <  2:410882/91‑1 (MQ=255)
                            aTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGt                                                                <  1:1392549/91‑1 (MQ=255)
                               gAGCAAGCAAATTACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGtttt                                                             >  1:2974723/1‑91 (MQ=255)
                               gAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGtttt                                                             >  1:963043/1‑91 (MQ=255)
                                agcaagcaAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGttttt                                                            >  1:2070522/1‑91 (MQ=255)
                                 gcaagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGtttttt                                                           >  2:98825/1‑91 (MQ=255)
                                 gcaagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGtttttt                                                           <  1:2300568/91‑1 (MQ=255)
                                 gcaagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGtttttt                                                           <  2:742843/91‑1 (MQ=255)
                                   aagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTg                                                         <  2:1716847/91‑1 (MQ=255)
                                   aagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTg                                                         <  1:425936/91‑1 (MQ=255)
                                   aagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTg                                                         <  2:1627873/91‑1 (MQ=255)
                                   aagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTg                                                         <  1:2140016/91‑1 (MQ=255)
                                    agcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGa                                                        >  2:1505504/1‑91 (MQ=255)
                                    agcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGa                                                        >  1:3413047/1‑91 (MQ=255)
                                    agcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGa                                                        >  1:3538580/1‑91 (MQ=255)
                                    agcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGa                                                        >  1:656094/1‑91 (MQ=255)
                                        aaTAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCtt                                                                   >  2:1475352/1‑76 (MQ=255)
                                          tAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGttt                                                  <  2:1213565/91‑1 (MQ=255)
                                          tAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGttt                                                  >  1:481916/1‑91 (MQ=255)
                                            aCTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTAt                                                <  1:2021147/91‑1 (MQ=255)
                                            aCTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTAt                                                >  1:1319732/1‑91 (MQ=255)
                                            aCTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTAt                                                <  2:2854131/91‑1 (MQ=255)
                                             cTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATg                                               <  1:286554/91‑1 (MQ=255)
                                              tAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGa                                              >  1:2253618/1‑91 (MQ=255)
                                               aaGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaa                                             <  1:2650486/91‑1 (MQ=255)
                                               aaGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaa                                             <  2:7832/91‑1 (MQ=255)
                                               aaGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaa                                             <  2:5679/91‑1 (MQ=255)
                                               aaGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaa                                             <  1:3269027/91‑1 (MQ=255)
                                                aGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGTTTTATGaat                                            <  1:1955864/91‑1 (MQ=255)
                                                 gATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaata                                           >  1:1567740/1‑91 (MQ=255)
                                                  aTGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaataa                                          >  1:3704203/1‑91 (MQ=255)
                                                  aTGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaataa                                          >  2:1052695/1‑91 (MQ=255)
                                                    gCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAataata                                        <  1:1450448/91‑1 (MQ=255)
                                                    gCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAataata                                        >  1:2599556/1‑91 (MQ=255)
                                                        gaCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGaaa                                    >  1:2558381/1‑91 (MQ=255)
                                                         aCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAat                                   <  1:3447006/91‑1 (MQ=255)
                                                          cAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAata                                  >  1:2417407/1‑91 (MQ=255)
                                                           aaTTTATAGTAGAACAGTATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatat                                 <  2:3704203/91‑1 (MQ=255)
                                                           aaTTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatat                                 >  2:3640783/1‑91 (MQ=255)
                                                             tttATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatatat                               <  1:222365/91‑1 (MQ=255)
                                                             tttATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatatat                               <  1:2890732/91‑1 (MQ=255)
                                                             tttATAGTAGAACAATATAAACCAAGATTTTCACATCACTAATTGTGAGAATCATGGTTTTTTTGAGGGTTTATGAATAATAGAAatatat                               >  2:1250310/1‑91 (MQ=255)
                                                              ttATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatatata                              >  1:1821152/1‑91 (MQ=255)
                                                              ttATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatatata                              >  1:1999226/1‑91 (MQ=255)
                                                               tataGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatatatat                             >  2:809640/1‑91 (MQ=255)
                                                               tataGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatatatat                             <  2:299389/91‑1 (MQ=255)
                                                                ataGTAGAACAATATAAACCAAGATTTTCACATAACTGATTGTGAGCATCTTGGTTTTTCTGAGGGTTTATGAATAATAGAAATATATAtt                            >  1:3095873/1‑91 (MQ=255)
                                                                 tagtagAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATAttt                           >  2:2255983/1‑91 (MQ=255)
                                                                 tagtagAACAATATAAACCAAGATTTTCACATTACTAATTGAGAGAATCTTGGCTATTTTGAGGGTTTATGAATAATAGAAATATATAttt                           >  2:3480803/1‑91 (MQ=255)
                                                                    tagAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATc                        >  2:3181104/1‑91 (MQ=255)
                                                                    tagAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATc                        >  1:1970473/1‑91 (MQ=255)
                                                                    tagAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATc                        <  2:2172305/91‑1 (MQ=255)
                                                                    tagAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATc                        >  1:3865463/1‑91 (MQ=255)
                                                                      gAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGt                      <  2:1381486/91‑1 (MQ=255)
                                                                        aCAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTaa                    <  2:2621117/91‑1 (MQ=255)
                                                                        aCAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTaa                    <  1:3181104/91‑1 (MQ=255)
                                                                         cAATATAAACCAAGATTTTCACATTGCTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAg                   <  1:3621161/91‑1 (MQ=255)
                                                                         cAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAg                   <  1:1604439/91‑1 (MQ=255)
                                                                          aaTATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAgg                  <  1:2244618/91‑1 (MQ=255)
                                                                          aaTATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAgg                  >  2:2244556/1‑91 (MQ=255)
                                                                           atatAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAggg                 >  1:1911140/1‑91 (MQ=255)
                                                                             ataAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAt               >  1:2762473/1‑91 (MQ=255)
                                                                             ataAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTATATGAATAATAGAAATATATATTTATCGTAAGGGAt               >  2:2691066/1‑91 (MQ=255)
                                                                             aaaaaCCAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAt               >  2:1345938/3‑91 (MQ=18)
                                                                                aaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAtttt            <  1:2560467/91‑1 (MQ=255)
                                                                                aaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAtttt            >  1:2511447/1‑91 (MQ=255)
                                                                                 aCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTta           >  2:2087492/1‑91 (MQ=255)
                                                                                  ccAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTtat          >  1:844317/1‑91 (MQ=21)
                                                                                    aaGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAg        >  2:2392057/1‑91 (MQ=21)
                                                                                     aGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGt       >  1:824081/1‑91 (MQ=255)
                                                                                      gATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGta      >  2:949956/1‑91 (MQ=255)
                                                                                      gATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAATTTATCGta      >  1:3490709/1‑91 (MQ=255)
                                                                                       aTTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaa     >  2:11936/1‑91 (MQ=255)
                                                                                       aTTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaa     >  2:2100352/1‑91 (MQ=255)
                                                                                        ttttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaat    <  2:1911140/91‑1 (MQ=255)
                                                                                         tttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaata   >  2:3687224/1‑91 (MQ=255)
                                                                                          |                                                                                           
CGATGAAGTTGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAA  >  USA300TCH1516_ALE/1421516‑1421696

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