Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 740,328:1 +G 10.5% intergenic (+555/+33) pitA_1 → / ← sle1_2 Low‑affinity inorganic phosphate transporter 1/N‑acetylmuramoyl‑L‑alanine amidase sle1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE740,3281.G10.5% 971.6 / 116.7 247intergenic (+555/+33)pitA_1/sle1_2Low‑affinity inorganic phosphate transporter 1/N‑acetylmuramoyl‑L‑alanine amidase sle1
Reads supporting (aligned to +/- strand):  ref base . (102/119);  new base G (12/14);  total (114/133)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.36e-01

Only 100 of 260 total aligned reads displayed.

AATATGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGTTCATA  >  USA300TCH1516_ALE/740237‑740418
                                                                                            |                                                                                          
 atatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGa                                                                                             >  1:1760224/1‑91 (MQ=255)
  tatgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGaa                                                                                           <  2:418244/91‑1 (MQ=255)
  tatgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGaa                                                                                           >  2:1703154/1‑91 (MQ=255)
   atgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAAt                                                                                          <  1:1296216/91‑1 (MQ=255)
   atgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAAt                                                                                          <  1:2687324/91‑1 (MQ=255)
      tatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGc                                                                                       >  1:3136304/1‑91 (MQ=255)
      tatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGc                                                                                       >  2:2472698/1‑91 (MQ=255)
       atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCt                                                                                      >  2:2074575/1‑91 (MQ=255)
       atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCt                                                                                      <  2:502416/91‑1 (MQ=255)
       atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCt                                                                                      <  1:1310746/91‑1 (MQ=255)
       atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCt                                                                                      <  2:219560/91‑1 (MQ=255)
       atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCt                                                                                      <  1:2561224/91‑1 (MQ=255)
        tgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTa                                                                                     >  2:1902410/1‑91 (MQ=255)
        tgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTa                                                                                     <  1:2702856/91‑1 (MQ=37)
          aTAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGt                                                                                   <  1:654878/91‑1 (MQ=37)
          aTAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGt                                                                                   <  2:1197998/91‑1 (MQ=37)
          aTACGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGt                                                                                   <  1:3056890/91‑1 (MQ=255)
           tAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGta                                                                                  >  1:2738335/1‑91 (MQ=255)
            aGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGtat                                                                                 <  1:2615841/91‑1 (MQ=255)
            aGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGttt                                                                                 <  1:2004877/91‑3 (MQ=25)
            aGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGttt                                                                                 <  1:113680/91‑3 (MQ=25)
             gggAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGtata                                                                                >  1:2047377/1‑91 (MQ=255)
             gggAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGtttt                                                                                <  1:3513597/91‑4 (MQ=25)
              ggAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGtatat                                                                               >  1:3575911/1‑91 (MQ=255)
              ggAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGttttt                                                                               <  2:2461802/91‑5 (MQ=25)
                 aGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATc                                                                            <  2:1395245/91‑1 (MQ=255)
                   ggTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAg                                                                          >  2:3429115/1‑91 (MQ=255)
                   ggTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAg                                                                          >  1:1851423/1‑91 (MQ=255)
                    gTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGt                                                                         <  1:10966/91‑1 (MQ=255)
                     tGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTa                                                                        <  2:1333119/91‑1 (MQ=255)
                       agaTAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACa                                                                      >  2:2640949/1‑91 (MQ=255)
                       agaTAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACa                                                                      >  1:1415631/1‑91 (MQ=255)
                         aTAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAat                                                                    <  1:2781794/91‑1 (MQ=255)
                            aTTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAattta                                                                 <  2:1168535/91‑1 (MQ=255)
                            aTTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAattta                                                                 >  1:569514/1‑91 (MQ=255)
                             ttttttGTACCAAAATGAAAATTTTTTTTAAAAATTAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatttat                                                                >  2:1768674/1‑91 (MQ=255)
                             ttttttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatttat                                                                <  2:3770660/91‑1 (MQ=255)
                             ttttttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatttat                                                                <  2:3841715/91‑1 (MQ=255)
                             ttttttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatttat                                                                <  2:569514/91‑1 (MQ=255)
                              tttttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatttatt                                                               <  1:2667789/91‑1 (MQ=255)
                              tttttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatttatt                                                               >  2:1241380/1‑91 (MQ=255)
                               ttttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatttattt                                                              >  1:3564662/1‑91 (MQ=255)
                                tttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAAtttattta                                                             >  2:2300579/1‑91 (MQ=255)
                                ttcgTACCAAAATGAAAATTTTTTTCAAAAACAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAAtttattta                                                             >  2:2961085/4‑91 (MQ=255)
                                 ttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaa                                                            >  2:3190031/1‑91 (MQ=255)
                                 ttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTACttaa                                                            >  2:1643488/1‑91 (MQ=255)
                                 ttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAACCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaa                                                            <  2:3487617/91‑1 (MQ=255)
                                  tGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaat                                                           >  2:2832030/1‑91 (MQ=255)
                                  tGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaat                                                           >  2:872105/1‑91 (MQ=255)
                                  tGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaat                                                           >  1:374143/1‑91 (MQ=255)
                                  tGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaat                                                           <  1:3636832/91‑1 (MQ=255)
                                  tGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTACttaat                                                           >  2:3145950/1‑91 (MQ=255)
                                      ccAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTtaattaat                                                       <  2:75601/91‑1 (MQ=255)
                                        aaaaTGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATgg                                                     <  1:1284652/91‑1 (MQ=255)
                                             gAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGaa                                                <  2:367244/91‑1 (MQ=255)
                                             gAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGaa                                                <  1:76967/91‑1 (MQ=255)
                                               aaaTTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATg                                              <  2:2929924/91‑1 (MQ=255)
                                                  tttttttCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGtat                                                                                 <  2:3618057/53‑1 (MQ=255)
                                                  tttttttCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCAt                                           <  2:3217178/91‑1 (MQ=255)
                                                     ttttCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGc                                        <  1:1121010/91‑1 (MQ=255)
                                                     ttttCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGc                                        <  1:2961047/91‑1 (MQ=255)
                                                        tCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAg                                     >  1:2107543/1‑91 (MQ=255)
                                                         cAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGa                                    >  1:2327107/1‑91 (MQ=255)
                                                          ttaaaTAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGaa                                   >  2:1552297/3‑91 (MQ=255)
                                                          aaaaaTAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGaa                                   <  1:3254142/91‑1 (MQ=255)
                                                          aaaaaTAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATACCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGaa                                   >  2:1726051/1‑91 (MQ=255)
                                                           aaaaTAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGaaa                                  <  1:2900414/91‑1 (MQ=255)
                                                            aaaTAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAAc                                 >  1:2489681/1‑91 (MQ=255)
                                                                aaaaaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCt                             >  1:438622/1‑91 (MQ=255)
                                                                 aaaaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTg                            <  1:806979/91‑1 (MQ=255)
                                                                 aaaaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTg                            >  1:1642269/1‑91 (MQ=255)
                                                                 aaaaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTg                            <  1:2568846/91‑1 (MQ=255)
                                                                 aaaaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTg                            >  2:2542461/1‑91 (MQ=255)
                                                                  tataCTAGCAATTCACATCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGa                           >  1:3481874/4‑91 (MQ=25)
                                                                   aaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGaa                          <  2:3828094/91‑1 (MQ=255)
                                                                    aaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAg                         >  2:78063/1‑91 (MQ=255)
                                                                    aaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAg                         >  2:215381/1‑91 (MQ=255)
                                                                     aCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGc                        <  2:686093/91‑1 (MQ=255)
                                                                     aCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGc                        <  2:383311/91‑1 (MQ=255)
                                                                     aCTAGCAATTCACATCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCGTAGCTAGAAACTTCTGAAGc                        >  1:1217072/1‑91 (MQ=25)
                                                                     aCTAGCAATTCACATCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGc                        >  2:3883249/1‑91 (MQ=37)
                                                                      cTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCt                       <  1:92859/91‑1 (MQ=255)
                                                                       tAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTg                      <  1:2267200/91‑1 (MQ=255)
                                                                         gCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGa                    <  2:1075321/91‑1 (MQ=255)
                                                                          cAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGaa                   <  2:3077907/91‑1 (MQ=255)
                                                                          cAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGaa                   >  2:2264495/1‑91 (MQ=255)
                                                                           aaTCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAAt                  >  1:3862193/1‑91 (MQ=255)
                                                                              ctcACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGt               <  2:3276194/91‑1 (MQ=255)
                                                                               tcACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTa              <  1:751275/91‑1 (MQ=255)
                                                                                cacaTGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTAc             <  2:1745078/91‑1 (MQ=255)
                                                                                    tgatgTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTa         >  2:1798359/1‑91 (MQ=255)
                                                                                      atgTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGt       <  2:2918287/91‑1 (MQ=255)
                                                                                       tgtgAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGtt      <  1:1703154/91‑1 (MQ=255)
                                                                                       tgtgAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGtt      <  1:2179020/91‑1 (MQ=255)
                                                                                        gtgAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGTTc     <  1:974005/91‑1 (MQ=255)
                                                                                         tgAATTGCTAGTATATATCAGTGCAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGTTCa    <  2:3132795/91‑1 (MQ=255)
                                                                                         tgAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGTTCa    >  1:2023270/1‑91 (MQ=255)
                                                                                          gAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGTTCat   <  1:1003414/91‑1 (MQ=255)
                                                                                          gAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGTTCat   <  1:1447766/91‑1 (MQ=255)
                                                                                           aaTTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGTTCata  <  1:99505/91‑1 (MQ=255)
                                                                                            |                                                                                          
AATATGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGTTCATA  >  USA300TCH1516_ALE/740237‑740418

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: