Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 581,901 G→C 6.4% V13L (GTG→CTG)  nusG_2 → Transcription termination/antitermination protein NusG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE581,9010GC 6.4% 699.7 / 23.8 203V13L (GTG→CTG) nusG_2Transcription termination/antitermination protein NusG
Reads supporting (aligned to +/- strand):  ref base G (93/97);  new base C (10/3);  total (103/100)
Fisher's exact test for biased strand distribution p-value = 8.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 211 total aligned reads displayed.

GATTTAGGAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCTTTAGAG  >  USA300TCH1516_ALE/581811‑581991
                                                                                          |                                                                                          
gATTTAGGAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAg                                                                                            <  2:3655703/91‑1 (MQ=255)
gATTTAGGAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAg                                                                                            >  2:2771695/1‑91 (MQ=255)
   ttAGGAATTACAGCATTGAAAAATTTATTATTGCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCACTGc                                                                                         <  2:1935933/91‑1 (MQ=255)
    tAGGAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGca                                                                                        >  2:229746/1‑91 (MQ=255)
     aGGAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGcat                                                                                       >  1:1274761/1‑91 (MQ=255)
     aGGAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGcat                                                                                       >  1:1276886/1‑91 (MQ=255)
     aGGAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGcat                                                                                       <  2:3342424/91‑1 (MQ=255)
     aGGAATTACAGCATCGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGcat                                                                                       <  2:1253767/91‑1 (MQ=255)
      ggAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGcata                                                                                      >  2:1675309/1‑91 (MQ=255)
      ggAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGcata                                                                                      >  1:2586356/1‑91 (MQ=255)
      ggAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGcata                                                                                      >  1:2838615/1‑91 (MQ=255)
      ggAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGcatac                                                                                     <  1:1338951/91‑1 (MQ=255)
      ggAATTACAGCAGTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGcata                                                                                      <  2:2003761/91‑1 (MQ=255)
       gAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGcatac                                                                                     <  2:3529785/91‑1 (MQ=255)
        aaTTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGcataca                                                                                    <  2:2222423/91‑1 (MQ=255)
         aTTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGcatacat                                                                                   <  1:3527243/91‑1 (MQ=255)
            aaaGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATAtt                                                                                <  2:891430/89‑1 (MQ=255)
            aCAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATAtt                                                                                >  1:842632/1‑91 (MQ=255)
            aCAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATAtt                                                                                >  1:1058094/1‑91 (MQ=255)
             cAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTc                                                                               >  1:3070702/1‑91 (MQ=255)
              aGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCt                                                                              <  1:1592206/91‑1 (MQ=255)
                cATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTgg                                                                            >  2:2588778/1‑91 (MQ=255)
                 aTTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGa                                                                           <  1:1123629/91‑1 (MQ=255)
                 aTTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGAAAAGCGTTGGTATGCAGCGCATACATATTCTGGa                                                                           >  2:3739310/1‑91 (MQ=255)
                  ttGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGat                                                                          >  1:1783822/1‑91 (MQ=255)
                    gAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGatat                                                                        <  1:1396677/91‑1 (MQ=255)
                    gAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGatat                                                                        <  1:1614042/91‑1 (MQ=255)
                    gAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTATGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGatat                                                                        <  1:2514603/91‑1 (MQ=255)
                     aaaaaTTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATg                                                                       >  2:2033198/1‑91 (MQ=255)
                     aaaaaTTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATg                                                                       <  2:2164103/91‑1 (MQ=255)
                     aaaaaTTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATg                                                                       >  2:3549205/1‑91 (MQ=255)
                      aaaaTTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAACGCGTTGGTATGCAGTGCATACATATTCTGGATATGa                                                                      >  2:3906276/1‑91 (MQ=255)
                      aaaaTTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGa                                                                      >  1:1420317/1‑91 (MQ=255)
                       aaaTTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGaa                                                                     <  1:321214/91‑1 (MQ=255)
                        aaTTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGaaa                                                                    >  2:660491/1‑91 (MQ=255)
                          tttATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAAt                                                                  <  2:857360/91‑1 (MQ=255)
                               tattTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGt                                                             <  1:2725166/91‑1 (MQ=255)
                                   tCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTaaa                                                         >  2:3453700/1‑91 (MQ=255)
                                     gTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTaaaaa                                                       >  2:2889570/1‑91 (MQ=255)
                                      ttAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAg                                                      <  1:3645991/91‑1 (MQ=255)
                                      ttAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAg                                                      <  1:3647697/91‑1 (MQ=255)
                                         gaggagTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAAt                                                   >  1:2885768/1‑91 (MQ=255)
                                         gaggagTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAAt                                                   >  1:2319042/1‑91 (MQ=255)
                                          aggagTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAAtt                                                  <  1:1710251/91‑1 (MQ=255)
                                          aggagTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAAtt                                                  >  2:1410914/1‑91 (MQ=255)
                                          aggagTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTGGGTAGGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAAtt                                                  <  2:2722036/91‑1 (MQ=255)
                                           ggagTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAAttt                                                 >  2:1431384/1‑91 (MQ=255)
                                            gagTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGGAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTa                                                <  1:3936698/91‑1 (MQ=255)
                                            gagTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTa                                                <  1:3223510/91‑1 (MQ=255)
                                            gagTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTa                                                <  1:3936523/91‑1 (MQ=255)
                                             agTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAg                                               <  2:678588/91‑1 (MQ=255)
                                             agTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAg                                               >  1:247678/1‑91 (MQ=255)
                                              gTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGa                                              >  1:2233606/1‑91 (MQ=255)
                                              gTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGa                                              >  2:3690391/1‑91 (MQ=255)
                                                  aGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGaaaaa                                          >  1:3344200/1‑91 (MQ=255)
                                                    aCATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAag                                        <  1:2629839/91‑1 (MQ=255)
                                                     cATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAaga                                       >  2:913520/1‑91 (MQ=255)
                                                     cATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAaga                                       <  2:375799/91‑1 (MQ=255)
                                                     cATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAaga                                       <  1:480448/91‑1 (MQ=255)
                                                       tGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGt                                     <  1:1459111/91‑1 (MQ=255)
                                                        gTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTa                                    >  2:2853214/1‑91 (MQ=255)
                                                         tCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAg                                   <  1:968348/91‑1 (MQ=255)
                                                         tCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAg                                   <  2:3245054/91‑1 (MQ=255)
                                                          cTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGa                                  >  1:504334/1‑91 (MQ=255)
                                                          cTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCGGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGa                                  <  1:3061796/91‑1 (MQ=255)
                                                           tGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTATGGATATGAAAATAAAGTTAAAAAAAAAAAAAAAAAAAGAGTAGaa                                 <  2:1114108/91‑1 (MQ=255)
                                                               gaagTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTa                             >  1:2349541/1‑91 (MQ=255)
                                                               gaagTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTa                             <  1:18022/91‑1 (MQ=255)
                                                               gaagTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTa                             <  1:476357/91‑1 (MQ=255)
                                                               gaagTTGGCGCAAAGCGTTGGTATGCAGTGCAGACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTa                             <  2:2398956/91‑1 (MQ=255)
                                                                 agTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATg                           <  1:1192863/91‑1 (MQ=255)
                                                                 agTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATg                           <  1:3453276/91‑1 (MQ=255)
                                                                 agTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAAAAAAAAAAAAAAAGAGTAGAATCTATg                           <  2:525350/91‑1 (MQ=255)
                                                                    tGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGCAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAat                        <  1:3140233/91‑1 (MQ=255)
                                                                    tGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAat                        >  1:2174640/1‑91 (MQ=255)
                                                                      gcgcAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAatat                      >  2:2405440/1‑91 (MQ=255)
                                                                      gcgcAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAatat                      >  1:67807/1‑91 (MQ=255)
                                                                       cgcAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATg                     <  1:1988687/91‑1 (MQ=255)
                                                                        gcAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGa                    >  2:3609186/1‑91 (MQ=255)
                                                                         cAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGAc                   <  1:3822058/91‑1 (MQ=255)
                                                                         cAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGAc                   <  1:310683/91‑1 (MQ=255)
                                                                           aaGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGCTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTg                 >  1:3141928/1‑91 (MQ=255)
                                                                            aGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGa                >  1:2763651/1‑91 (MQ=255)
                                                                              cGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAAc              >  2:2809120/1‑91 (MQ=255)
                                                                              cGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAAc              >  2:44240/1‑91 (MQ=255)
                                                                              cGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAAAAAAAAAAAAAGAGTAGAATCTATGAATATGACTGAAc              <  2:1590771/91‑1 (MQ=255)
                                                                               gTTGGAATGCAGAGAATAAATATTCTGGATATGAAAAAAAAGTAAAAAAAAAATAAGAAAAAAGAGTAGAATCTATGAATATGACTGAACa             <  2:852770/91‑1 (MQ=255)
                                                                                 tGGTATGCAGTGCATACATATTCTGGATATGCAAATAAAGTTAAACAGAATTCAGAAAAAAGAGTAGAATCTATGAATATTACTGAACaaa           >  2:760553/1‑91 (MQ=255)
                                                                                 atgTATGCACTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACaaa           >  1:3625650/3‑91 (MQ=255)
                                                                                  ggTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAAt          >  1:2966092/1‑91 (MQ=255)
                                                                                   gTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATc         >  2:1614253/1‑91 (MQ=255)
                                                                                    tATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCt        >  1:3751433/1‑91 (MQ=255)
                                                                                    tATGCACTGCATACATATTCTGGATATGAAAATAAAGTTAAAATGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCt        >  2:3059794/1‑91 (MQ=255)
                                                                                    tATGCACTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCt        >  2:3147181/1‑91 (MQ=255)
                                                                                    tATGCACTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCt        >  2:1648934/1‑91 (MQ=255)
                                                                                     aTGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCtt       <  1:63492/91‑1 (MQ=255)
                                                                                      tGCACTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCttt      >  2:2542229/1‑91 (MQ=255)
                                                                                      tGCACTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCttt      >  2:1592921/1‑91 (MQ=255)
                                                                                       gCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCTTTa     >  1:3289629/1‑91 (MQ=255)
                                                                                         aGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCTTTaga   >  2:3087004/1‑91 (MQ=255)
                                                                                          |                                                                                          
GATTTAGGAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCTTTAGAG  >  USA300TCH1516_ALE/581811‑581991

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: