Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,421,610:1 +A 6.1% intergenic (+46/+94) 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,6101.A 6.1% 983.1 / 51.5 230intergenic (+46/+94)pepF1_2/phoUOligoendopeptidase F, plasmid/Phosphate‑specific transport system accessory protein PhoU
Reads supporting (aligned to +/- strand):  ref base . (117/99);  new base A (7/7);  total (124/106)
Fisher's exact test for biased strand distribution p-value = 7.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 239 total aligned reads displayed.

AAGTTGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTT  >  USA300TCH1516_ALE/1421521‑1421700
                                                                                          |                                                                                          
aaGTTGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTcaca                                                                                           <  1:3758586/91‑1 (MQ=255)
  gTTGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACAtt                                                                                         >  1:237078/1‑91 (MQ=255)
  gTTGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACAtt                                                                                         <  2:3756717/91‑1 (MQ=255)
    tGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTAc                                                                                       <  1:2126890/91‑1 (MQ=255)
    tGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTAc                                                                                       <  1:3144761/91‑1 (MQ=255)
     gAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACt                                                                                      >  2:388111/1‑91 (MQ=255)
     gAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACt                                                                                      >  2:3676975/1‑91 (MQ=255)
     gAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACt                                                                                      >  2:3281496/1‑91 (MQ=255)
        aaaTTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAAt                                                                                   >  1:759935/1‑91 (MQ=255)
           ttAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtgt                                                                                <  1:1927127/91‑1 (MQ=255)
            tAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtgtg                                                                               <  2:3330748/91‑1 (MQ=255)
             aaCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTga                                                                              >  1:3016727/1‑91 (MQ=255)
             aaCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTga                                                                              >  1:3016502/1‑91 (MQ=255)
             aaCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTga                                                                              >  1:254383/1‑91 (MQ=255)
               cAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTgaga                                                                            >  1:1056877/1‑91 (MQ=255)
                aGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGaa                                                                           <  2:3040818/91‑1 (MQ=255)
                 gATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTAc                                                                                       >  1:628239/1‑78 (MQ=255)
                 gATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTAc                                                                                       <  2:628239/78‑1 (MQ=255)
                 gATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAAt                                                                          >  1:1261902/1‑91 (MQ=255)
                  aTGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATc                                                                         <  1:3799758/91‑1 (MQ=255)
                   tGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCt                                                                        >  2:1426929/1‑91 (MQ=255)
                   tGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCt                                                                        <  1:2689913/91‑1 (MQ=255)
                   tGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCt                                                                        >  2:2486104/1‑91 (MQ=255)
                    gAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCtt                                                                       <  1:1852807/91‑1 (MQ=255)
                    gAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTCTCACAATTAGTAATGTGAAAATCtt                                                                       <  2:3791712/91‑1 (MQ=21)
                     aaaTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTg                                                                      >  2:1181162/1‑91 (MQ=255)
                      aaTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTgg                                                                     <  2:2021413/91‑1 (MQ=255)
                      aaTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTgg                                                                     >  2:1162201/1‑91 (MQ=255)
                      aaTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTCTCACAATTAGTAATGTGAAAATCTTgg                                                                     <  1:662624/91‑1 (MQ=21)
                       aTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGt                                                                    >  1:1585678/1‑91 (MQ=255)
                       aTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGt                                                                    <  1:1392549/91‑1 (MQ=255)
                        ttGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGtt                                                                   <  2:1059870/91‑1 (MQ=255)
                          gAGCAAGCAAATTACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGtttt                                                                 >  1:2974723/1‑91 (MQ=255)
                          gAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGtttt                                                                 >  1:961022/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:1627873/91‑1 (MQ=255)
                              aagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTg                                                             <  2:1716847/91‑1 (MQ=255)
                               agcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGa                                                            >  1:3413047/1‑91 (MQ=255)
                               agcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGa                                                            >  2:1505504/1‑91 (MQ=255)
                               agcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGa                                                            >  1:3538580/1‑91 (MQ=255)
                                     tAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGttt                                                      <  2:1213565/91‑1 (MQ=255)
                                     tAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGttt                                                      >  1:481916/1‑91 (MQ=255)
                                     tAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGttt                                                      >  1:3233401/1‑91 (MQ=255)
                                      aaCTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTa                                                     >  2:274988/1‑91 (MQ=255)
                                       aCTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTAt                                                    >  2:265745/1‑91 (MQ=255)
                                       aCTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTAt                                                    <  1:2021147/91‑1 (MQ=255)
                                       aCTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTAt                                                    >  1:1319732/1‑91 (MQ=255)
                                        cTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATg                                                   <  1:286554/91‑1 (MQ=255)
                                         tAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGa                                                  >  1:2253618/1‑91 (MQ=255)
                                          aaGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaa                                                 <  1:2650486/91‑1 (MQ=255)
                                          aaGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaa                                                 <  1:3269027/91‑1 (MQ=255)
                                          aaGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaa                                                 <  2:5679/91‑1 (MQ=255)
                                          aaGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaa                                                 <  2:7832/91‑1 (MQ=255)
                                            gATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaata                                               >  1:1567740/1‑91 (MQ=255)
                                             aTGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaataa                                              >  2:1052695/1‑91 (MQ=255)
                                             aTGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaataa                                              >  1:3704203/1‑91 (MQ=255)
                                               gCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAataata                                            <  1:1450448/91‑1 (MQ=255)
                                               gCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAataata                                            >  1:2599556/1‑91 (MQ=255)
                                                    aCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAat                                       <  1:3447006/91‑1 (MQ=255)
                                                     cAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAata                                      >  1:2417407/1‑91 (MQ=255)
                                                     cAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAata                                      >  2:3702482/1‑91 (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:1999158/1‑91 (MQ=255)
                                                         ttATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatatata                                  >  1:1821152/1‑91 (MQ=255)
                                                          tataGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatatatat                                 <  2:299389/91‑1 (MQ=255)
                                                          tataGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatatatat                                 >  2:809640/1‑91 (MQ=255)
                                                           ataGTAGAACAATATAAACCAAGATTTTCACATAACTGATTGTGAGCATCTTGGTTTTTCTGAGGGTTTATGAATAATAGAAATATATAtt                                >  1:3095873/1‑91 (MQ=255)
                                                               tagAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATc                            >  1:3865463/1‑91 (MQ=255)
                                                               tagAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATc                            >  1:1970473/1‑91 (MQ=255)
                                                                 gAACAATGTAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGt                          <  2:3456622/91‑1 (MQ=255)
                                                                 gAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGt                          <  2:1381486/91‑1 (MQ=255)
                                                                   aCAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTaa                        <  2:2621117/91‑1 (MQ=255)
                                                                    cAATATAAACCAAGATTTTCACATTGCTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAg                       <  1:3621161/91‑1 (MQ=255)
                                                                     aaTATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAgg                      <  1:2244618/91‑1 (MQ=255)
                                                                     aaTATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAgg                      >  2:2244618/1‑91 (MQ=255)
                                                                        ataAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAt                   >  1:2762473/1‑91 (MQ=255)
                                                                        aaaaaCCAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAt                   >  2:1345938/3‑91 (MQ=18)
                                                                          aaaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAttt                 <  2:2177965/91‑1 (MQ=255)
                                                                           aaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAtttt                <  1:2560467/91‑1 (MQ=255)
                                                                           aaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAtttt                >  1:2511447/1‑91 (MQ=255)
                                                                            aCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTta               >  2:2087492/1‑91 (MQ=255)
                                                                               aaGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAg            <  2:1266983/91‑1 (MQ=255)
                                                                               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:2100352/1‑91 (MQ=255)
                                                                                   ttttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaat        <  2:1911140/91‑1 (MQ=255)
                                                                                      tCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAAc     >  1:9006/1‑91 (MQ=255)
                                                                                        gcaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGATTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACtt   <  2:3842895/90‑1 (MQ=255)
                                                                                        acaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACtt   >  2:3558765/1‑91 (MQ=255)
                                                                                        acaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACtt   >  1:1592152/1‑91 (MQ=255)
                                                                                         caTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACttt  <  1:981237/91‑1 (MQ=255)
                                                                                         caTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACtta  <  2:2201414/91‑2 (MQ=255)
                                                                                          |                                                                                          
AAGTTGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTT  >  USA300TCH1516_ALE/1421521‑1421700

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: