Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,517,782 T→A 11.2% 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,7820TA11.2% 500.2 / 44.7 160L950L (CTA→CTTnarGRespiratory nitrate reductase 1 alpha chain
Reads supporting (aligned to +/- strand):  ref base T (81/61);  new base A (9/9);  total (90/70)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 164 total aligned reads displayed.

AGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAAC  >  USA300TCH1516_ALE/2517692‑2517872
                                                                                          |                                                                                          
aGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACtt                                                                                            <  1:3259962/91‑1 (MQ=255)
 gAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTTTGATAATTTACCATTCGTAGCAGATGAAATACTTa                                                                                           >  2:1073550/1‑91 (MQ=255)
 gAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTa                                                                                           <  1:787550/91‑1 (MQ=255)
 gAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTa                                                                                           >  1:2238300/1‑91 (MQ=255)
 gAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTa                                                                                           >  1:1613800/1‑91 (MQ=255)
   aaTATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGCAATACTAAGt                                                                                         <  1:3032636/91‑1 (MQ=255)
      atCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTTTTTTTTTTTTTCATTCGTAGCAGATGAAATACTTAGTAtt                                                                                      >  2:1560957/1‑91 (MQ=255)
      atCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTAtt                                                                                      >  2:3150191/1‑91 (MQ=255)
       tCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTATCATTCGTAGCAGATGAAATACTTAGTATTg                                                                                     >  1:1775618/1‑91 (MQ=255)
       tCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTg                                                                                     >  2:2668846/1‑91 (MQ=255)
        cTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTAtttc                                                                                    <  2:1930537/91‑3 (MQ=38)
        cTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTAtttc                                                                                    <  2:2968448/91‑3 (MQ=38)
        cTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGc                                                                                    <  1:402928/91‑1 (MQ=255)
        cTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGc                                                                                    <  1:403337/91‑1 (MQ=255)
        cTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGc                                                                                    <  1:363764/91‑1 (MQ=255)
        cTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGc                                                                                    <  2:831438/91‑1 (MQ=255)
         ttttAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCa                                                                                   >  1:1323241/1‑91 (MQ=255)
          tttAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCAt                                                                                  <  1:727170/91‑1 (MQ=255)
          tttAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCAt                                                                                  >  1:3069822/1‑91 (MQ=255)
            tAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTATTTCATct                                                                                <  1:1737437/91‑2 (MQ=37)
             aaCGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTATTTCATctg                                                                               <  1:403613/91‑3 (MQ=37)
             aaCGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAg                                                                               <  1:3867079/91‑1 (MQ=255)
             aaCGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAg                                                                               >  1:291443/1‑91 (MQ=255)
              aCGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTATTTCATctgc                                                                              <  1:715999/91‑4 (MQ=25)
                tgCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTa                                                                            >  2:1031553/2‑91 (MQ=255)
                 gCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTATTTCATCTGCTAc                                                                           <  1:906308/91‑1 (MQ=25)
                 gCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTATTTCATCTGCTAc                                                                           <  1:906258/91‑1 (MQ=25)
                  cATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTATCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACa                                                                          >  2:3414611/1‑91 (MQ=255)
                  cATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACa                                                                          >  2:2179606/1‑91 (MQ=255)
                  cATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACa                                                                          >  2:1854316/1‑91 (MQ=255)
                  cATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACa                                                                          >  2:3361103/1‑91 (MQ=255)
                  cATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACa                                                                          >  2:2494173/1‑91 (MQ=255)
                   aTTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAt                                                                         >  1:1700113/1‑91 (MQ=255)
                    ttCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAtt                                                                        <  2:2887923/91‑1 (MQ=255)
                    ttCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAtt                                                                        >  1:3423790/1‑91 (MQ=255)
                    ttCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAtt                                                                        >  1:3423805/1‑91 (MQ=255)
                      ccAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTAc                                                                      <  2:1905257/91‑1 (MQ=255)
                      ccAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTAc                                                                      >  1:925244/1‑91 (MQ=255)
                       cAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACg                                                                     >  2:645123/1‑91 (MQ=255)
                        aGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGTTATTTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGt                                                                    >  2:1481574/1‑91 (MQ=255)
                          tttGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGc                                                                  <  1:2349036/91‑1 (MQ=255)
                          tttGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGc                                                                  <  2:2163068/91‑1 (MQ=255)
                          tttGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGc                                                                  <  2:936215/91‑1 (MQ=255)
                           ttGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCt                                                                 >  1:1387906/1‑91 (MQ=255)
                            tGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTg                                                                >  2:1566952/1‑91 (MQ=255)
                             gTTCTTCGAGATCTTCATATGATTTTTGTGTTATTTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGt                                                               >  2:130497/1‑91 (MQ=255)
                             gTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGc                                                                              <  1:2680916/76‑1 (MQ=255)
                             gTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGt                                                               <  2:2310775/91‑1 (MQ=255)
                                cttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATc                                                            <  2:389699/91‑1 (MQ=255)
                                 ttcGAGATCTTCATATGATTTTTGTTATTTTTTATCATTTGTAGCAGATTAAATACTTAGTATTTCTTCATCTACATTACGTGCTGTATCa                                                           >  2:3620431/1‑91 (MQ=255)
                                 ttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCa                                                           <  1:102840/91‑1 (MQ=255)
                                   cGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAAt                                                         >  2:2179167/1‑91 (MQ=255)
                                      gaTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACg                                                      >  1:2359765/1‑91 (MQ=255)
                                      gaTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACg                                                      >  2:2976636/1‑91 (MQ=255)
                                      gaTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACg                                                      <  2:2008346/91‑1 (MQ=255)
                                       aTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGt                                                     >  2:76176/1‑91 (MQ=255)
                                         cTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTAc                                                                           <  1:2297688/67‑1 (MQ=255)
                                            cATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACg                                                <  1:745525/91‑1 (MQ=255)
                                            cATATAATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACg                                                <  2:173125/91‑1 (MQ=255)
                                             atatGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGa                                               <  1:1661531/91‑1 (MQ=255)
                                             atatGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGa                                               >  1:187602/1‑91 (MQ=255)
                                               atGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGAtt                                             >  2:2683949/1‑91 (MQ=255)
                                                tGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTc                                            <  1:190619/91‑1 (MQ=255)
                                                   tttttGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGct                                         <  1:2686428/91‑1 (MQ=255)
                                                        gtgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTcac                                    <  1:982790/91‑1 (MQ=255)
                                                        gtgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTcac                                    >  1:779216/1‑91 (MQ=255)
                                                        gtgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTcac                                    <  1:982443/91‑1 (MQ=255)
                                                         tgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTcaca                                   >  2:2787411/1‑91 (MQ=255)
                                                             aaTTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAat                               <  1:898643/91‑1 (MQ=255)
                                                              aTTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAatc                              >  2:91245/1‑91 (MQ=255)
                                                               tttACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAatca                             >  2:3608573/1‑91 (MQ=255)
                                                                ttACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAatcat                            <  2:2629801/91‑1 (MQ=255)
                                                                 tACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAatcatc                           <  2:878240/91‑1 (MQ=255)
                                                                   ccATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATcatcat                         >  2:3286130/1‑91 (MQ=255)
                                                                   ccATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATcatcat                         >  2:780935/1‑91 (MQ=255)
                                                                      ttCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGctct                                       <  1:1016991/74‑1 (MQ=255)
                                                                      ttCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTg                      <  1:2126599/91‑1 (MQ=255)
                                                                       tCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGt                     >  2:1095903/1‑91 (MQ=255)
                                                                        tgTAGCTGATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTa                    >  2:2239666/2‑91 (MQ=25)
                                                                        cGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTa                    <  1:503071/91‑1 (MQ=255)
                                                                        cGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTa                    <  1:502073/91‑1 (MQ=255)
                                                                         gTAGCTGATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTAt                   >  2:1144664/1‑91 (MQ=25)
                                                                          tAGCTGATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATc                  >  2:2540245/1‑91 (MQ=25)
                                                                           agctgATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCa                 >  2:3402483/5‑91 (MQ=25)
                                                                           agctgATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCa                 >  2:3554662/5‑91 (MQ=25)
                                                                           agctgATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCa                 >  1:2749930/5‑91 (MQ=25)
                                                                           aGCTGATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCa                 >  2:312362/1‑91 (MQ=25)
                                                                           aGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCa                 >  2:707178/1‑91 (MQ=255)
                                                                            gCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCAc                >  1:141054/1‑91 (MQ=255)
                                                                              aGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTc              <  1:2502179/91‑1 (MQ=255)
                                                                               gATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTcc             <  1:1488426/91‑1 (MQ=255)
                                                                                aTGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCa            <  1:2827484/91‑1 (MQ=255)
                                                                                  gAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACg          >  2:3743122/1‑91 (MQ=255)
                                                                                  gAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACg          >  1:1384204/1‑91 (MQ=255)
                                                                                  gAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACg          >  1:3855087/1‑91 (MQ=255)
                                                                                    aaTACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGCATCATCATTTGTATCACTCCACGTa        >  2:3522007/1‑91 (MQ=255)
                                                                                      tACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTAcc      >  1:1165685/1‑91 (MQ=255)
                                                                                      tACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTAcc      >  1:1155694/1‑91 (MQ=255)
                                                                                        cTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTa    >  1:635816/1‑91 (MQ=255)
                                                                                         ttAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTaa   <  2:3017289/91‑1 (MQ=255)
                                                                                          |                                                                                          
AGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAAC  >  USA300TCH1516_ALE/2517692‑2517872

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: