Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,447,045 T→A 6.3% 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.3% 844.2 / 33.3 252I64F (ATC→TTC) cobTAerobic cobaltochelatase subunit CobT
Reads supporting (aligned to +/- strand):  ref base T (121/114);  new base A (13/3);  total (135/117)
Fisher's exact test for biased strand distribution p-value = 3.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

Only 100 of 263 total aligned reads displayed.

CTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTAACTTGTG  >  USA300TCH1516_ALE/1446955‑1447134
                                                                                          |                                                                                         
cTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGat                                                                                           <  1:2579463/91‑1 (MQ=255)
 ttGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGata                                                                                          <  1:2840454/91‑1 (MQ=255)
  tGTACTACTTCATTTACAATTTGTACATCCATCATCTGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATATCGaaat                                                                                         <  2:2753798/91‑1 (MQ=255)
   gTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATg                                                                                        >  2:325082/1‑91 (MQ=255)
   gTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATATCGAAATg                                                                                        <  2:445184/91‑1 (MQ=255)
    tactacTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATgg                                                                                       <  2:3262250/91‑1 (MQ=255)
     actactTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGc                                                                                      <  1:3707086/91‑1 (MQ=255)
     actactTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGc                                                                                      >  1:872185/1‑91 (MQ=255)
     actactTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGc                                                                                      <  2:2000554/91‑1 (MQ=255)
      ctactTCATTTACGATTTGTACATCCATCATATGATAACCGTAGGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCt                                                                                     <  1:3949280/91‑1 (MQ=255)
      ctactGCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCt                                                                                     <  2:1822956/91‑1 (MQ=255)
       tactTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCtt                                                                                    >  2:144714/1‑91 (MQ=255)
       tactTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCtt                                                                                    <  1:202000/91‑1 (MQ=255)
       tactTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCtt                                                                                    >  2:2836639/1‑91 (MQ=255)
        actTCATTTACAATTTGTACATCCATCATATGATAACCGTCAGCAGCTAACATAACATCTGTTCTTGGACCAGCCATTTCCATATGGCTTg                                                                                   >  2:3667987/1‑91 (MQ=255)
          ttCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGt                                                                                 <  1:1009413/91‑1 (MQ=255)
           tCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTc                                                                                <  1:737888/91‑1 (MQ=255)
           tCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTc                                                                                >  1:1362690/1‑91 (MQ=255)
            cATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCg                                                                               >  2:2451902/1‑91 (MQ=255)
              tttACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTcaca                                                                             <  2:1793696/91‑4 (MQ=255)
               ttACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGAtg                                                                            >  2:3751559/1‑91 (MQ=255)
               ttACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGAtg                                                                            <  1:2756624/91‑1 (MQ=255)
                tACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGAtgt                                                                           <  1:2930687/91‑1 (MQ=255)
                tACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGAtgt                                                                           <  2:162991/91‑1 (MQ=255)
                  cAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGa                                                                         >  2:2042045/1‑91 (MQ=255)
                   aaTTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGAc                                                                        <  2:3531048/91‑1 (MQ=255)
                     tttGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCa                                                                      <  1:375280/91‑1 (MQ=255)
                      ttGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCaa                                                                     <  1:2808193/91‑1 (MQ=255)
                       tGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCaaa                                                                    <  1:84742/91‑1 (MQ=255)
                        gTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCaaaa                                                                   <  2:1474137/91‑1 (MQ=255)
                         tACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAAc                                                                  >  1:2213845/1‑91 (MQ=255)
                          aCATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAAcc                                                                 >  1:2457159/1‑91 (MQ=255)
                          aCATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAAcc                                                                 <  1:1716539/91‑1 (MQ=255)
                           catccatcATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCa                                                                <  1:1918498/91‑1 (MQ=255)
                             tccatcATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAAc                                                              >  1:687668/1‑91 (MQ=255)
                              ccatcATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACt                                                             >  1:1311106/1‑91 (MQ=255)
                              ccatcACATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACt                                                             >  1:3069478/1‑91 (MQ=255)
                                atcatATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCg                                                           >  2:252457/1‑91 (MQ=255)
                                atcatATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCg                                                           <  1:3398736/91‑1 (MQ=255)
                                 tcatATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATATGACCAAAACCAACTCGt                                                          >  2:2677311/1‑91 (MQ=255)
                                  catATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTa                                                         >  2:3561486/1‑91 (MQ=255)
                                     atGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATc                                                      >  1:866839/1‑91 (MQ=255)
                                     atGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATc                                                      <  2:400253/91‑1 (MQ=255)
                                     atGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATc                                                      >  2:2967262/1‑91 (MQ=255)
                                      tGATAACCGTAAGCAGCTAACATAACACCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCa                                                     <  2:2819573/91‑1 (MQ=255)
                                        aTAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCaaa                                                   >  1:155133/1‑91 (MQ=255)
                                         tAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCaaaa                                                  >  2:963725/1‑91 (MQ=255)
                                          aaCCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAAc                                                 <  2:2562595/91‑1 (MQ=255)
                                            ccGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACAt                                               <  1:3197165/91‑1 (MQ=255)
                                             cGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACAtt                                              >  2:745077/1‑91 (MQ=255)
                                                 agcagcTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGa                                          >  1:1325586/1‑91 (MQ=255)
                                                 agcagcTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGa                                          >  1:1982675/1‑91 (MQ=255)
                                                 agcagcTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGa                                          <  1:2166948/91‑1 (MQ=255)
                                                 agcagcTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGa                                          >  1:975886/1‑91 (MQ=255)
                                                  gcagcTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCACAACCAACTCGTAATCAACACATTTTgca                                         >  2:2620868/1‑89 (MQ=255)
                                                  gcagcTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGaa                                         <  1:1792592/91‑1 (MQ=255)
                                                  gcagcTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGaa                                         >  2:2253296/1‑91 (MQ=255)
                                                    agcTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACa                                       <  1:3717352/91‑1 (MQ=255)
                                                     gcTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAg                                      >  1:573482/1‑91 (MQ=255)
                                                     gcTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAg                                      <  1:276596/91‑1 (MQ=255)
                                                       taacataacaTCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGa                                    <  1:485212/91‑1 (MQ=255)
                                                       taacataacaTCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGa                                    <  2:1859994/91‑1 (MQ=255)
                                                       taacataacaTCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGa                                    >  1:3018842/1‑91 (MQ=255)
                                                        aacataacatCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGAt                                   >  1:2677834/1‑91 (MQ=255)
                                                         acataacatCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGAtt                                  >  2:800832/1‑91 (MQ=255)
                                                         acataacatCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGAtt                                  >  1:1245488/1‑91 (MQ=255)
                                                          cataacatCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATta                                 <  1:2720184/91‑1 (MQ=255)
                                                             aacatCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAg                              >  2:3838444/1‑91 (MQ=255)
                                                               catCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTa                            <  1:2309958/91‑1 (MQ=255)
                                                               catCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTa                            >  1:1574442/1‑91 (MQ=255)
                                                                 tCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAg                          <  1:761524/91‑1 (MQ=255)
                                                                 tCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGGCGATGTGACCAAAACCAACTCGTAATCAAAACAGTTTGAACAGGATTATAGTAAg                          <  1:1838918/91‑1 (MQ=255)
                                                                 tCTATTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCACCTCGTAATCAACACATTTTGAACAGGACTATAGTAAg                          >  2:3505231/1‑91 (MQ=255)
                                                                  cTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAgg                         >  2:3630659/1‑91 (MQ=255)
                                                                   tGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCCAAACCAACTCGTAATCAAAACCTTTTGAACAGGATTATAGTAAggg                        >  2:441454/1‑91 (MQ=255)
                                                                   tGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAggg                        <  2:1298368/91‑1 (MQ=255)
                                                                      ttttAAACCAGCCATTTCGATATGGCTTGATCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGAAAGGGAAt                     >  2:3438169/1‑91 (MQ=255)
                                                                       tttAAACCAGCCATTTCGATATGTATTGGTCGATGTGACCAAAACCAACTTGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAAtt                    <  2:3679528/91‑1 (MQ=255)
                                                                        ttAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAAttt                   <  1:148328/91‑1 (MQ=255)
                                                                         tAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAAtttt                  <  2:2299721/91‑1 (MQ=255)
                                                                           aaCCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTg                <  1:1686204/91‑1 (MQ=255)
                                                                            aCCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACCTTTTGAACAGGATTATAGTAAGGGAATTTTTgg               >  2:300799/1‑91 (MQ=255)
                                                                               aGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGAtt            <  1:2400591/91‑1 (MQ=255)
                                                                               aGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGAtt            <  1:2133609/91‑1 (MQ=255)
                                                                               aGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGAtt            <  1:1146556/91‑1 (MQ=255)
                                                                                  catatCGAAATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTa         >  2:3117453/5‑91 (MQ=255)
                                                                                   aTTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTaa        <  2:2873574/91‑1 (MQ=255)
                                                                                   aTTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTaa        >  1:1132750/1‑91 (MQ=255)
                                                                                    tttCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTAAc       >  2:3411115/1‑91 (MQ=255)
                                                                                    tatCGAAATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGCTTATAGTAAGGGAATTTTTGGATTTTAAc       >  2:1056237/3‑91 (MQ=255)
                                                                                    tatCGAAATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGCTTTTAAc       >  1:344669/3‑91 (MQ=255)
                                                                                     ttCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTAACt      <  1:2197751/91‑1 (MQ=255)
                                                                                      tCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTAACtt     <  1:1537128/91‑1 (MQ=255)
                                                                                      tCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTAACtt     <  1:1093454/91‑1 (MQ=255)
                                                                                      tCGAAATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTAACtt     >  1:3307695/1‑91 (MQ=255)
                                                                                        gATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTAACTtgt   >  1:2474598/1‑91 (MQ=255)
                                                                                        gATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACACTTTGAACAGGATTATAGTAAGGGACTTTTTGGATTTTAACTtgt   >  1:3143713/1‑91 (MQ=255)
                                                                                         atatGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTAACTtgtg  <  1:3793776/91‑1 (MQ=255)
                                                                                         atatGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTAACTtgtg  >  2:639393/1‑91 (MQ=255)
                                                                                         atatGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTAACTtgtg  <  1:3793821/91‑1 (MQ=255)
                                                                                          |                                                                                         
CTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAGGATTATAGTAAGGGAATTTTTGGATTTTAACTTGTG  >  USA300TCH1516_ALE/1446955‑1447134

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