Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 739,250 T→A 7.8% I162N (ATC→AAC)  pitA_1 → Low‑affinity inorganic phosphate transporter 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE739,2500TA 7.8% 862.4 / 44.7 266I162N (ATC→AAC) pitA_1Low‑affinity inorganic phosphate transporter 1
Reads supporting (aligned to +/- strand):  ref base T (136/108);  new base A (4/17);  total (141/125)
Fisher's exact test for biased strand distribution p-value = 1.27e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 267 total aligned reads displayed.

AAGGTTTCACAAAAATTATTATTGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATCATTCTCTCACGGTA  >  USA300TCH1516_ALE/739160‑739339
                                                                                          |                                                                                         
 aGGTTTCACAAAAATTATTATTGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATc                                                                                          >  2:2823590/1‑91 (MQ=255)
 aGGTTTCACAAAAATTATTATTGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATc                                                                                          >  2:1994173/1‑91 (MQ=255)
  ggTTTCACAAAAATTATTATTGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCt                                                                                         >  2:2234653/1‑91 (MQ=255)
   gTTTCACAAAAATTATTATTGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCtt                                                                                        <  1:2807039/91‑1 (MQ=255)
    tttCACAAAAATTATTATTGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGATAACttt                                                                                       <  2:208895/91‑1 (MQ=255)
     ttCACAAAAATTATTATTGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCACTTTTTAAAGTTATCTTTa                                                                                      >  2:2045971/1‑91 (MQ=255)
     ttCACAAAAATTATTATTGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTa                                                                                      >  2:2045739/1‑91 (MQ=255)
     ttCACAAAAATTATTATTGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGATAACTTTa                                                                                      <  2:3690330/91‑1 (MQ=255)
      tCACAAAAATTATTATTGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCttt                                                                                       >  1:438932/1‑89 (MQ=255)
      tCACAAAAATTATTATTGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTaa                                                                                     <  2:1549330/91‑1 (MQ=255)
      tCACAAAAATTATTATTGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGATAACTTTaa                                                                                     <  2:2939582/91‑1 (MQ=255)
       cacaAAAATTATTATTGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTaaa                                                                                    >  1:2166684/1‑91 (MQ=255)
       cacaAAAATTATTATTGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTaaa                                                                                    >  1:2988596/1‑91 (MQ=255)
       cacaAAAATTATTATTGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGATAACTTTaaa                                                                                    <  1:2831588/91‑1 (MQ=255)
        acaAAAATTATTATTGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTaaaa                                                                                   >  2:728646/1‑91 (MQ=255)
        acaAAAATTATTATTGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGATAACTTTaaaa                                                                                   <  1:1061393/91‑1 (MQ=255)
          aaaaaTTATTATTGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGATAACTTTAAAAAt                                                                                 <  1:1968717/91‑1 (MQ=255)
            aaaTTATTATTGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGc                                                                               >  1:1352490/1‑91 (MQ=255)
            aaaTTATTATTGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGATAACTTTAAAAAttg                                                                               <  1:1606796/91‑3 (MQ=39)
              attattattGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCaa                                                                             <  1:3350933/91‑1 (MQ=255)
                tattattGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGc                                                                               >  2:268590/1‑87 (MQ=255)
                 attattGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAAtt                                                                          <  1:2802269/91‑1 (MQ=255)
                 attattGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAAtt                                                                          >  1:3514349/1‑91 (MQ=255)
                 attattGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAAtt                                                                          >  1:745402/1‑91 (MQ=255)
                        tATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAag                                                                   >  1:374235/1‑91 (MQ=255)
                           tAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGc                                                                <  2:2336871/91‑1 (MQ=255)
                           tAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGc                                                                <  1:2455193/91‑1 (MQ=255)
                           tAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGc                                                                <  1:522792/91‑1 (MQ=255)
                              tCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGaa                                                             >  1:1256527/1‑91 (MQ=255)
                              tCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGaa                                                             <  1:2851067/91‑1 (MQ=255)
                               cGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAAt                                                            >  1:3199372/1‑91 (MQ=255)
                                gTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATc                                                           >  2:3623601/1‑91 (MQ=255)
                                 tttCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCg                                                          >  2:2037291/1‑91 (MQ=255)
                                   tCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTa                                                        >  2:1247907/1‑91 (MQ=255)
                                     aCCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAAc                                                      >  2:1099266/1‑91 (MQ=255)
                                     aCCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAAc                                                      <  1:317612/91‑1 (MQ=255)
                                      ccGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACt                                                     >  1:1030551/1‑91 (MQ=255)
                                       cGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACtt                                                    <  2:1840257/91‑1 (MQ=255)
                                       cGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACtt                                                    <  2:915115/91‑1 (MQ=255)
                                           tATCGCATTTTGCGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGa                                                >  1:2981002/1‑91 (MQ=255)
                                             tCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGAtt                                              >  2:69945/1‑91 (MQ=255)
                                              cGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGAttt                                             <  1:1054941/91‑1 (MQ=255)
                                              cGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGAttt                                             >  1:1242050/1‑91 (MQ=255)
                                              cGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGAttt                                             <  1:2609798/91‑1 (MQ=255)
                                              cGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGCAACTTTAGAttt                                             >  2:1638893/1‑91 (MQ=255)
                                               gCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGAtttt                                            <  1:823719/91‑1 (MQ=255)
                                               gCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGAtttt                                            >  2:3351637/1‑91 (MQ=255)
                                                 aTTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTc                                          <  2:834762/91‑1 (MQ=255)
                                                  ttttGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTcc                                         >  2:3668534/1‑91 (MQ=255)
                                                  ttttGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTcc                                         >  2:3047339/1‑91 (MQ=255)
                                                  ttttGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTcc                                         <  2:1344682/91‑1 (MQ=255)
                                                  ttttGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACATTAGATTTTTcc                                         >  2:3231765/1‑91 (MQ=255)
                                                      gtgtTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAAt                                     >  2:3149570/1‑91 (MQ=255)
                                                      gtgtTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAAt                                     >  2:2551041/1‑91 (MQ=255)
                                                        gtTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTCAGATTTTTCCAAAttt                                   >  2:3703840/1‑91 (MQ=255)
                                                         ttGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAAtttt                                  >  1:6594/1‑91 (MQ=255)
                                                         ttGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAAtttt                                  >  2:279537/1‑91 (MQ=255)
                                                         ttGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAAtttt                                  >  1:2480050/1‑91 (MQ=255)
                                                         ttGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTCTCCAAAttct                                  >  2:1595999/1‑89 (MQ=255)
                                                           ggTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCACATTTAACAAGAGCGAATCGTAACTTTAGATTCTTCCAAATATTcc                                >  2:3155885/1‑90 (MQ=255)
                                                           ggTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTca                                <  2:45680/91‑1 (MQ=255)
                                                             tttCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTcaca                              <  2:3318493/91‑1 (MQ=255)
                                                             tttCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTcaca                              <  1:2699054/91‑1 (MQ=255)
                                                             tttCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTcaca                              <  1:572753/91‑1 (MQ=255)
                                                             tttCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTcaca                              <  2:3394503/91‑1 (MQ=255)
                                                               tCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCAcagc                            >  1:3268795/1‑91 (MQ=255)
                                                                cTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCAcagca                           >  1:173877/1‑91 (MQ=255)
                                                                  tGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACagcagc                         <  2:1536527/91‑1 (MQ=255)
                                                                  tGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACagcagc                         >  1:1780198/1‑91 (MQ=255)
                                                                   gATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGAttttt                                           <  2:1410941/72‑1 (MQ=255)
                                                                    aTGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGt                       >  2:449881/1‑91 (MQ=255)
                                                                    aTGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGt                       >  1:193693/1‑91 (MQ=255)
                                                                     tGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGtt                      <  1:2128877/91‑1 (MQ=255)
                                                                     tGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGtt                      >  2:3515701/1‑91 (MQ=255)
                                                                      gTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTa                     <  1:1548502/91‑1 (MQ=255)
                                                                         ttCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACaa                  >  1:950723/1‑91 (MQ=255)
                                                                           cAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATc                <  1:2761486/91‑1 (MQ=255)
                                                                           cAATTTTTAAAGTTATCTTTAAAAATGCAAATATAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATc                >  2:3801137/1‑91 (MQ=255)
                                                                            aaTTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATCa               >  2:2293940/1‑91 (MQ=255)
                                                                             aTTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATCAt              >  1:3065358/1‑91 (MQ=255)
                                                                               ttttAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATCATtc            <  1:892665/91‑1 (MQ=255)
                                                                                tttAAAGATAACTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATCATtct           >  1:2923647/1‑91 (MQ=255)
                                                                                 ttAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATCATtctc          <  1:1801437/91‑1 (MQ=255)
                                                                                 ttAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATCATtctc          <  1:3759271/91‑1 (MQ=255)
                                                                                 ttAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATCATtctc          <  1:2874165/91‑1 (MQ=255)
                                                                                  tAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATCATtctct         >  1:1837021/1‑91 (MQ=255)
                                                                                   aaaGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATCATtctctc        <  1:3465894/91‑1 (MQ=255)
                                                                                   aaaGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATCATtctctc        >  1:1443652/1‑91 (MQ=255)
                                                                                   aaaGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATCATtctctc        >  1:1462407/1‑91 (MQ=255)
                                                                                    aaGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATCATTCTCTCa       <  1:3154684/91‑1 (MQ=255)
                                                                                     aGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATCATTCTCTCAc      <  1:2234653/91‑1 (MQ=255)
                                                                                     aGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATCATTCTCTCAc      >  2:2474870/1‑91 (MQ=255)
                                                                                     aGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATCATTCTCTCAc      >  1:2346433/1‑91 (MQ=255)
                                                                                     aGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATCATTCTCTCAc      >  1:2058796/1‑91 (MQ=255)
                                                                                       ttATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATCATTCTCTCACgg    <  2:275688/91‑1 (MQ=255)
                                                                                        tATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATCATTCTCTCACGGt   <  2:477231/91‑1 (MQ=255)
                                                                                        tATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATCATTCTCTCACGGt   <  1:1432411/91‑1 (MQ=255)
                                                                                         aTCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATCATTCTCTCACGGTa  >  2:289407/1‑91 (MQ=255)
                                                                                         aTCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATCATTCTCTCACGGTa  <  1:3706885/91‑1 (MQ=255)
                                                                                         aTCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATCATTCTCTCACGGTa  <  1:2556622/91‑1 (MQ=255)
                                                                                          |                                                                                         
AAGGTTTCACAAAAATTATTATTGTATTAATCGTTTCACCGATTATCGCATTTTGTGTTGGTTTCTTGATGTATTCAATTTTTAAAGTTATCTTTAAAAATGCAAATTTAACAAGAGCGAATCGTAACTTTAGATTTTTCCAAATTTTCACAGCAGCGTTACAATCATTCTCTCACGGTA  >  USA300TCH1516_ALE/739160‑739339

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: