Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,421,606 T→C 9.3% 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 9.3% 662.1 / 40.0 223intergenic (+42/+98)pepF1_2/phoUOligoendopeptidase F, plasmid/Phosphate‑specific transport system accessory protein PhoU
Reads supporting (aligned to +/- strand):  ref base T (106/95);  new base C (15/6);  total (122/101)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 245 total aligned reads displayed.

CGATGAAGTTGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAA  >  USA300TCH1516_ALE/1421516‑1421696
                                                                                          |                                                                                           
cGATGAAGTTGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGAttt                                                                                             <  2:567602/91‑1 (MQ=255)
 gATGAAGTTGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGAtttt                                                                                            <  1:3118747/91‑1 (MQ=255)
     aaGTTGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTcaca                                                                                       <  1:1040380/91‑1 (MQ=255)
      aGTTGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACAt                                                                                      >  2:1597410/1‑91 (MQ=255)
       gTTGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACAtt                                                                                     >  1:1620390/1‑91 (MQ=255)
         tGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTAc                                                                                   <  1:3266967/91‑1 (MQ=255)
          gAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACt                                                                                  >  2:1573633/1‑91 (MQ=255)
          gAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACt                                                                                  <  1:3607987/91‑1 (MQ=255)
           aaaaaTTAAGAGATGAGATTGAGCAAGCAAAGAGCGAAGAGGAGAGACAATTTATAGTAGAACAAGATAAACCAAGATTTTCACATTACTa                                                                                 <  2:3273185/91‑1 (MQ=255)
            aaaaTTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTaa                                                                                >  2:3594422/1‑91 (MQ=255)
              aaTTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAAtt                                                                              >  2:2469107/1‑91 (MQ=255)
              aaTTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAAtt                                                                              <  1:459555/91‑1 (MQ=255)
               aTTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtg                                                                             >  1:537807/1‑91 (MQ=255)
                ttAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtgt                                                                            <  2:206763/91‑1 (MQ=255)
                ttAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtgt                                                                            <  2:1954013/91‑1 (MQ=255)
                ttAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtgt                                                                            <  1:1505755/91‑1 (MQ=255)
                 tAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtgtg                                                                           <  1:1485731/91‑1 (MQ=255)
                 tAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtgtg                                                                           <  1:3901370/91‑1 (MQ=255)
                 tAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTCTCACAATTAGTAAtgtg                                                                            <  2:311077/91‑4 (MQ=37)
                 tAACAGATGAAATGGAGAAAGCAAATAAATAAGATGCGAGACAATTGATAGTAGAACAATATAAACCAAGATTTTCACATTACTAAGtgtg                                                                           <  2:239089/91‑1 (MQ=255)
                  aaCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTga                                                                          >  1:1740393/1‑91 (MQ=255)
                   tcAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTgag                                                                         <  2:3927884/90‑1 (MQ=255)
                   aCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTgag                                                                         >  2:1849303/1‑91 (MQ=255)
                   aCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTgag                                                                         >  1:3499566/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)
                      gATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAAt                                                                      <  2:657978/91‑1 (MQ=255)
                       aTGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATc                                                                     <  2:2545167/91‑1 (MQ=255)
                       aTGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATc                                                                     <  2:83355/91‑1 (MQ=255)
                        tGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCt                                                                    <  1:2117145/91‑1 (MQ=255)
                        tGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTCTCACAATTAGTAATGTGAAAATCt                                                                    <  2:1745933/91‑1 (MQ=21)
                        tGAAATTAAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATACTCACAATTAGTAATGTGAAAATCt                                                                    <  2:643981/91‑1 (MQ=21)
                           aaTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTgg                                                                 <  2:3170265/91‑1 (MQ=255)
                              tGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGttt                                                              >  2:3463679/1‑91 (MQ=255)
                               gAGCAAGCAAATAACTAAGATGCGAGACAAGTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGtttt                                                             <  1:2802995/91‑1 (MQ=255)
                                 gcaagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGtttttt                                                           >  1:3024938/1‑91 (MQ=255)
                                 gcaagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGtttttt                                                           <  1:3432821/91‑1 (MQ=255)
                                  caagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGttttttt                                                          >  1:2152884/1‑91 (MQ=255)
                                  caagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGttttttt                                                          >  2:3564908/1‑91 (MQ=255)
                                        aaTAACTACGATGCGAGACAATTTATAGTCGAACAATATAAAACAAGATATTCACAACAATAATTGTACGAATCCTGGGTCTTCTGAgggg                                                    >  2:422870/1‑90 (MQ=255)
                                            aCTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTAt                                                >  1:1185240/1‑91 (MQ=255)
                                            aCTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTGTGAGGGTTTAt                                                >  1:3585281/1‑91 (MQ=255)
                                             cTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGTGTTTATg                                               >  2:1362820/1‑91 (MQ=255)
                                              tAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGa                                              >  1:670385/1‑91 (MQ=255)
                                                aGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaat                                            >  2:2492960/1‑91 (MQ=255)
                                                aGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaat                                            <  1:1580552/91‑1 (MQ=255)
                                                aGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaat                                            <  2:1276230/91‑1 (MQ=255)
                                                aGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaat                                            <  2:2890769/91‑1 (MQ=255)
                                                  aTGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaataa                                          <  1:673997/91‑1 (MQ=255)
                                                   tGCGAGACAATTTATAGTAGAACAATATAAACTAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaataat                                         <  1:3370748/91‑1 (MQ=255)
                                                       agaCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGaa                                     >  1:2576916/1‑91 (MQ=255)
                                                       agaCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGaa                                     <  2:136560/91‑1 (MQ=255)
                                                        gaCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGaaa                                    <  1:3204035/91‑1 (MQ=255)
                                                           aaTTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatat                                 <  1:2252011/91‑1 (MQ=255)
                                                           aaTTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatat                                 >  1:229456/1‑91 (MQ=255)
                                                            aTTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatata                                >  1:1328877/1‑91 (MQ=255)
                                                              ttATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatatata                              <  1:1535582/91‑1 (MQ=255)
                                                               tataGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTACGAATAATAGAAatatatat                             >  2:109623/1‑91 (MQ=255)
                                                                ataGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATAtt                            >  2:2673720/1‑91 (MQ=255)
                                                                ataGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATAtt                            >  1:1408209/1‑91 (MQ=255)
                                                                 tagtagAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATAttt                           >  1:3928823/1‑91 (MQ=255)
                                                                  agtagAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGATTAATAGAAATATATATTTa                          >  2:3072020/1‑91 (MQ=255)
                                                                  agtagAACAATATAAACCAAGATTTTCACATTACTAACTGTGAGAATCTTGCTTTTTTTGAGGGTGTATGAATAATAGAAATATACATTTa                          >  1:3680982/1‑91 (MQ=255)
                                                                     agAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCg                       >  2:1033660/1‑91 (MQ=255)
                                                                     agAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCg                       >  1:3227540/1‑91 (MQ=255)
                                                                     agAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGCTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCg                       >  2:2327535/1‑91 (MQ=255)
                                                                       aaCAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATTTTTATCGTa                     >  1:3802335/1‑91 (MQ=255)
                                                                        aCAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTaa                    <  2:367644/91‑1 (MQ=255)
                                                                        aCAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTaa                    <  1:2957199/91‑1 (MQ=255)
                                                                          aaTATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAgg                  <  1:1709035/91‑1 (MQ=255)
                                                                           atatAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAggg                 <  2:379880/91‑1 (MQ=255)
                                                                             ataAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAt               >  2:3319844/1‑91 (MQ=255)
                                                                             ataAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAt               >  2:901977/1‑91 (MQ=255)
                                                                             ataAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAt               <  1:1068753/91‑1 (MQ=255)
                                                                             ataAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAt               <  1:3729485/91‑1 (MQ=255)
                                                                               aaaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAttt             <  2:3738528/91‑1 (MQ=255)
                                                                               aaaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAttt             <  2:2122839/91‑1 (MQ=255)
                                                                                aaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAtttt            >  2:1781863/1‑91 (MQ=255)
                                                                                aaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTATTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAgttt            >  2:2142904/1‑91 (MQ=255)
                                                                                aaCCAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTGAGGGTTTATGAATAATAGACATATATATTTATCGTAAGGGATTTta           >  1:2871837/1‑91 (MQ=18)
                                                                                aaCCAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTta           >  1:1002238/1‑91 (MQ=18)
                                                                                 aCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTta           <  2:1661632/91‑1 (MQ=255)
                                                                                 aCCAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAtttga           >  1:819584/1‑89 (MQ=21)
                                                                                  ccAAGATTTTCACATTACTAATTGTTAGAATCTTGGGTTTTTTGAGGGTTTATGAATAAGAGAAATATATATTTATCGTGAGGGATTTtat          <  2:220172/91‑1 (MQ=255)
                                                                                  ccAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTtat          >  1:1912027/1‑91 (MQ=21)
                                                                                  ccAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTtat          >  1:2513427/1‑91 (MQ=21)
                                                                                  ccAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTtat          >  1:2913691/1‑91 (MQ=21)
                                                                                     aGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGCGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGt       >  2:3003398/1‑91 (MQ=21)
                                                                                     aGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTCAGGGATTTTATAGt       >  2:2627261/1‑91 (MQ=21)
                                                                                       aTTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaa     >  1:634837/1‑91 (MQ=255)
                                                                                       aTTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaa     <  1:2390302/91‑1 (MQ=255)
                                                                                       aTTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaa     <  1:1847301/91‑1 (MQ=255)
                                                                                        ttttCACTTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaat    <  2:3510081/91‑1 (MQ=255)
                                                                                        ttttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaat    >  1:731550/1‑91 (MQ=255)
                                                                                        ttttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaat    <  2:1876211/91‑1 (MQ=255)
                                                                                        ttttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaat    >  1:3400666/1‑91 (MQ=255)
                                                                                         tttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTCTGACTAACAGAAACATACACTTATCGTAAGGGATTTTATAGtaata   >  2:3402216/1‑91 (MQ=255)
                                                                                          ttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaataa  >  2:3395781/1‑91 (MQ=255)
                                                                                          ttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaataa  >  1:1281474/1‑91 (MQ=255)
                                                                                          |                                                                                           
CGATGAAGTTGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAA  >  USA300TCH1516_ALE/1421516‑1421696

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: