Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 581,901 G→C 6.9% 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.9% 655.2 / 26.6 188V13L (GTG→CTG) nusG_2Transcription termination/antitermination protein NusG
Reads supporting (aligned to +/- strand):  ref base G (100/75);  new base C (10/3);  total (110/78)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 193 total aligned reads displayed.

GATTTAGGAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCTTTAGAG  >  USA300TCH1516_ALE/581811‑581991
                                                                                          |                                                                                          
gATTTAGGAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAg                                                                                            >  1:1637842/1‑91 (MQ=255)
gATTTAGGAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAg                                                                                            >  2:887903/1‑91 (MQ=255)
  tttAGGAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTg                                                                                          >  2:718217/1‑91 (MQ=255)
   ttAGGAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCACTGc                                                                                         <  1:464175/91‑1 (MQ=255)
     aGGAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGcat                                                                                       <  1:137298/91‑1 (MQ=255)
     aGGAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGcat                                                                                       <  1:764663/91‑1 (MQ=255)
     aGGAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGcat                                                                                       <  1:764544/91‑1 (MQ=255)
      ggAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGcata                                                                                      <  1:437669/91‑1 (MQ=255)
      ggAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGcata                                                                                      <  1:2143040/91‑1 (MQ=255)
      ggAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGcata                                                                                      <  1:2141935/91‑1 (MQ=255)
          ttACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGcatacata                                                                                  <  2:3477781/91‑1 (MQ=255)
           tACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACatat                                                                                 >  1:3706410/1‑91 (MQ=255)
           tACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACatat                                                                                 >  2:3567680/1‑91 (MQ=255)
             cAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTc                                                                               >  1:2182562/1‑91 (MQ=255)
             cAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTc                                                                               >  1:3290760/1‑91 (MQ=255)
               gCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTg                                                                             >  2:160359/1‑91 (MQ=255)
                    gAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGatat                                                                        <  2:2301233/91‑1 (MQ=255)
                     aaaaaTTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATg                                                                       <  1:2307758/91‑1 (MQ=255)
                      aaaaTTTATTATTTGGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGa                                                                      <  1:1401888/91‑1 (MQ=255)
                      aaaaTTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGa                                                                      >  2:3584096/1‑91 (MQ=255)
                      aaaaTTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGa                                                                      >  2:2599755/1‑91 (MQ=255)
                      aaaaTTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGa                                                                      >  2:2043394/1‑91 (MQ=255)
                      aaaaTTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGa                                                                      >  2:2535766/1‑91 (MQ=255)
                         aTTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGaaaa                                                                   >  1:2632479/1‑91 (MQ=255)
                          tttATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAAt                                                                  >  1:3739175/1‑91 (MQ=255)
                              ttattTCGTTAGAGGAGTGAAGACATGTCTGCAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAg                                                              >  1:1895001/1‑91 (MQ=255)
                              ttattTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAg                                                              >  1:2885506/1‑91 (MQ=255)
                                attTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGtt                                                            <  2:201716/91‑1 (MQ=255)
                                   tCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTaaa                                                         <  1:3036615/91‑1 (MQ=255)
                                   tCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTaaa                                                         <  2:3204838/91‑1 (MQ=255)
                                   tCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTaaa                                                         <  2:1545629/91‑1 (MQ=255)
                                    cGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTaaaa                                                        <  1:1166737/91‑1 (MQ=255)
                                      ttAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAg                                                      <  2:1274757/91‑1 (MQ=255)
                                        agagGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGaa                                                    <  2:3770273/91‑1 (MQ=255)
                                          aggagTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAAtt                                                  >  1:1996346/1‑91 (MQ=255)
                                          aagagTGAAGACATGTCTGAAGAAGTTGTCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAAtt                                                  <  2:1273082/89‑1 (MQ=255)
                                            gagTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTa                                                <  1:2963878/91‑1 (MQ=255)
                                            gagTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTa                                                <  2:1860668/91‑1 (MQ=255)
                                              gTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGa                                              >  2:300063/1‑91 (MQ=255)
                                              gTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGa                                              <  1:250668/91‑1 (MQ=255)
                                                gAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGaaa                                            >  2:47011/1‑91 (MQ=255)
                                                 aaGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGaaaa                                           >  1:3762461/1‑91 (MQ=255)
                                                 aaGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGaaaa                                           <  1:1578023/91‑1 (MQ=255)
                                                   gACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGaaaaaa                                         <  1:1531805/91‑1 (MQ=255)
                                                    aCATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAag                                        <  2:2923542/91‑1 (MQ=255)
                                                    aCATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAag                                        <  1:893892/91‑1 (MQ=255)
                                                       tGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAGGTTAAAAAAAAAAAAAAAAAAAGAGt                                     <  2:3739175/91‑1 (MQ=255)
                                                       tGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGt                                     <  2:3290760/91‑1 (MQ=255)
                                                       tGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGt                                     <  1:2535766/91‑1 (MQ=255)
                                                        gTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTa                                    <  2:3786160/91‑1 (MQ=255)
                                                        gTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTa                                    >  1:2803647/1‑91 (MQ=255)
                                                        gTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTa                                    >  2:2501859/1‑91 (MQ=255)
                                                         tCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAg                                   >  1:2302669/1‑91 (MQ=255)
                                                          cTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGa                                  >  1:2173170/1‑91 (MQ=255)
                                                          cTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGa                                  >  2:760798/1‑91 (MQ=255)
                                                           tGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGaa                                 <  2:3173049/91‑1 (MQ=255)
                                                            gaagaaGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAAt                                <  1:2386512/91‑1 (MQ=255)
                                                             aagaagTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGt                                     >  2:1156981/1‑85 (MQ=255)
                                                             aagaagTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATc                               >  1:679703/1‑91 (MQ=255)
                                                             aagaagTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATc                               >  1:3012936/1‑91 (MQ=255)
                                                               gaagTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTa                             <  2:2729985/91‑1 (MQ=255)
                                                                aagTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTAt                            >  2:3315849/1‑91 (MQ=255)
                                                                 agTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATg                           >  2:3677245/1‑91 (MQ=255)
                                                                   ttGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGaa                         >  2:1985361/1‑91 (MQ=255)
                                                                   ttGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAAAAAAAAAAAAAAAGAGTAGAATCTATGaa                         <  2:1220489/91‑1 (MQ=255)
                                                                     ggCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAata                       >  2:900818/1‑91 (MQ=255)
                                                                      gcgcAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAatat                      >  1:1141897/1‑91 (MQ=255)
                                                                      gcgcAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAatat                      >  1:2289457/1‑91 (MQ=255)
                                                                      gcgcAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAatat                      >  1:3582304/1‑91 (MQ=255)
                                                                        gcAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGa                    >  2:2350613/1‑91 (MQ=255)
                                                                           aaGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTg                 >  1:274134/1‑91 (MQ=255)
                                                                            aGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGa                >  1:1013248/1‑91 (MQ=255)
                                                                            aGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGa                <  2:2688381/91‑1 (MQ=255)
                                                                            aGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGa                >  2:2712621/1‑91 (MQ=255)
                                                                            aGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGa                >  2:577055/1‑91 (MQ=255)
                                                                            aGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGa                >  2:1232584/1‑91 (MQ=255)
                                                                              cGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAAc              <  2:2182562/91‑1 (MQ=255)
                                                                              cGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAAc              <  1:3186729/91‑1 (MQ=255)
                                                                              cGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAAc              <  2:1831686/91‑1 (MQ=255)
                                                                              cGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAAc              >  2:1542537/1‑91 (MQ=255)
                                                                               gTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACa             >  2:3113389/1‑91 (MQ=255)
                                                                               atatgTATGCACTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACa             >  1:2810152/5‑91 (MQ=255)
                                                                                 tGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACaaa           >  2:79619/1‑91 (MQ=255)
                                                                                 tGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACaaa           >  2:1039950/1‑91 (MQ=255)
                                                                                 tGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACaaa           >  1:3479627/1‑91 (MQ=255)
                                                                                 atgTATGCACTGCATACCTATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACaaa           >  1:171529/3‑91 (MQ=255)
                                                                                  tgTATGCACTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAAt          >  2:3448648/2‑91 (MQ=255)
                                                                                  ggTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGTATATGACTGAACAAAt          >  2:1970844/1‑91 (MQ=255)
                                                                                  ggTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAAt          >  2:3260306/1‑91 (MQ=255)
                                                                                  ggTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAAt          >  2:1122075/1‑91 (MQ=255)
                                                                                    tATGCACTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCt        >  1:3383052/1‑91 (MQ=255)
                                                                                     aTGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTg                 <  2:146157/81‑1 (MQ=255)
                                                                                     aTGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCtt       >  1:2907697/1‑91 (MQ=255)
                                                                                     aTGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCtt       >  2:1110161/1‑91 (MQ=255)
                                                                                     aTGCACTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTCGAAAAAAGAGTAGAATCTATGAATATGACAGAACAAATCtt       >  2:3406997/1‑91 (MQ=255)
                                                                                     aTGCACTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCtt       >  2:322900/1‑91 (MQ=255)
                                                                                      tGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGATAAAAGAGTAGAATCTATGAATATGACTGAACAAATCttt      >  2:831807/1‑91 (MQ=255)
                                                                                      tGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCttt      <  2:51347/91‑1 (MQ=255)
                                                                                       gCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCTTTa     >  2:1715243/1‑91 (MQ=255)
                                                                                          gTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCTTTagag  <  1:2046530/91‑1 (MQ=255)
                                                                                          |                                                                                          
GATTTAGGAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCTTTAGAG  >  USA300TCH1516_ALE/581811‑581991

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: