Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 964,580 T→C 9.0% F709L (TTT→CTT)  addB → ATP‑dependent helicase/deoxyribonuclease subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE964,5800TC 9.0% 167.8 / 2.2 64F709L (TTT→CTT) addBATP‑dependent helicase/deoxyribonuclease subunit B
Reads supporting (aligned to +/- strand):  ref base T (22/36);  new base C (3/3);  total (25/39)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGATAAGGAGATCAGCCCATTTTTAAATCAAATTCAATCATTGTTCAACCAATTGGAAATTACTAACATTCCTCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGATTATTTAATGTCAGCATTAACGTTTGACAATGAA  >  USA300TCH1516_ALE/964442‑964717
                                                                                                                                          |                                                                                                                                         
gatgatAAGGAGATCAGCCCATTTTTAAATCAAATTCAATCATTGTTCAACCAATTGGAAATTACTAACATTCCTCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTAttt                                                                                                                                         >  1:736971/1‑141 (MQ=255)
 atgatAAGGAGATCAGCCCATTTTTAAATCAAATTCAATCATTGTTCAACCAATTGGAAATTACTAACATTCCTCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTg                                                                                                                                        <  2:124546/141‑1 (MQ=255)
   gatAAGGAGATCACCCCATTTTTAAATCAAATACAACCCTTGTACAAACAATTGTAAATTACTAACATTACTCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGaa                                                                                                                                      <  1:61140/141‑1 (MQ=255)
        ggAGATCAGCCTATTTTTAAATCAAATTCAATCATTGTTCGACCAATCGTGATTGACTAACATTCCTCTATACCATGAAGTTAACCCATTGTCACGAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCAtt                                                                                                                                 <  2:323040/141‑1 (MQ=255)
        ggAGATCAGCCCATTTTTAAATCAAATTCAATCATTGTTCAACCAATTGGAAATTACTAACATTCCTCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCAtt                                                                                                                                 >  1:486835/1‑141 (MQ=255)
          gtaTCAGCCCATTTTTAAATCAAATTCAATCATTGTTCAACCAATTGGAAATTACTAACATTCCTCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGc                                                                                                                               <  1:127778/139‑1 (MQ=255)
          agaTCAGCCCATTTTTAAATCAAATTCAATCATTGTTCATCCAATTGGAAATTACTAACATTCCTCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGc                                                                                                                               <  2:333829/141‑1 (MQ=255)
          agaTCAGCCCATTTTTAAATCAAATTCAATCATTGTTCAACCAATTGGAAATTACTAACATTCCTCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGc                                                                                                                               >  1:107001/1‑141 (MQ=255)
           gaTCAGCCCATTTTTAAATCAAATTCAATCATTGTTCAACCAATTGGAAATTACTAACATTCCTCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCg                                                                                                                              <  2:504070/141‑1 (MQ=255)
                gCCCATTTTTAAATCAAATTCAATCATTGTTCAACCAATTGGAAATTACTAACATTCCTCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCtt                                                                                                                         >  2:228396/1‑141 (MQ=255)
                 tccATTTTTAAATCAAATTCAATCATTGTTCAACCAATTGGAAATTACTAACATTCCTCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTg                                                                                                                        <  1:340530/140‑1 (MQ=255)
                  ccATTTTTAAATCAAATTCAATCATTGTTCAACCAATTGGCAATTACTAACATTCCTCAATACCCTGAATTTAGCCCATTGTCACTAATGCAACATTCTAAACAAACCAAAATTACATTACTTGAAGCATTGCGGGCTTgt                                                                                                                       >  1:121221/1‑140 (MQ=255)
                  ccATTTTTAAATCAAATTCAATCATTGTTCAACCAATTGGAAATTACTAACATTCCTCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTACTTGAAGCATTGCGTGCTTgg                                                                                                                       >  1:55058/1‑141 (MQ=255)
                  ccATTTTTAAATCAAATTCAATCATTGTTCAACCAATTGGAAATTACTAACATTCCTCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTACTTGAAGCATTGCGTGCTTgg                                                                                                                       >  1:750684/1‑141 (MQ=255)
                     tttttAAATCAAATTCAATCATTGTTCAACCAATTGGAAATTACTAACATTCCTCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTta                                                                                                                    >  1:683535/1‑141 (MQ=255)
                      ttttAAATCAAATTCAATCATTGTTCAACCAATTGGAAATTACTAACATTCCTCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTtat                                                                                                                   <  1:632842/141‑1 (MQ=255)
                          aaaTCAAATTCAATCATTGTTCAACCAATTGGAAATTACTCACATCCCCCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATatga                                                                                                               <  1:735544/141‑1 (MQ=255)
                          aaaTCAAATTCAATCATTGTTCAACCAATTGGAAATTACTAACATTCCTCAATACCATGAAGTGAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTGATatga                                                                                                               >  1:604664/1‑141 (MQ=255)
                           aaTCAAATTCAATCATTGTTCAACCAATTGGAAATTACTAACATTCCTCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATatgat                                                                                                              <  2:311368/141‑1 (MQ=255)
                                             ttaaagcAATTGGAATTTACTAACTATCCTCCATACCATGAAGTTAGCCCATTGTGACGCATGCAACATGCTAAGCAAACCAACATTACAGTATTTGAAGCATGGCGTTCTTGGTTATATGATGAAATTGTGGCTGATAGt                                                                                            <  2:327399/135‑1 (MQ=255)
                                              tCAACCAATTGGAAATTACTAACATTCCTCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCAGTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGtt                                                                                           <  2:340609/141‑1 (MQ=255)
                                                 aCCCATTGGAAATTACTAACATTCCTCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGt                                                                                        <  1:121718/141‑1 (MQ=255)
                                                      ttGGAAATTACTAACATTCCTCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTTTGGCTGATAGTTGGTTAGAt                                                                                   <  1:480664/141‑1 (MQ=255)
                                                          aaaTTACTAACATTCCTCAATTCCATGAAGTTAACTCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACACTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCtt                                                                               <  1:589836/141‑1 (MQ=255)
                                                                 tAACATTCCCCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGt                                                                        <  1:392483/141‑1 (MQ=255)
                                                                  aaCATTCCTCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTa                                                                       >  1:79655/1‑141 (MQ=255)
                                                                       tCCTCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATCGTTGGTTAGCTGCCTATCAAGTAATTag                                                                  >  2:122194/1‑141 (MQ=255)
                                                                       tCCTCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTCATCAAGTAATTag                                                                  >  2:567336/1‑141 (MQ=255)
                                                                         cTCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATACTTATCAAGTAATTagaa                                                                >  2:26066/1‑140 (MQ=255)
                                                                            aaTACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTATGGATa                                                             >  2:453241/1‑141 (MQ=255)
                                                                             aTACCATGAAGTTAACCCATTGTCACTCATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAg                                                            >  2:133344/1‑141 (MQ=255)
                                                                                ccATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGa                                                         <  2:736971/141‑1 (MQ=255)
                                                                                        ttAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTaa                                                 <  1:1988/141‑1 (MQ=255)
                                                                                         tAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTaaa                                                >  2:666104/1‑141 (MQ=255)
                                                                                         tAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTaaa                                                >  2:122594/1‑141 (MQ=255)
                                                                                          aaCCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGCGACTGATAGTTGGTTAGATGCATATCAAGTAATTGGAGATAGCGATCATTTAAAt                                               >  2:432368/1‑141 (MQ=255)
                                                                                             ccAGTGTCCCTAATGCAACATGGTAAGCAAACCAAAATTACATTATTCGAAGCATGGCGTGCTTGGTTATATGATGAAATTGTGGCTGTACGCTGGGTGGCGGCTTATCAGCTAATTAGAGATAGCGAACGTTTAGAGCaa                                            >  2:120210/1‑141 (MQ=255)
                                                                                              cATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAg                                           <  2:309765/141‑1 (MQ=255)
                                                                                                 tGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGTGATCATTTAAATCAAGGtt                                        <  1:473291/141‑1 (MQ=255)
                                                                                                  gTCACTAACGCAACATGCTAAGCAAACCAAAATTACATTATTTTAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGATATAGCGATCATTTAAATCAAGGttt                                       <  1:366412/141‑1 (MQ=255)
                                                                                                    cACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAg                                     >  2:258984/1‑141 (MQ=255)
                                                                                                     aCTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGa                                    >  2:641078/1‑141 (MQ=255)
                                                                                                     aCTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGa                                    <  2:680058/141‑1 (MQ=255)
                                                                                                     aCTAAGTCAACATGCTAAGCAAACCAAAATTACAGTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATTAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGa                                    <  1:44141/141‑1 (MQ=255)
                                                                                                      cTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTCATTAGAGATAGCGATCATTTAAATCAAGGTTTAGat                                   >  1:607787/1‑141 (MQ=255)
                                                                                                         aTGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGattat                                <  1:228396/141‑1 (MQ=255)
                                                                                                               cATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGATTATTTAATg                          >  1:667850/1‑141 (MQ=255)
                                                                                                                   cTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGATTATTTAATGTCAg                      <  1:722786/141‑1 (MQ=255)
                                                                                                                    tAAGCAAACCCAAATTACATTATTGGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGATTATTTAATGTCAGc                     <  2:261069/141‑1 (MQ=255)
                                                                                                                    tAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGATTATTTAATGTCAGc                     <  2:438117/141‑1 (MQ=255)
                                                                                                                      aGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAATTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGATTATTTAATGTCAGCAt                   <  1:508848/141‑1 (MQ=255)
                                                                                                                      aGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGATTATTTAATGTCAGCAt                   <  2:482424/141‑1 (MQ=255)
                                                                                                                      aGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGATTATTTAATGTCAGCAt                   <  1:130868/141‑1 (MQ=255)
                                                                                                                        cAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGATTATTTAATGTCAGCATTa                 <  1:461402/141‑1 (MQ=255)
                                                                                                                          aaCCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGATTATTTAATGTCAGCATTAAc               <  2:196170/141‑1 (MQ=255)
                                                                                                                           aCCAAAATTACATTACTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGATTATTTAATGTCAGCATTAACg              <  2:55058/141‑1 (MQ=255)
                                                                                                                           aCCAAAATTACATTACTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGATTATTTAATGTCAGCATTAACg              <  2:750684/141‑1 (MQ=255)
                                                                                                                           aCCAAAATTACATTACTGGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGATTATTTAATGTCAGCATTAACg              <  2:121221/141‑1 (MQ=255)
                                                                                                                             cAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTTTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGATTATTTAATGTCAGCATTAACGtt            >  2:394066/1‑141 (MQ=255)
                                                                                                                             cAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGATTATTTAATGTCAGCATTAACGtt            >  2:194641/1‑141 (MQ=255)
                                                                                                                             cAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGATTATTTAATGTCAGCATTAACGtt            >  1:34711/1‑141 (MQ=255)
                                                                                                                               aaaTTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGATTATTTAATGTCAGCATTAACGtttt          <  2:134871/141‑2 (MQ=255)
                                                                                                                                 aTTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGATTATTTAATGTCAGCATTAACGTTTGAc        <  1:77648/141‑1 (MQ=255)
                                                                                                                                    aCATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGATTATTTAATGTCAGCATTAACGTTTGACAAt     <  2:25613/141‑1 (MQ=255)
                                                                                                                                     cATTATTTGAGGCATTGCGGGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGATTATTTAATGTGAGCATTAACGGTTGTCAATg    <  2:267109/141‑1 (MQ=255)
                                                                                                                                       ttattTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGATTATTTAATGTCAGCATTAACGTTTGACAATCaa  >  2:753439/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                         
GATGATAAGGAGATCAGCCCATTTTTAAATCAAATTCAATCATTGTTCAACCAATTGGAAATTACTAACATTCCTCAATACCATGAAGTTAACCCATTGTCACTAATGCAACATGCTAAGCAAACCAAAATTACATTATTTGAAGCATTGCGTGCTTGGTTATATGATGAAATTGTGGCTGATAGTTGGTTAGATGCTTATCAAGTAATTAGAGATAGCGATCATTTAAATCAAGGTTTAGATTATTTAATGTCAGCATTAACGTTTGACAATGAA  >  USA300TCH1516_ALE/964442‑964717

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