Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,126,293 A→C 5.1% intergenic (+26/‑52) USA300HOU_RS10850 → / → USA300HOU_RS10855 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,126,2930AC 5.1% 914.4 / 29.5 276intergenic (+26/‑52)USA300HOU_RS10850/USA300HOU_RS10855hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (111/149);  new base C (2/12);  total (113/163)
Fisher's exact test for biased strand distribution p-value = 4.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 284 total aligned reads displayed.

AATACTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCATTAATAAATTAAA  >  USA300TCH1516_ALE/2126203‑2126382
                                                                                          |                                                                                         
aaTACTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACaa                                                                                           >  2:3007418/1‑91 (MQ=255)
aaTACTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACaa                                                                                           <  1:1545126/91‑1 (MQ=255)
 aTACTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAg                                                                                          <  2:172372/91‑1 (MQ=255)
  tACTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAgt                                                                                         <  2:1896726/91‑1 (MQ=255)
  tACTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAgt                                                                                         <  2:33681/91‑1 (MQ=255)
    cTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAACAACGTCTACAAgtgt                                                                                       >  2:3730334/1‑91 (MQ=255)
    cTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAgtgt                                                                                       <  2:872030/91‑1 (MQ=255)
    cTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAgtgt                                                                                       >  1:3568084/1‑91 (MQ=255)
      tGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACACTTGTAg                                                                                     <  1:269604/91‑1 (MQ=255)
       ggAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACACTTGTAGa                                                                                    <  1:3636590/91‑1 (MQ=255)
       ggAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGa                                                                                    <  1:2645429/91‑1 (MQ=255)
        gAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACACTTGTAGAc                                                                                   <  1:891157/91‑1 (MQ=255)
        gAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGAc                                                                                   >  2:473878/1‑91 (MQ=255)
         aagttatttCACAAAAAAGAGATGGCGATAAGTTAACCAAAGGTGATGCATTATGATTTTACTAATCAAAAAACGTCTACACTTGTAGACg                                                                                  <  2:3423184/88‑1 (MQ=255)
         aaTTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACg                                                                                  >  2:585545/1‑91 (MQ=255)
         aaTTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCAGTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACg                                                                                  <  2:2394855/91‑1 (MQ=255)
          atttatttCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGt                                                                                 <  1:530231/91‑1 (MQ=255)
            ttatttCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTg                                                                               <  2:3297810/91‑1 (MQ=255)
            ttatttCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTg                                                                               <  2:510183/91‑1 (MQ=255)
            ttatttCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTATAATAATAAAAAAACGTCTACAAGTGTAGACGTTg                                                                               <  2:382123/91‑1 (MQ=255)
             tatttCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGa                                                                              <  2:2309564/91‑1 (MQ=255)
              atttCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGaa                                                                             <  1:785602/91‑1 (MQ=255)
                ttCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAAtg                                                                           <  1:3031496/91‑1 (MQ=255)
                ttCAAAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAGACGTCTACAAGTGTAGACGTTGAAtg                                                                           <  2:1562896/91‑1 (MQ=255)
                tgcacaAAAAATAGATGGCGATAAGTTAACCAAAGAATATACATTATGATATTACTAATAAAAAAACGTCTACAAGTGTAGACGTTGAAtg                                                                           <  2:2329003/89‑1 (MQ=38)
                  cacaAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAAtggt                                                                         <  1:1125365/91‑1 (MQ=255)
                  cacaAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTAATAAAAAAAAAACGTCTACAAGTGTAGACGTTGAAtggt                                                                         <  2:150084/91‑1 (MQ=255)
                       aaaaTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTga                                                                    <  2:1259234/91‑1 (MQ=255)
                        aaaTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTgag                                                                   <  1:3442409/91‑1 (MQ=255)
                        aaaTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTgag                                                                   <  2:671114/91‑1 (MQ=255)
                         caTAGATGGCGATAAGTTAGCCGAAGAAGGTTCGTTATGCTTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTgaga                                                                  <  1:252062/90‑1 (MQ=25)
                            aGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAgtg                                                               >  1:656831/1‑91 (MQ=255)
                             gATGGCGATAAGTTAACTAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAgtgt                                                              >  2:1928197/1‑91 (MQ=255)
                             gATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAgtgt                                                              <  1:2887140/91‑1 (MQ=255)
                             gATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAgtgt                                                              >  1:2203450/1‑91 (MQ=255)
                             gATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAgtgt                                                              >  2:1929900/1‑91 (MQ=255)
                              aTGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGtgtg                                                             >  2:217295/1‑91 (MQ=255)
                               tGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGa                                                            >  1:1035043/1‑91 (MQ=255)
                                ggCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAg                                                           >  2:1006722/1‑91 (MQ=255)
                                 gCGATAAGTTAACCAAAGGAGATTCATTATGATTTTACGAATCAAAAAACGTCTACAAGTGTAGACGTTGAGTGGTGGTGAGAGTGTGAGc                                                          <  2:415863/91‑1 (MQ=255)
                                 gCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGc                                                          <  2:917032/91‑1 (MQ=255)
                                 gCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGc                                                          <  1:2105817/91‑1 (MQ=255)
                                 gCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGc                                                          >  1:3279473/1‑91 (MQ=255)
                                  cGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCg                                                         <  1:2961052/91‑1 (MQ=255)
                                   gATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCga                                                        >  1:1809052/1‑91 (MQ=255)
                                    aTAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCgag                                                       <  1:3842257/91‑1 (MQ=255)
                                     tAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCgaga                                                      >  2:2885089/1‑91 (MQ=255)
                                      aaGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGaa                                                     <  1:1584016/91‑1 (MQ=255)
                                       aGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAAt                                                    >  1:2562723/1‑91 (MQ=255)
                                        gTTAACCAAAGAAGATTCATTATGATTTTACGAATCAAAAAACGGCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATa                                                   <  1:3154933/91‑1 (MQ=255)
                                         ttAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATaa                                                  <  1:2372364/91‑1 (MQ=255)
                                         ttAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATaa                                                  >  1:2757865/1‑91 (MQ=255)
                                          gaaCCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATaaa                                                 >  1:1356574/2‑91 (MQ=255)
                                           aaCCAAAGAAGATTCATATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAgg                                               <  2:2752283/91‑1 (MQ=255)
                                               aaaGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAATTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGgaga                                            <  2:3544724/91‑1 (MQ=255)
                                                aagaagATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGaa                                           <  1:203338/91‑1 (MQ=255)
                                                 agaagaTTCATTATGATTTTACTAATCAACAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGaaa                                          >  2:1106284/1‑91 (MQ=255)
                                                 agaagaTTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGaaa                                          >  2:3073552/1‑91 (MQ=255)
                                                  gaagaTTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAat                                         >  2:3555475/1‑91 (MQ=255)
                                                   aagaTTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAatg                                        >  2:2530013/1‑91 (MQ=255)
                                                    agaTTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAatga                                       >  2:3219211/1‑91 (MQ=255)
                                                    agaTTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAatga                                       <  1:2414899/91‑1 (MQ=255)
                                                     gaTTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAatgat                                      <  1:1390034/91‑1 (MQ=255)
                                                     gaTTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAatgat                                      <  1:1866978/91‑1 (MQ=255)
                                                       ttCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAAtgatga                                    >  2:1208312/1‑91 (MQ=255)
                                                       ttCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAAtgatga                                    <  1:3353049/91‑1 (MQ=255)
                                                        tCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGAc                                   <  1:95017/91‑1 (MQ=255)
                                                         cATTATGATTTTACTAATCAAAAAACGTCTATAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACg                                  >  2:1372391/1‑91 (MQ=255)
                                                         cATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACg                                  >  1:3696076/1‑91 (MQ=255)
                                                          aTTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGc                                 >  1:14366/1‑91 (MQ=255)
                                                           ttATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCa                                <  1:1267691/91‑1 (MQ=255)
                                                               gATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCataat                            >  1:3761289/1‑91 (MQ=255)
                                                                 ttttCCTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAatat                          >  1:3125429/1‑91 (MQ=255)
                                                                  tttACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAAtata                         <  2:1856772/91‑1 (MQ=255)
                                                                   ttACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAg                        <  2:960565/91‑1 (MQ=255)
                                                                      cTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGaaa                     <  1:2719651/91‑1 (MQ=255)
                                                                        aaTCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGaaaaa                   >  2:1405187/1‑91 (MQ=255)
                                                                            aaaaaaCGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCa               >  1:3560695/1‑91 (MQ=255)
                                                                             aaaaaCGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCAt              >  1:3507043/1‑91 (MQ=255)
                                                                                aaCGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCATtaa           <  1:922872/91‑1 (MQ=255)
                                                                                  cGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCATtaata         <  1:3894599/91‑1 (MQ=255)
                                                                                  cGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCATtaata         >  1:1477157/1‑91 (MQ=255)
                                                                                  cGTCTACAAGTCTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCATtaata         <  1:3590444/91‑1 (MQ=255)
                                                                                   gTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCATtaataa        <  2:3690912/91‑1 (MQ=255)
                                                                                    tCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCATTAATaaa       <  1:1503820/91‑1 (MQ=255)
                                                                                    tCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCATTAATaaa       >  1:1229476/1‑91 (MQ=255)
                                                                                     cTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCATTAATAAAt      >  2:2082644/1‑91 (MQ=255)
                                                                                     cTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCATTAATAAAt      >  1:730910/1‑91 (MQ=255)
                                                                                     cTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCATTAATAAAt      >  1:1066996/1‑91 (MQ=255)
                                                                                     cTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCATTAATAAAt      >  1:1003143/1‑91 (MQ=255)
                                                                                     cTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCATTAATAAAt      >  1:3480834/1‑91 (MQ=255)
                                                                                      tACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGCATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCATTAATAAAtt     >  2:2423272/1‑91 (MQ=255)
                                                                                      tACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCATTAATAAAtt     >  1:681593/1‑91 (MQ=255)
                                                                                      tACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCATTAATAAAtt     >  1:1704615/1‑91 (MQ=255)
                                                                                       aCAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCATTAATAAATTa    <  1:2176914/91‑1 (MQ=255)
                                                                                       aCAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCATTAATAAATTa    >  2:2060697/1‑91 (MQ=255)
                                                                                       aCAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCATTAATAAATTa    <  2:1671761/91‑1 (MQ=255)
                                                                                       aCAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCATTAATAAATTa    <  1:871133/91‑1 (MQ=255)
                                                                                        cAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCATTAATAAATTaa   >  2:2828845/1‑91 (MQ=255)
                                                                                         aaGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCATTAATAAATTaaa  <  2:1290000/91‑1 (MQ=255)
                                                                                          |                                                                                         
AATACTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCGAGAATAAAGGAGAAATGATGACGCATAATATAGAAAAACGCATTAATAAATTAAA  >  USA300TCH1516_ALE/2126203‑2126382

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: