Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,744,424 A→G 11.2% L83S (TTA→TCA)  apt ← Adenine phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,744,4240AG11.2% 193.2 / 14.3 76L83S (TTA→TCA) aptAdenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (35/32);  new base G (6/3);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

TTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATT  >  USA300TCH1516_ALE/1744288‑1744555
                                                                                                                                        |                                                                                                                                   
ttCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAAttt                                                                                                                                 >  2:854842/1‑141 (MQ=255)
  cAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTcc                                                                                                                               <  2:889085/141‑1 (MQ=255)
   aaTCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGCCCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTccc                                                                                                                              <  1:630704/141‑1 (MQ=255)
     tCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTGATTTCCctt                                                                                                                            <  2:351253/141‑1 (MQ=255)
     ccGTACCACCAGTAGCTAATAAATCATCTGTAATTAACTCACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCctt                                                                                                                            <  2:575941/140‑1 (MQ=255)
       gTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCcttct                                                                                                                          >  2:683442/1‑141 (MQ=255)
       atACCACCAGTAGCTAGTAAATTATATGTAATTAACACACGTTGACCTGGTTTAATTGACTATTCGTGAATTGTTAATACGTTTGTACCATATTCTATGTCATACTCATAACGAATGACTTCACGAGGTGATTTCCcttct                                                                                                                          >  2:280350/2‑141 (MQ=255)
        tACCACCAGTATCTAATAAATCATCTGTAATTAATACACTTGAACGTGGTTGAATTTCCTCTTTTTGCATTAATATACCGTGTTTAGCATATTCTACGACCTTCTACTCACGAATGACTTCACGAGTTAATTTCtcttctg                                                                                                                         <  2:418249/141‑2 (MQ=255)
         accaccAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTGATTTCCCTTCttt                                                                                                                        <  2:741201/141‑1 (MQ=255)
         accaccAGTAGCTAATAAATCATACGTAATTAACACACGTTGACCTGCTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTGATTTCCCTTCttt                                                                                                                        <  2:864309/141‑1 (MQ=255)
              caGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTac                                                                                                                   <  2:407857/141‑2 (MQ=255)
               aGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAAc                                                                                                                  <  1:872725/141‑1 (MQ=255)
               aGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAAc                                                                                                                  <  2:409745/141‑1 (MQ=255)
                 tAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAg                                                                                                                <  2:315548/141‑1 (MQ=255)
                  aGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAgg                                                                                                               <  1:427027/141‑1 (MQ=255)
                   gCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGt                                                                                                              >  1:262485/1‑141 (MQ=255)
                       ataaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCaa                                                                                                          <  1:416044/141‑1 (MQ=255)
                                 tGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATccc                                                                                                <  2:148339/141‑1 (MQ=255)
                                   tAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCa                                                                                              >  1:943703/1‑141 (MQ=255)
                                   tAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCa                                                                                              >  1:869307/1‑141 (MQ=255)
                                      ttAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTg                                                                                           <  2:135982/141‑1 (MQ=255)
                                        aaCACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGATGCAAAGCCAATCCCCATTGaa                                                                                         >  1:1000887/1‑141 (MQ=255)
                                        aaCACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCATTTCTAACAGATGCAAAGCCAATCCCCATTGaa                                                                                         >  2:619287/1‑141 (MQ=255)
                                        aaCACACGTTGACCTGGGTGAATTGCCTCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGaa                                                                                         >  2:286288/1‑141 (MQ=255)
                                                 tGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACa                                                                                >  1:29181/1‑141 (MQ=255)
                                                  gACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTGATTTCCCTTCTTTTCTAACAGGTGCAACGCGCATGCCCATTGAATAAGTTACAg                                                                               >  1:40698/1‑141 (MQ=255)
                                                  gACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAg                                                                               >  2:1046529/1‑141 (MQ=255)
                                                  gACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATGTAATAAGCTACAg                                                                               >  2:616257/1‑141 (MQ=255)
                                                   aCCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAgg                                                                              <  1:986452/141‑1 (MQ=255)
                                                   aCCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAgg                                                                              <  1:114995/141‑1 (MQ=255)
                                                     cTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTATCAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGAc                                                                            <  1:491287/141‑1 (MQ=255)
                                                      tGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACa                                                                           <  1:863891/141‑1 (MQ=255)
                                                      tGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACa                                                                           <  2:1035160/141‑1 (MQ=255)
                                                              ttGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACCGGTGCAAAGCCAATCCCCATTGAATAGACTACGTGCCAGCCaagaa                                                                   >  2:873600/1‑138 (MQ=255)
                                                              ttGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGa                                                                   <  2:971475/141‑1 (MQ=255)
                                                                 cATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTTATACGCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATaa                                                                >  1:759984/1‑141 (MQ=255)
                                                                 cATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGATGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATaa                                                                <  2:464928/141‑1 (MQ=255)
                                                                 cATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATaa                                                                <  2:450656/141‑1 (MQ=255)
                                                                           cATTGTTAAAACATTTGTACCATGGTCGAGGTCATACTCATAACGAATGACTTCACGAGGTGATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACgcgc                                                      >  2:612905/1‑141 (MQ=255)
                                                                                 tAAAACATTTGTACCATATACTACATAAGGTTCATAACGAATGACTTCACGAGGTGATTTACATTCTTTTCTAACAGATGCAGAGCCAATCCCCATTGAATAAGTTAAAGCACAGCCAATGATAAAGCAACGCGATTAAgg                                                >  2:507793/1‑141 (MQ=255)
                                                                                    aaCATTTGTACCATATTCTAGGTCATAAGCATAACGAATGACTTCAAGAGCTGAGTTCCCTGCTTTTCTAAAAGGTGCAAAGCAGATCCCCCTTGCATAAGCTACAGGACAGCCAATGATCAAGCCCACCGCTTCAGGTcc                                             >  2:351759/1‑141 (MQ=255)
                                                                                      cATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTCCAGGTCCTa                                           >  1:728075/1‑141 (MQ=255)
                                                                                      cATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAAGCGCTACAGGACCTa                                           >  1:1026168/1‑141 (MQ=255)
                                                                                         ttGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCACTTCTTTTCTAACAGGTGCACAGCCAATCCCCATTGAATAAGCTACAGAACAGCCAATGATAAGGCCACGCGCTTCAGGCCCTACaa                                        >  2:31363/1‑141 (MQ=255)
                                                                                            tACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACgt                                     >  1:884359/1‑140 (MQ=255)
                                                                                            tACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCACATTGAATAAGCTACAGGACGGCCACTGATAAAGCCACGCGCTTCAGGTCCTACAACGa                                     >  2:1003896/1‑141 (MQ=255)
                                                                                            tACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCACATTGAATAAGCTACAGGACAGCCAATGCTAAAGCCACGCGCTTCACGTCCTCCAACGa                                     >  1:539519/1‑141 (MQ=255)
                                                                                                atatTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGTCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATc                                 <  1:223502/141‑1 (MQ=255)
                                                                                                atatTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATc                                 <  1:922729/141‑1 (MQ=255)
                                                                                                atatTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATc                                 <  1:167476/141‑1 (MQ=255)
                                                                                                agatTCTAGGTCATACTCATGACGAATGCCTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGCATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATc                                 <  1:997631/139‑1 (MQ=255)
                                                                                                acaatgtAGGGCAGACTGATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATc                                 <  1:797533/135‑1 (MQ=255)
                                                                                                 tatTCTAGGTCATACTCATAACGAATGACTTCACGAGGTGATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACCACGATATCa                                >  1:958387/1‑141 (MQ=255)
                                                                                                   ttCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAAc                              <  2:869307/141‑1 (MQ=255)
                                                                                                     cTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAt                            >  1:899848/1‑141 (MQ=255)
                                                                                                     cTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAt                            >  1:991812/1‑141 (MQ=255)
                                                                                                     cTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATAAACAt                            >  1:122716/1‑141 (MQ=255)
                                                                                                        ggTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtctc                         >  2:763323/1‑141 (MQ=255)
                                                                                                        ggTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAATTCCGACAACAAAACCAACGtctc                         >  2:1110513/1‑141 (MQ=255)
                                                                                                         gTCATACTCATAACCAAAGTCATGACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGAAAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATctct                        >  2:786192/1‑141 (MQ=255)
                                                                                                            aTACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACCCGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTc                     <  1:31363/141‑1 (MQ=255)
                                                                                                            aTACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTc                     >  2:693736/1‑141 (MQ=255)
                                                                                                             tACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTCACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGGCCCGCTCTTCAGGTAATAAAAAAATAAAAACATCTCTGTCt                    <  2:454733/141‑1 (MQ=255)
                                                                                                                 cATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGCTTTCAACATCTCTGTCTTTTg                >  2:174568/1‑141 (MQ=255)
                                                                                                                 cATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTg                >  1:663634/1‑141 (MQ=255)
                                                                                                                 cATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGACTTTTg                >  2:537557/1‑141 (MQ=255)
                                                                                                                 cATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTAACTAAGCTACAGGACAGCCAATGATAAAGAAACGCGCTTCATGTCCTACAACGTTATCAACATATCTGTCTTTTg                >  2:214301/1‑141 (MQ=255)
                                                                                                                 cATAACGAATGACTTCACGAGGTAATTTCCCTTATTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCGTCAGGTCCTACACCGATCTCACCACCCGTGTCTATTg                >  2:216435/1‑141 (MQ=255)
                                                                                                                   tAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAAGGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCg              >  2:315942/1‑141 (MQ=255)
                                                                                                                      cGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTAt           >  1:285698/1‑141 (MQ=255)
                                                                                                                         aTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAGTCCCCATCGAATAAGCTACAGGACATCCAATGATAAAGCCACGCGCTTCAGCTCCTACAACGATCTCTAATGCCCTGTCTTTTGCGTATTCt        <  1:814707/141‑1 (MQ=255)
                                                                                                                         aTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCt        <  2:29181/141‑1 (MQ=255)
                                                                                                                         aTGACTGCACGAGGTAATTTCCCTTCGTTTCTAACAGTTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCt        <  2:264397/141‑1 (MQ=255)
                                                                                                                           gACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCAGGGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTAc      >  1:767257/1‑141 (MQ=255)
                                                                                                                               tCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAAtt  <  2:46622/141‑1 (MQ=255)
                                                                                                                               tCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGACTTTTGCGCATTCTACAAtt  >  1:556177/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                   
TTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATT  >  USA300TCH1516_ALE/1744288‑1744555

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: ‑