Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 724,452 G→A 5.9% E197K (GAA→AAA)  feuC_1 → Iron‑uptake system permease protein FeuC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE724,4520GA 5.9% 653.8 / 17.2 196E197K (GAA→AAA) feuC_1Iron‑uptake system permease protein FeuC
Reads supporting (aligned to +/- strand):  ref base G (101/83);  new base A (2/10);  total (103/93)
Fisher's exact test for biased strand distribution p-value = 1.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 202 total aligned reads displayed.

TATGGTGGCTCAATTACAATTATGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTAAATCGAATACACTAAATATTATTCAT  >  USA300TCH1516_ALE/724362‑724541
                                                                                          |                                                                                         
tatGGTGGCTCAATTACAATTATGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACg                                                                                           >  1:2366527/1‑91 (MQ=255)
tatGGTGGCTCAATTACAATTATGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACg                                                                                           >  1:3555574/1‑91 (MQ=255)
  tggtggCTCAATTACAATTATGTCAAAATTTGATGATAAGTAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGaa                                                                                         >  2:3190727/1‑91 (MQ=255)
  tggtggCTCAATTACAATTATGTCAAAATTTGATGATAAGAAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGaa                                                                                         >  2:3189012/1‑91 (MQ=255)
   ggtggCTCAATTACAATTATGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACAAAt                                                                                        <  2:2356726/91‑1 (MQ=255)
   ggtggCTCAATTACAATTATGTCAAAATTTGATGATAAGCAATCTGATGTCGTCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAAt                                                                                        <  2:1178311/91‑1 (MQ=255)
    gtggCTCACTTACAATTATGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTCATATTTGGGGTGACGAATg                                                                                       <  1:2104409/91‑1 (MQ=255)
    gtggCTCAATTACAATTATGTCAAAATTTGATGATAAGCAATCTGCTTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATg                                                                                       >  1:1680817/1‑91 (MQ=255)
     tggCTCAATTACAATTATGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATgg                                                                                      <  1:3773676/91‑1 (MQ=255)
      ggCTCAATTACAATTATGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGc                                                                                     >  1:2381823/1‑91 (MQ=255)
       gCTCAATTACAATTATGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGcc                                                                                    >  1:2678425/1‑91 (MQ=255)
        cTCAATTACAATTATGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCa                                                                                   <  2:3559717/91‑1 (MQ=255)
        cTCAATTACAATTATGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACAAATGGCCa                                                                                   <  1:2751877/91‑1 (MQ=255)
        cTCAATTACAATTATGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACAAATGGCCa                                                                                   <  2:2780625/91‑1 (MQ=255)
         tCAATTACAATTATGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCAt                                                                                  >  1:934172/1‑91 (MQ=255)
         tCAATTACAATTATGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCAt                                                                                  >  1:2218709/1‑91 (MQ=255)
         tCAATTACAATTATGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACAAATGGCCAt                                                                                  >  1:3711586/1‑91 (MQ=255)
          caattacaattaTGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCAtt                                                                                 >  1:2196272/1‑91 (MQ=255)
            attacaattaTGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTg                                                                               <  1:1131609/91‑1 (MQ=255)
            attacaattaTGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACAAATGGCCAttcg                                                                               <  1:1696018/91‑3 (MQ=255)
             ttacaattaTGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGt                                                                              >  2:2174202/1‑91 (MQ=255)
             ttacaattaTGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGt                                                                              >  1:261293/1‑91 (MQ=255)
             ttacaattaTGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCCTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGt                                                                              <  2:2951285/91‑1 (MQ=255)
               acaattaTGTCTAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCa                                                                            >  2:2666627/1‑91 (MQ=255)
               acaattaTGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCa                                                                            <  2:3862772/91‑1 (MQ=255)
                caattaTGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCAt                                                                           >  2:1432128/1‑91 (MQ=255)
                caattaTGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCAt                                                                           >  2:1560185/1‑91 (MQ=255)
                 aattaTGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCAtt                                                                          >  1:1994737/1‑91 (MQ=255)
                  attaTGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGCATGGCCATTTGTCATTg                                                                         >  1:3327957/1‑91 (MQ=255)
                  attaTGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTg                                                                         <  1:2556493/91‑1 (MQ=255)
                   ttaTGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGc                                                                        <  1:389865/91‑1 (MQ=255)
                    taTGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCa                                                                       >  1:2271748/1‑91 (MQ=255)
                     aTGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCAt                                                                      <  2:142912/91‑1 (MQ=255)
                      tGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCAtt                                                                     <  2:1072844/91‑1 (MQ=255)
                      tGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCAtt                                                                     <  1:1461018/91‑1 (MQ=255)
                      tGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCAtt                                                                     <  1:47652/91‑1 (MQ=255)
                       gTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCAttt                                                                    <  1:674090/91‑1 (MQ=255)
                         cAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCAttttt                                                                  >  2:3334314/1‑91 (MQ=255)
                         cAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCAttttt                                                                  >  2:1458638/1‑91 (MQ=255)
                          aaaaTTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTa                                                                 <  1:976201/91‑1 (MQ=255)
                           aaaTTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTAc                                                                <  2:2046573/91‑1 (MQ=255)
                             aTTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCg                                                              <  1:2263247/91‑1 (MQ=255)
                                tgatgaTAAGCAATCCGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTgg                                                           >  1:3306879/1‑91 (MQ=255)
                                 gatgatAAGCAATCTGATTTCATCGCTGCTTGGTTTGTAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTggg                                                          >  1:931480/1‑91 (MQ=255)
                                 gatgatAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTggg                                                          >  1:1550445/1‑91 (MQ=255)
                                 gatgatAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTggg                                                          >  1:1260188/1‑91 (MQ=255)
                                 gatgatAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTggg                                                          >  1:1486589/1‑91 (MQ=255)
                                  atgatAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGgt                                                         >  1:2570892/1‑91 (MQ=255)
                                  atgatAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGgt                                                         >  1:3203519/1‑91 (MQ=255)
                                  atgatAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGgt                                                         >  1:643877/1‑91 (MQ=255)
                                   tgatAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGgtg                                                        >  1:3723719/1‑91 (MQ=255)
                                   tgatAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGgtg                                                        >  1:3564176/1‑91 (MQ=255)
                                   tgatAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGgtg                                                        >  2:980064/1‑91 (MQ=255)
                                      tAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTg                                                     <  1:3551682/91‑1 (MQ=255)
                                      tAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTg                                                     <  2:1774873/91‑1 (MQ=255)
                                      tAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTg                                                     <  2:1914982/91‑1 (MQ=255)
                                       aaGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGa                                                    <  2:1227762/91‑1 (MQ=255)
                                         gCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGatt                                                  <  2:574160/91‑1 (MQ=255)
                                            aTCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGattatt                                               >  1:2620271/1‑91 (MQ=255)
                                             tCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGattatta                                              >  1:1175316/1‑91 (MQ=255)
                                              cTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGattattat                                             >  2:1788887/1‑91 (MQ=255)
                                              cTGATTTCATCGCTGCTTGGTTTGCAGGTAATACTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGattattat                                             <  1:3448172/91‑1 (MQ=255)
                                               tGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGattattatt                                            >  2:2386650/1‑91 (MQ=255)
                                                   ttCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCctt                                        >  1:940131/1‑91 (MQ=255)
                                                    tCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCctta                                       <  1:2653513/91‑1 (MQ=255)
                                                     cATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCcttac                                      <  2:592166/91‑1 (MQ=255)
                                                      aTCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCcttact                                     >  2:3511265/1‑91 (MQ=255)
                                                       tCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCcttactt                                    <  1:2174202/91‑1 (MQ=255)
                                                       tCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCcttactt                                    >  1:2872163/1‑91 (MQ=255)
                                                          ctgctTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTtacttact                                 <  1:1301310/91‑1 (MQ=255)
                                                          ctgctTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTtacttact                                 <  2:2182766/91‑1 (MQ=255)
                                                          ctgctTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTtacttact                                 >  2:1970728/1‑91 (MQ=255)
                                                           tgctTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTa                                >  1:2459996/1‑91 (MQ=255)
                                                            gctTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTAt                               <  1:2292410/91‑1 (MQ=255)
                                                             ctTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTAtt                              <  2:2079909/91‑1 (MQ=255)
                                                               tGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTa                            >  2:3045303/1‑91 (MQ=255)
                                                                 gTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTaaa                          >  2:3544360/1‑91 (MQ=255)
                                                                   ttGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTAAATc                        <  2:412149/91‑1 (MQ=255)
                                                                   ttGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTACTCCTTACTTACTATTTAAATc                        <  2:2532982/91‑1 (MQ=255)
                                                                    tGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTAAATCg                       <  1:3813435/91‑1 (MQ=255)
                                                                    tGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTAAATCg                       >  2:2258533/1‑91 (MQ=255)
                                                                    tGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTAAATCg                       >  1:2120079/1‑91 (MQ=255)
                                                                     gCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTAAATCGa                      >  1:633826/1‑91 (MQ=255)
                                                                     gCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTAAATCGa                      >  1:783835/1‑91 (MQ=255)
                                                                      cAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTAAATCGaa                     >  1:3844873/1‑91 (MQ=255)
                                                                       aGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTGTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTAAATCGAAt                    >  2:484997/1‑91 (MQ=255)
                                                                         gTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTAAATCGAATac                  >  1:1365040/1‑91 (MQ=255)
                                                                           aaTATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTAAATCGAATacac                >  2:933495/1‑91 (MQ=255)
                                                                           aaTATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTAAATCGAATacac                >  1:1899490/1‑91 (MQ=255)
                                                                           aaTATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTAAATCGAATacac                >  2:1200981/1‑91 (MQ=255)
                                                                            atatTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTAAATCGAATACACt               <  2:1680817/91‑1 (MQ=255)
                                                                                 tGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTAAATCGAATACACTAAata          >  2:2314843/1‑91 (MQ=255)
                                                                                    ggTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTAAATCGAATACACTAAAtatta       <  2:2247304/91‑1 (MQ=255)
                                                                                     ctGACGAACTGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTAAAACGAGTACACTAAAtattat      <  2:5438/90‑1 (MQ=255)
                                                                                      tGACAAATGGCCATTCGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTAAATCGAATACACTAAATattatt     >  1:192611/1‑91 (MQ=255)
                                                                                       gACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTAAATCGAATACACTAAATATTATTc    >  1:768621/1‑91 (MQ=255)
                                                                                        aCGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTAAATCGAATACACTAAATATTATTCa   >  2:1068305/1‑91 (MQ=255)
                                                                                        aCGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTAAATCGAATACACTAAATATTATTCa   <  2:726777/91‑1 (MQ=255)
                                                                                        aCGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTAAATCGAATACACTAAATATTATTCa   <  2:3306879/91‑1 (MQ=255)
                                                                                        aCGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTCTTGATTATTATTCATTACTTACTATTTAAATCGAATACACTAAATATAATTCa   >  1:3611779/1‑91 (MQ=255)
                                                                                          |                                                                                         
TATGGTGGCTCAATTACAATTATGTCAAAATTTGATGATAAGCAATCTGATTTCATCGCTGCTTGGTTTGCAGGTAATATTTGGGGTGACGAATGGCCATTTGTCATTGCATTTTTACCGTGGGTGTTGATTATTATTCCTTACTTACTATTTAAATCGAATACACTAAATATTATTCAT  >  USA300TCH1516_ALE/724362‑724541

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: