Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,517,782 T→A 12.1% L950L (CTA→CTT narG ← Respiratory nitrate reductase 1 alpha chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,517,7820TA12.1% 552.2 / 56.4 183L950L (CTA→CTTnarGRespiratory nitrate reductase 1 alpha chain
Reads supporting (aligned to +/- strand):  ref base T (71/89);  new base A (7/15);  total (78/105)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 190 total aligned reads displayed.

AGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAAC  >  USA300TCH1516_ALE/2517692‑2517872
                                                                                          |                                                                                          
  aaaTATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAg                                                                                          <  1:3396829/91‑1 (MQ=255)
  aaaTATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAg                                                                                          <  1:1847724/91‑1 (MQ=255)
  aaaTATCTTTTAACGACATTCCAGTTTGTTCTTCGAGATCTTCATATGCTTTTTGTGATACTTTACCATTCGTAGCAGATGAAATACTTAg                                                                                          >  2:1579412/1‑91 (MQ=255)
   aaTATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGCAATACTAAGt                                                                                         <  1:752097/91‑1 (MQ=255)
   aaTATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGt                                                                                         <  2:3394889/91‑1 (MQ=255)
   aaTATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGt                                                                                         >  2:726169/1‑91 (MQ=255)
    atatCTTTTAACGGCTTTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTCCCATTCGTAGCAGATGAAATACTTAGTa                                                                                        <  2:827964/91‑1 (MQ=255)
     tatCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTAt                                                                                       >  2:2064904/1‑91 (MQ=255)
     tatCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTAt                                                                                       <  2:1260840/91‑1 (MQ=255)
      atCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTAtt                                                                                      >  1:2328096/1‑91 (MQ=255)
       tCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTAttt                                                                                     <  1:2892364/91‑2 (MQ=38)
       tCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGCAATACTAAGTAttt                                                                                     <  2:3046473/91‑2 (MQ=255)
        cTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTAtttc                                                                                    <  1:2342482/91‑3 (MQ=38)
        cTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGc                                                                                    <  2:1006879/91‑1 (MQ=255)
          tttAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCAt                                                                                  >  2:1758167/1‑91 (MQ=255)
          tttAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCAt                                                                                  <  2:2714821/91‑1 (MQ=255)
           ttAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTATTTCATc                                                                                 <  1:3257492/91‑1 (MQ=38)
           ttAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTATTTCATc                                                                                 <  1:1165544/91‑1 (MQ=38)
           ttAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTATTTCATc                                                                                 <  1:907110/91‑1 (MQ=37)
             aaCGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTTTTTTATTTTTCCATTCGTAGCAGATGAAATACTTAGTATTGCATCAg                                                                               >  2:2369488/1‑91 (MQ=255)
              aCGGCATTCCAGTTTGTTCTTCGAGATCTTCGTATGATTTTTGTGATAATTTAGCATGCGTAGCAGATGAAATACTTAGTATTGCATCAGc                                                                              <  1:3157735/91‑1 (MQ=255)
               cGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTATTTCATCTGCt                                                                             <  2:2603314/91‑1 (MQ=25)
               cGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCt                                                                             <  2:2418777/91‑1 (MQ=255)
                ggCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTa                                                                            <  1:1902760/91‑1 (MQ=255)
                ggCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTa                                                                            <  2:2677416/91‑1 (MQ=255)
                 gCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTAc                                                                           <  1:1783855/91‑1 (MQ=255)
                  cATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACa                                                                          >  1:2157444/1‑91 (MQ=255)
                  cATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACa                                                                          >  1:1672946/1‑91 (MQ=255)
                    ttCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAtt                                                                        <  2:470859/91‑1 (MQ=255)
                    ttCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAtt                                                                        <  1:104777/91‑1 (MQ=255)
                    ttCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAtt                                                                        <  1:3582912/91‑1 (MQ=255)
                    ttCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAtt                                                                        >  2:3487059/1‑91 (MQ=255)
                      ccAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTAc                                                                      <  2:2542254/91‑1 (MQ=255)
                       cAGTTGGTTCTTCGAGATCTTCATATGATTGTTGTGATGATTTAGCATTGGTAGCAGATGAGATACTTAGTATTGCATCAGCTACATTACg                                                                     <  2:591350/91‑1 (MQ=255)
                        aGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGt                                                                    >  1:443100/1‑91 (MQ=255)
                          tttGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGc                                                                  >  2:2317227/1‑91 (MQ=255)
                            tGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTg                                                                >  1:1272400/1‑91 (MQ=255)
                            tGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTg                                                                >  1:1745314/1‑91 (MQ=255)
                              ttcttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTa                                                              <  1:1191144/91‑1 (MQ=255)
                               tcttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTAt                                                             <  1:3455753/91‑1 (MQ=255)
                               tcttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTAt                                                             <  2:3221637/91‑1 (MQ=255)
                                cttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATc                                                            <  1:3845333/91‑1 (MQ=255)
                                  tcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCaa                                                          <  1:1364718/91‑1 (MQ=255)
                                    gagaTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATa                                                        <  1:3663749/91‑1 (MQ=255)
                                       aTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGt                                                     <  1:1651330/91‑1 (MQ=255)
                                       aTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAACTACATTACGTGCTGTATCAATACGt                                                     <  2:2200769/91‑1 (MQ=255)
                                        tCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTg                                                    >  1:2163110/1‑91 (MQ=255)
                                         cTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTgg                                                   >  2:143906/1‑91 (MQ=255)
                                         cTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTgg                                                   <  2:644491/91‑1 (MQ=255)
                                         cTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTgg                                                   >  1:3542505/1‑91 (MQ=255)
                                            cATATGATTTTTGTGATAATTTACCATTCGTATCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACg                                                <  1:1234099/91‑1 (MQ=255)
                                            cATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACg                                                <  2:3231986/91‑1 (MQ=255)
                                              tatGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCGGCTACATTACGTGCTGTATCAATACGTGGACGAt                                              <  2:3619143/91‑1 (MQ=255)
                                              tatGATTATTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGAt                                              <  1:1953284/91‑1 (MQ=255)
                                               atGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGAtt                                             >  2:2088317/1‑91 (MQ=255)
                                                 gATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCg                                           <  2:3641994/91‑1 (MQ=255)
                                                    ttttGTGATAATTTACCATTCGTAGCAGATGAAATACTGAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGctc                                        <  2:2957818/91‑1 (MQ=255)
                                                     tttGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTTCATTACGTGCTGTATCAATACGTGGACGATTCGctct                                       <  1:3340429/91‑1 (MQ=255)
                                                     tttGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGctct                                       >  2:3537280/1‑91 (MQ=255)
                                                     tttGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGctct                                       <  1:2953973/91‑1 (MQ=255)
                                                     tttGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGctct                                       <  1:666547/91‑1 (MQ=255)
                                                       tgtgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTca                                     <  2:2926893/91‑1 (MQ=255)
                                                       tgtgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTca                                     <  1:1235559/91‑1 (MQ=255)
                                                        gtgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTcac                                    >  1:634149/1‑91 (MQ=255)
                                                        gtgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTcac                                    <  1:3535159/91‑1 (MQ=255)
                                                          gATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAg                                  <  2:3894786/91‑1 (MQ=255)
                                                          gATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAg                                  >  1:2617120/1‑91 (MQ=255)
                                                           ttAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCTCAGa                                 >  2:3127895/2‑91 (MQ=255)
                                                           aTAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGa                                 <  1:2780704/91‑1 (MQ=255)
                                                           aTAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGa                                 >  2:3127790/1‑91 (MQ=255)
                                                            tAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGaa                                <  1:1037870/91‑1 (MQ=255)
                                                             aaTTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAat                               <  2:3844989/91‑1 (MQ=255)
                                                               tttACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAatca                             <  1:2003992/91‑1 (MQ=255)
                                                               tttACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAatca                             >  1:2052680/1‑91 (MQ=255)
                                                                ttACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAatcat                            >  2:1289828/1‑91 (MQ=255)
                                                                ttACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAatcat                            >  1:218425/1‑91 (MQ=255)
                                                                    cATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGTATCattatt                        >  2:1640115/1‑91 (MQ=255)
                                                                    cATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGCATCATCAtt                        >  2:2481983/1‑91 (MQ=255)
                                                                    cATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCAtt                        <  2:1147632/91‑1 (MQ=255)
                                                                      ttCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTg                      >  2:1946329/1‑91 (MQ=255)
                                                                       tCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGt                     >  1:243572/1‑91 (MQ=255)
                                                                       tCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGt                     >  1:242829/1‑91 (MQ=255)
                                                                         gTAGCTGATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTAt                   >  1:1322373/1‑91 (MQ=25)
                                                                            gCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTTCTGTATCAATACGTTGGCGATTCGCTCTCACAGAATCATCATTTGTATccc                >  2:3280950/1‑89 (MQ=255)
                                                                              aGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCAGTc              >  2:1202678/1‑91 (MQ=255)
                                                                              aGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTc              <  2:901614/91‑1 (MQ=255)
                                                                               gATGCAATACTAAGTATTTCATCTGCCACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTcc             >  2:2238965/1‑91 (MQ=38)
                                                                                aTGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCa            <  1:697701/91‑1 (MQ=255)
                                                                                  gcaaTACTAAGTATTTCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACg          >  1:3754757/3‑91 (MQ=255)
                                                                                  gcaaTACTAAGTATTTCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACg          >  1:2534709/3‑91 (MQ=255)
                                                                                       acgtAGTATTGCATCAGCTACATTACGGGCTGTAGCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCt     <  2:3762962/88‑1 (MQ=255)
                                                                                       aCTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCt     <  2:2577761/91‑1 (MQ=255)
                                                                                       aCTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCt     <  2:3763016/91‑1 (MQ=255)
                                                                                       aCTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCt     >  1:299123/1‑91 (MQ=255)
                                                                                       aCTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCt     >  1:298663/1‑91 (MQ=255)
                                                                                       aCTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCt     >  1:2041804/1‑91 (MQ=255)
                                                                                       aCTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCACACAGAATCATCATTTGTATCACACCACGTACCt     <  2:1107046/91‑1 (MQ=255)
                                                                                        cTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTa    <  1:139368/91‑1 (MQ=255)
                                                                                         ttAGTATTGGATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTaa   <  1:3787223/91‑1 (MQ=255)
                                                                                          tAGCATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAAc  <  1:3062861/91‑1 (MQ=255)
                                                                                          |                                                                                          
AGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAAC  >  USA300TCH1516_ALE/2517692‑2517872

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: