Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,447,045 T→A 6.1% I64F (ATC→TTC)  cobT ← Aerobic cobaltochelatase subunit CobT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,447,0450TA 6.1% 779.5 / 30.1 230I64F (ATC→TTC) cobTAerobic cobaltochelatase subunit CobT
Reads supporting (aligned to +/- strand):  ref base T (109/107);  new base A (11/3);  total (120/110)
Fisher's exact test for biased strand distribution p-value = 5.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

Only 100 of 238 total aligned reads displayed.

CTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTAACTTGTGT  >  USA300TCH1516_ALE/1446955‑1447135
                                                                                          |                                                                                          
cTTGTACTACTCCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGat                                                                                            <  2:2593164/91‑1 (MQ=255)
 ttGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGata                                                                                           <  1:2514016/91‑1 (MQ=255)
  tGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATATCGaaat                                                                                          <  1:2886531/91‑1 (MQ=255)
   gTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATg                                                                                         >  1:408253/1‑91 (MQ=255)
    tactacTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATgg                                                                                        <  1:3540031/91‑1 (MQ=255)
     actactTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGc                                                                                       <  1:111659/91‑1 (MQ=255)
     actactTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGc                                                                                       <  2:64702/91‑1 (MQ=255)
     actactTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATATCGAAATGGc                                                                                       <  1:69033/91‑1 (MQ=255)
       tactTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCtt                                                                                     <  2:10330/91‑1 (MQ=255)
       tactTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCtt                                                                                     <  2:1430392/91‑1 (MQ=255)
       tactTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCtt                                                                                     <  2:2613250/91‑1 (MQ=255)
        actTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTg                                                                                    <  1:2760317/91‑1 (MQ=255)
        actTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTg                                                                                    >  2:229047/1‑91 (MQ=255)
          ttCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGt                                                                                  <  2:204312/91‑1 (MQ=255)
           tCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATCTCGATATGGCTTGGTc                                                                                 >  2:2983091/1‑91 (MQ=255)
             aTTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGa                                                                               <  1:3212613/91‑1 (MQ=255)
              tttACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTcaca                                                                              <  1:792741/91‑4 (MQ=255)
                tACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGAtgt                                                                            <  1:1620920/91‑1 (MQ=255)
                 aCAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGAtgtg                                                                           <  2:831869/91‑1 (MQ=255)
                 aCAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGAtgtg                                                                           <  2:378727/91‑1 (MQ=255)
                  cAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGa                                                                          >  2:3694684/1‑91 (MQ=255)
                   aaTTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGAc                                                                         >  2:1464705/1‑91 (MQ=255)
                    aTTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGAcc                                                                        >  2:1130185/1‑91 (MQ=255)
                     tttGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCGGCCATTTCGATATGGCTTGGTCGATGTGACCa                                                                       <  2:3810666/91‑1 (MQ=255)
                     tttGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCa                                                                       <  1:3694384/91‑1 (MQ=255)
                      ttGTACATCCATCATATGATAACCGTAAGGAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCaa                                                                      <  1:3114696/91‑1 (MQ=255)
                      ttGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCaa                                                                      <  2:2097129/91‑1 (MQ=255)
                       tGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCaaa                                                                     >  2:3547475/1‑91 (MQ=255)
                        gTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCaaaa                                                                    <  1:1607183/91‑1 (MQ=255)
                         tACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAAc                                                                   >  1:1222799/1‑91 (MQ=255)
                         tACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAAc                                                                   <  2:1201751/91‑1 (MQ=255)
                          aCATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAAcc                                                                  >  2:3652413/1‑91 (MQ=255)
                           catccatcATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCa                                                                 >  2:2105587/1‑91 (MQ=255)
                            atccatcATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCaa                                                                >  2:1709239/1‑91 (MQ=255)
                                atcatATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCg                                                            <  1:1710477/91‑1 (MQ=255)
                                 tcatATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGt                                                           >  1:2167527/1‑91 (MQ=255)
                                      tGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCa                                                      <  1:3846697/91‑1 (MQ=255)
                                      tGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCa                                                      <  1:3846888/91‑1 (MQ=255)
                                       gATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCaa                                                     <  2:2923492/91‑1 (MQ=255)
                                        aTAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCaaa                                                    >  2:1231390/1‑91 (MQ=255)
                                         tAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCaaaa                                                   >  1:516542/1‑91 (MQ=255)
                                          aaCCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAAc                                                  >  2:2065192/1‑91 (MQ=255)
                                            ccGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACAt                                                <  2:1421525/91‑1 (MQ=255)
                                            ccGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACAt                                                <  1:1102295/91‑1 (MQ=255)
                                              gTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACAttt                                              <  1:1454780/91‑1 (MQ=255)
                                              gTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACAttt                                              <  1:3359447/91‑1 (MQ=255)
                                              gTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACAtt                                               >  2:2451327/1‑90 (MQ=255)
                                                aaGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTg                                            <  1:194824/91‑1 (MQ=255)
                                                 agcagcTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGa                                           <  1:2000580/91‑1 (MQ=255)
                                                 agcagcTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGa                                           >  2:2820252/1‑91 (MQ=255)
                                                 agcagcTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGa                                           >  2:3381451/1‑91 (MQ=255)
                                                 agcagcTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGa                                           <  2:1545836/91‑1 (MQ=255)
                                                  gcagcTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGaa                                          >  2:3152747/1‑91 (MQ=255)
                                                    agcTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACa                                        <  1:3783216/91‑1 (MQ=255)
                                                     gcTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAg                                       >  1:2059779/1‑91 (MQ=255)
                                                     gcTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAg                                       >  2:3096836/1‑91 (MQ=255)
                                                      cTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAgg                                      >  1:1411478/1‑91 (MQ=255)
                                                       taacatatcaTCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGa                                     <  1:2445567/91‑1 (MQ=255)
                                                       taacataacaTCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGa                                     <  2:3282115/91‑1 (MQ=255)
                                                       taacataacaTCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGa                                     >  1:3632285/1‑91 (MQ=255)
                                                        aacataacatCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGAt                                    <  1:1801852/91‑1 (MQ=255)
                                                        aacataacatCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGAt                                    >  1:2998650/1‑91 (MQ=255)
                                                         acataacatCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGAtt                                   <  2:1978857/91‑1 (MQ=255)
                                                           ataacatCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACAtt                                               <  1:1556513/77‑1 (MQ=255)
                                                           ataacatCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATtat                                 >  2:3781230/1‑91 (MQ=255)
                                                           ataacatCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATtat                                 >  2:5561/1‑91 (MQ=255)
                                                            taacatCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATtata                                >  2:240648/1‑91 (MQ=255)
                                                            taacatCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATtata                                >  1:2116223/1‑91 (MQ=255)
                                                               catCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTa                             >  1:1256171/1‑91 (MQ=255)
                                                                atCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGGaa                            <  1:739875/91‑1 (MQ=255)
                                                                 tCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAg                           <  2:2746333/91‑1 (MQ=255)
                                                                  cTGTTTTTAAACCAGCCATTTCGATATTGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAgg                          <  2:1618278/91‑1 (MQ=255)
                                                                    gTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGa                        <  2:516542/91‑1 (MQ=255)
                                                                     tttttAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGaa                       >  2:582594/1‑91 (MQ=255)
                                                                      ttttAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAAt                      <  1:837500/91‑1 (MQ=255)
                                                                        ttAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAAttt                    <  2:1950288/91‑1 (MQ=255)
                                                                        ttAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAAttt                    <  2:3216429/91‑1 (MQ=255)
                                                                           aaCCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTg                 <  2:46195/91‑1 (MQ=255)
                                                                            ccaaGCCATATCGAAATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTgg                >  1:935706/4‑91 (MQ=255)
                                                                            aCCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTgg                <  2:1827112/91‑1 (MQ=255)
                                                                            aCCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTgg                >  1:2747957/1‑91 (MQ=255)
                                                                             ccAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGa               >  2:3706976/1‑91 (MQ=255)
                                                                              cAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAAGCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTCTGGAt              >  1:3612747/1‑91 (MQ=255)
                                                                              cAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGAt              >  1:886723/1‑91 (MQ=255)
                                                                              cAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGAt              >  1:2274701/1‑91 (MQ=255)
                                                                                gCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGAttt            >  1:3390498/1‑91 (MQ=255)
                                                                                gCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGAttt            >  1:3176562/1‑91 (MQ=255)
                                                                                gCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGAttt            >  2:1768854/1‑91 (MQ=255)
                                                                                 ccATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGAtttt           >  1:2247556/1‑91 (MQ=255)
                                                                                 ccATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGAtttt           <  1:397287/91‑1 (MQ=255)
                                                                                 ccATATCGAAATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGAtttt           >  1:3529899/1‑91 (MQ=255)
                                                                                  cATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTa          >  1:2745425/1‑91 (MQ=255)
                                                                                  cATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTa          >  2:2206841/1‑91 (MQ=255)
                                                                                  cATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTa          >  1:1234857/1‑91 (MQ=255)
                                                                                    tttCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTAAc        >  2:460417/1‑91 (MQ=255)
                                                                                    tatCGAAATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTAAc        >  1:1066243/3‑91 (MQ=255)
                                                                                     atCGAAATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTAACt       >  1:2487798/2‑91 (MQ=255)
                                                                                      tCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTAACtt      <  1:1786886/91‑1 (MQ=255)
                                                                                      tCGAAATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTAACtt      >  1:28636/1‑91 (MQ=255)
                                                                                      tCGAAATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTAACtt      >  2:2172222/1‑91 (MQ=255)
                                                                                          |                                                                                          
CTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTAACTTGTGT  >  USA300TCH1516_ALE/1446955‑1447135

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: