Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,447,040 T→A 6.0% E65D (GAA→GAT cobT ← Aerobic cobaltochelatase subunit CobT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,447,0400TA 6.0% 792.2 / 28.9 233E65D (GAA→GATcobTAerobic cobaltochelatase subunit CobT
Reads supporting (aligned to +/- strand):  ref base T (110/108);  new base A (3/11);  total (114/119)
Fisher's exact test for biased strand distribution p-value = 5.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 240 total aligned reads displayed.

TTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTAACT  >  USA300TCH1516_ALE/1446950‑1447130
                                                                                          |                                                                                          
tttATCTTGTACTACCTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAACTAACATAACATCTGTTTTTAAACCAGCCAtt                                                                                            >  1:3309454/1‑91 (MQ=255)
 ttATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCAttt                                                                                           >  2:1647846/1‑91 (MQ=255)
   aTCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCg                                                                                         >  2:3060685/1‑91 (MQ=255)
    tCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTACACCAGCCATTTCGa                                                                                        >  2:1193443/1‑91 (MQ=255)
     cTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATATCGaa                                                                                       <  2:867323/91‑2 (MQ=255)
     cTTGTACTACTCCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGat                                                                                       <  2:2593164/91‑1 (MQ=255)
      ttGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGata                                                                                      <  2:1834492/91‑1 (MQ=255)
         tactacTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATgg                                                                                   <  1:3540031/91‑1 (MQ=255)
          actactTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGc                                                                                  <  1:111659/91‑1 (MQ=255)
            tactTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCtt                                                                                <  2:1430392/91‑1 (MQ=255)
            tactTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCtt                                                                                <  2:2613250/91‑1 (MQ=255)
            tactTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATATCGAAATGGCtg                                                                                <  1:2453250/91‑2 (MQ=255)
             actTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTg                                                                               >  2:229047/1‑91 (MQ=255)
              ctTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTgg                                                                              <  1:405414/91‑1 (MQ=255)
               ttCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGt                                                                             <  2:204312/91‑1 (MQ=255)
                tCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTc                                                                            >  2:56112/1‑91 (MQ=255)
                tCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATCTCGATATGGCTTGGTc                                                                            >  2:2983091/1‑91 (MQ=255)
                 cATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCg                                                                           >  1:1133775/1‑91 (MQ=255)
                  aTTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTcac                                                                          <  2:3720234/91‑3 (MQ=255)
                  aTTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGa                                                                          <  1:3212613/91‑1 (MQ=255)
                   tttACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTcaca                                                                         <  1:792741/91‑4 (MQ=255)
                   tttACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTcaca                                                                         <  1:3882100/91‑4 (MQ=255)
                    ttACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGAtg                                                                        >  1:281465/1‑91 (MQ=255)
                     tACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGAtgt                                                                       <  2:1414782/91‑1 (MQ=255)
                       cAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGa                                                                     >  2:3694684/1‑91 (MQ=255)
                       cAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGa                                                                     >  2:493081/1‑91 (MQ=255)
                         aTTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGAcc                                                                   >  2:1130185/1‑91 (MQ=255)
                         aTTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGAcc                                                                   <  1:306704/91‑1 (MQ=255)
                           ttGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCaa                                                                 <  2:2097129/91‑1 (MQ=255)
                            tGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCaaa                                                                >  1:417400/1‑91 (MQ=255)
                            tGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCaaa                                                                >  2:3547475/1‑91 (MQ=255)
                            tGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCaaa                                                                <  1:2006996/91‑1 (MQ=255)
                             gTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCaaaa                                                               <  1:1607183/91‑1 (MQ=255)
                             gTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCaaaa                                                               <  1:3666518/91‑1 (MQ=255)
                                catccatcATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCa                                                            <  1:2983279/91‑1 (MQ=255)
                                catccatcATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCa                                                            >  2:2105587/1‑91 (MQ=255)
                                 atccatcATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCaa                                                           <  2:1826458/91‑1 (MQ=255)
                                  tccatcATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCa                                                                  >  1:2261339/1‑83 (MQ=255)
                                    catcatATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTc                                                        <  2:1726615/91‑1 (MQ=255)
                                     atcatATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCg                                                       <  1:1710477/91‑1 (MQ=255)
                                         tatGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAAt                                                   >  1:379102/1‑91 (MQ=255)
                                         tatGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAAt                                                   <  1:542750/91‑1 (MQ=255)
                                          atGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATc                                                  >  2:3782253/1‑91 (MQ=255)
                                           tGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCa                                                 <  1:3846727/91‑1 (MQ=255)
                                           tGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGCAATCa                                                 >  2:3147488/1‑91 (MQ=255)
                                            gATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCaa                                                >  2:670669/1‑91 (MQ=255)
                                            gATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCaa                                                >  2:1081522/1‑91 (MQ=255)
                                             aTAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCaaa                                               >  2:3373424/1‑91 (MQ=255)
                                             aTAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCaaa                                               >  2:1231390/1‑91 (MQ=255)
                                               aaCCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAAc                                             >  2:2065192/1‑91 (MQ=255)
                                                 ccGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACAt                                           <  2:1957872/91‑1 (MQ=255)
                                                 ccGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACAt                                           <  1:1102295/91‑1 (MQ=255)
                                                   gTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACAttt                                         >  1:1226254/1‑91 (MQ=255)
                                                   gTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACAttt                                         <  1:1454780/91‑1 (MQ=255)
                                                     aaGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTg                                       <  1:194824/91‑1 (MQ=255)
                                                      agcagcTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGa                                      >  2:2820252/1‑91 (MQ=255)
                                                          gcTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAg                                  >  1:2059779/1‑91 (MQ=255)
                                                          gcTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAg                                  >  1:725768/1‑91 (MQ=255)
                                                          gcTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAg                                  <  2:2417438/91‑1 (MQ=255)
                                                          gcTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAg                                  >  2:1985115/1‑91 (MQ=255)
                                                           cTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAgg                                 >  1:1411612/1‑91 (MQ=255)
                                                            taacatatcaTCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGa                                <  1:2445567/91‑1 (MQ=255)
                                                            taacataacaTCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGa                                >  1:1640321/1‑91 (MQ=255)
                                                            taacataacaTCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGa                                <  1:1923562/91‑1 (MQ=255)
                                                            taacataacaTCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGa                                <  2:3282024/91‑1 (MQ=255)
                                                                ataacatCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACAtt                                          <  1:1556513/77‑1 (MQ=255)
                                                                ataacatCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATtat                            >  2:5354/1‑91 (MQ=255)
                                                                 taacatCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATtata                           >  1:2116223/1‑91 (MQ=255)
                                                                 taacatCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATtata                           >  2:240648/1‑91 (MQ=255)
                                                                  aacatCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAg                          >  1:2502694/1‑91 (MQ=255)
                                                                   acatCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGt                         >  1:1799145/1‑91 (MQ=255)
                                                                      tCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAg                      <  2:2746333/91‑1 (MQ=255)
                                                                      tCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAg                      <  1:766471/91‑1 (MQ=255)
                                                                       cTGTTTTTAAACCAGCCATTTCGATATTGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAgg                     <  2:1618278/91‑1 (MQ=255)
                                                                       cTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAgg                     <  1:3172571/91‑1 (MQ=255)
                                                                       cTGTTTTTAAACCAGCCATTTCGATATGGCCTGGTCGATGTGACCAAAACCAACTCGTAATCAACACATTTTGAACAGGATTATAGTAAgg                     >  2:387210/1‑91 (MQ=255)
                                                                        tGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAggg                    >  1:1336006/1‑91 (MQ=255)
                                                                         gTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGa                   >  2:3309244/1‑91 (MQ=255)
                                                                         gTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGa                   <  1:3511858/91‑1 (MQ=255)
                                                                          tttttAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGaa                  >  2:582712/1‑91 (MQ=255)
                                                                          tttttAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGaa                  <  2:3507933/91‑1 (MQ=255)
                                                                           ttttAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAAt                 <  1:837500/91‑1 (MQ=255)
                                                                            tttAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAAtt                <  2:3566007/91‑1 (MQ=255)
                                                                                aaCCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTg            <  2:2822133/91‑1 (MQ=255)
                                                                                aaCCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTg            <  2:46060/91‑1 (MQ=255)
                                                                                 ccaaGCCATATCGAAATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTgg           >  1:935706/4‑91 (MQ=255)
                                                                                 aCCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTgg           >  1:2747957/1‑91 (MQ=255)
                                                                                  ccAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGa          >  2:3168851/1‑91 (MQ=255)
                                                                                   cAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAAGCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTCTGGAt         >  1:3612747/1‑91 (MQ=255)
                                                                                   cAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGAt         >  1:886723/1‑91 (MQ=255)
                                                                                    aGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGAtt        >  1:2871244/1‑91 (MQ=255)
                                                                                     gCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGAttt       >  2:1460341/1‑91 (MQ=255)
                                                                                     gCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGAttt       >  1:3390497/1‑91 (MQ=255)
                                                                                      ccATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGAtttt      <  2:3788541/91‑1 (MQ=255)
                                                                                       cATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTa     >  1:2745425/1‑91 (MQ=255)
                                                                                       cATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTa     >  1:1235501/1‑91 (MQ=255)
                                                                                       cATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTa     >  2:2206841/1‑91 (MQ=255)
                                                                                        aTTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTaa    <  2:1513740/91‑1 (MQ=255)
                                                                                         tttCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTAAc   <  1:2215066/91‑1 (MQ=255)
                                                                                          ttCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGCATTTTAACt  >  1:3720047/1‑91 (MQ=255)
                                                                                          |                                                                                          
TTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTAACT  >  USA300TCH1516_ALE/1446950‑1447130

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: