Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,421,610:1 +A 8.3% 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 8.3% 963.3 / 78.5 240intergenic (+46/+94)pepF1_2/phoUOligoendopeptidase F, plasmid/Phosphate‑specific transport system accessory protein PhoU
Reads supporting (aligned to +/- strand):  ref base . (117/103);  new base A (15/5);  total (132/108)
Fisher's exact test for biased strand distribution p-value = 6.51e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 258 total aligned reads displayed.

AAGTTGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTT  >  USA300TCH1516_ALE/1421521‑1421700
                                                                                          |                                                                                          
aaGTTGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTcaca                                                                                           <  1:1040380/91‑1 (MQ=255)
 aGTTGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACAt                                                                                          >  2:1597515/1‑91 (MQ=255)
  gTTGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACAtt                                                                                         >  1:1620390/1‑91 (MQ=255)
  gTTGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACAtt                                                                                         >  1:3815762/1‑91 (MQ=255)
    tGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTAc                                                                                       <  1:3266967/91‑1 (MQ=255)
    tGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTAc                                                                                       >  2:3863574/1‑91 (MQ=255)
     gAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACt                                                                                      >  2:1573633/1‑91 (MQ=255)
     gAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACt                                                                                      >  1:1461984/1‑91 (MQ=255)
     gAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACt                                                                                      <  1:3607936/91‑1 (MQ=255)
      aaaaaTTAAGAGATGAGATTGAGCAAGCAAAGAGCGAAGAGGAGAGACAATTTATAGTAGAACAAGATAAACCAAGATTTTCACATTACTa                                                                                     <  2:3273185/91‑1 (MQ=255)
       aaaaTTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTaa                                                                                    >  2:3594511/1‑91 (MQ=255)
         aaTTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAAtt                                                                                  <  2:2372484/91‑1 (MQ=255)
         aaTTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAAtt                                                                                  <  1:459555/91‑1 (MQ=255)
           ttAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtgt                                                                                <  2:2874761/91‑1 (MQ=255)
           ttAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtgt                                                                                >  2:2496608/1‑91 (MQ=255)
            tAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtgtg                                                                               >  2:3454342/1‑91 (MQ=255)
            tAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtgtg                                                                               <  1:1485731/91‑1 (MQ=255)
            tAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTCTCACAATTAGTAAtgtg                                                                                <  2:311077/91‑4 (MQ=37)
              aCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTgag                                                                             >  2:2474950/1‑91 (MQ=255)
              aCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTgag                                                                             >  1:3640930/1‑91 (MQ=255)
              aCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTgag                                                                             >  2:2016187/1‑91 (MQ=255)
              aCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTgag                                                                             >  2:1339456/1‑91 (MQ=255)
              aCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTgag                                                                             >  2:650010/1‑91 (MQ=255)
                aGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCCCATTACTAATTGTGAGaa                                                                           >  2:2783442/1‑91 (MQ=255)
                aGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGaa                                                                           >  2:190714/1‑91 (MQ=255)
                  aTGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATc                                                                         <  2:2545167/91‑1 (MQ=255)
                   tGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCt                                                                        <  1:2117145/91‑1 (MQ=255)
                   tGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTCTCACAATTAGTAATGTGAAAATCt                                                                        <  2:644094/91‑1 (MQ=21)
                      aaTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTgg                                                                     <  2:1757071/91‑1 (MQ=255)
                       aTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGt                                                                    <  2:3381195/91‑1 (MQ=255)
                        ttGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTCTCACAATTAGTAATGTGAAAATCTTGGtt                                                                   <  1:3223794/91‑1 (MQ=21)
                         tGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGCATCTTGGttt                                                                  >  2:2178888/1‑91 (MQ=255)
                         tGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGttt                                                                  >  2:3463679/1‑91 (MQ=255)
                          gAGCAAGCAAATAACTAAGATGCGAGACAAGTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGtttt                                                                 <  1:2802995/91‑1 (MQ=255)
                             caagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGttttttt                                                              <  2:488742/91‑1 (MQ=255)
                              aagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTg                                                             <  2:1857342/91‑1 (MQ=255)
                                    aTAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGtt                                                       <  2:2086903/91‑1 (MQ=255)
                                     tAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGttt                                                      <  2:3933107/91‑1 (MQ=255)
                                     tAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGttt                                                      >  1:1751573/1‑91 (MQ=255)
                                      aaCTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTa                                                     <  2:2426724/91‑1 (MQ=255)
                                       aCTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTAt                                                    <  1:192735/91‑1 (MQ=255)
                                        cTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGTGTTTATg                                                   >  2:1362820/1‑91 (MQ=255)
                                         tAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGa                                                  >  1:811382/1‑91 (MQ=255)
                                         tAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGa                                                  >  1:670385/1‑91 (MQ=255)
                                         tAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGa                                                  >  2:699608/1‑91 (MQ=255)
                                         tAAGATGCGAGACAACTTATAGTAGAACAATCTAAACCAAGTTTTTCACATTACTAATTGTGAGAATCATGGTTTTTTTGAGGGTTTATGa                                                  >  2:3009068/1‑91 (MQ=255)
                                           aGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaat                                                <  1:1580552/91‑1 (MQ=255)
                                           aGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaat                                                <  2:1276230/91‑1 (MQ=255)
                                           aGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaat                                                <  1:3002787/91‑1 (MQ=255)
                                           aGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaat                                                <  1:2882028/91‑1 (MQ=255)
                                             aTGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaataa                                              <  1:673997/91‑1 (MQ=255)
                                             aTGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaataa                                              <  1:3774040/91‑1 (MQ=255)
                                                 gagaCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGa                                          >  2:2652518/1‑91 (MQ=255)
                                                  agaCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGaa                                         >  1:2154306/1‑91 (MQ=255)
                                                  agaCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGaa                                         <  2:136560/91‑1 (MQ=255)
                                                  agaCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGaa                                         >  1:2576916/1‑91 (MQ=255)
                                                      aaTTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatat                                     <  1:2252011/91‑1 (MQ=255)
                                                       aTTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatata                                    >  1:1328877/1‑91 (MQ=255)
                                                        tttATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatatat                                   >  1:539964/1‑91 (MQ=255)
                                                          tataGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTACGAATAATAGAAatatatat                                 >  2:109623/1‑91 (MQ=255)
                                                           ataGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATAtt                                >  2:2673720/1‑91 (MQ=255)
                                                            tagtagAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATAttt                               >  1:3928823/1‑91 (MQ=255)
                                                                agAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCg                           >  2:1033660/1‑91 (MQ=255)
                                                                 gAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGt                          <  1:3709743/91‑1 (MQ=255)
                                                                  aaCAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTa                         <  1:354487/91‑1 (MQ=255)
                                                                   aCAATATAAACCAATATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTaa                        <  2:2128450/91‑1 (MQ=255)
                                                                   aCAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTaa                        <  1:2957199/91‑1 (MQ=255)
                                                                     aaTATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAgg                      <  1:1709035/91‑1 (MQ=255)
                                                                     aaTATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAgg                      <  2:3592974/91‑1 (MQ=255)
                                                                        ataAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAt                   >  2:901977/1‑91 (MQ=255)
                                                                        ataAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAt                   <  1:1068753/91‑1 (MQ=255)
                                                                          aaaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAttt                 >  2:732346/1‑91 (MQ=255)
                                                                          aaaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAttt                 <  2:3738257/91‑1 (MQ=255)
                                                                           aaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAtttt                >  2:1781863/1‑91 (MQ=255)
                                                                           aaCCAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTta               >  1:1002238/1‑91 (MQ=18)
                                                                            aCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTta               >  2:3136506/1‑91 (MQ=255)
                                                                            aCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTta               <  2:1661632/91‑1 (MQ=255)
                                                                            aCCAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTta               >  1:3440359/1‑91 (MQ=21)
                                                                             ccAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTtat              >  1:2913691/1‑91 (MQ=21)
                                                                             ccAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTtat              >  1:2513427/1‑91 (MQ=21)
                                                                                aGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGCGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGt           >  2:3003398/1‑91 (MQ=21)
                                                                                  aTTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaa         <  1:1111600/91‑1 (MQ=255)
                                                                                  aTTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaa         <  1:2390139/91‑1 (MQ=255)
                                                                                  aTTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaa         <  1:1847301/91‑1 (MQ=255)
                                                                                  aTTCTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaa         <  1:3191232/91‑1 (MQ=255)
                                                                                   ttttCACGTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaat        <  2:3510106/91‑1 (MQ=255)
                                                                                   ttttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaat        >  1:3400667/1‑91 (MQ=255)
                                                                                   ttttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaat        <  1:1515054/91‑1 (MQ=255)
                                                                                   ttttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaat        >  1:3400636/1‑91 (MQ=255)
                                                                                    tttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaata       >  2:375319/1‑91 (MQ=255)
                                                                                     ttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaataa      >  2:3395781/1‑91 (MQ=255)
                                                                                     ttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaataa      >  1:1281474/1‑91 (MQ=255)
                                                                                      tCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAAc     <  1:650010/91‑1 (MQ=255)
                                                                                       cacaTTACTAATTGTGAGAATCTTGGTTTTTTTTATTTTTTATGAATAATAGAAATATACATTTATCGTAAGGGATTTTATAGTAATAACt    >  2:2822803/1‑91 (MQ=255)
                                                                                       cacaTTACTAATTGTGAGAATCTTGGTTTTTTTGTGGGTTTATGAATAATGGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACt    <  2:22218/91‑1 (MQ=255)
                                                                                       cacaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACt    >  1:217503/1‑91 (MQ=255)
                                                                                        acaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACtt   >  1:2344572/1‑91 (MQ=255)
                                                                                        acaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACtt   >  1:1530486/1‑91 (MQ=255)
                                                                                        acaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACtt   >  1:1251581/1‑91 (MQ=255)
                                                                                        acaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACtt   <  1:694893/91‑1 (MQ=255)
                                                                                          |                                                                                          
AAGTTGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTT  >  USA300TCH1516_ALE/1421521‑1421700

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