Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 479,096 T→C 5.5% F41F (TTT→TTC USA300HOU_RS02280 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE479,0960TC 5.5% 233.8 / 3.8 87F41F (TTT→TTCUSA300HOU_RS02280hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (38/41);  new base C (2/3);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.02e-01

TATAGGAGGTGAAAGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGATATGTTACTTGAAC  >  USA300TCH1516_ALE/478957‑479235
                                                                                                                                           |                                                                                                                                           
tataGGATGTGAAAGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTg                                                                                                                                            <  2:456326/141‑1 (MQ=255)
tataGGAGGTGAAAGGTCTGACAACACAGTTAAATATCCATTCAGTCGTTGAACATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTg                                                                                                                                            <  2:742715/141‑1 (MQ=255)
tataGGAGGTGAAAGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCGGCACGCAATCCATGGGAAGGATTCGAAGCGGATACGTTTg                                                                                                                                            >  1:336877/1‑141 (MQ=255)
tataGGAGGTGAAAGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTg                                                                                                                                            >  1:541615/1‑141 (MQ=255)
tagaggaggTGAAAGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTg                                                                                                                                            <  1:186537/138‑1 (MQ=255)
 ataGGAGGTGAAAGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGa                                                                                                                                           <  1:334142/141‑1 (MQ=255)
   aggaggTGAAAGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAg                                                                                                                                         <  1:106886/141‑1 (MQ=255)
    ggaggTGAAAGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTAAACCATTATCACTAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGa                                                                                                                                        >  2:663102/1‑141 (MQ=255)
     gaggTGAAAGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTGATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGAt                                                                                                                                       >  1:283623/1‑141 (MQ=255)
          gAAAGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGc                                                                                                                                  >  1:607933/1‑141 (MQ=255)
            tatgTATGTCAATACAGTTAATGATCAATTCAGTGATTGAAAATGCGAAGCGGGTCATTACACCATTATCACCAATTTCGATTTTTGCAGCACCAAATCCATGTGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCaa                                                                                                                                <  1:339361/138‑1 (MQ=255)
            aaGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCaa                                                                                                                                >  1:682083/1‑141 (MQ=255)
            aaGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATGAGAAGAGGATACGTTTGCAGATGTCGCaa                                                                                                                                >  2:421143/1‑141 (MQ=255)
             aGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGTATACATTTGAAGATGTCGCaaa                                                                                                                               <  2:564199/141‑1 (MQ=255)
             aGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCaaa                                                                                                                               >  1:162523/1‑141 (MQ=255)
             aGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGGGTTATTACACCATTATCACCAATTTCGATGTTTGCAGCACGAAATCCATGGGAAGGAATAGAAGCGGATACGTTTGAAGAGGTAGCaaa                                                                                                                               >  2:104455/1‑141 (MQ=255)
                 aTGACAACACGGTTAAATATCACTTGAGTCATTGAAAATGCGAAACGTGTTATTAGACCATTATCACCAAGTTCGTGTTTTGCAGCACGAAATCGATGGTAAGGATTAGAAGCGGATCGGTTTAAAGATGTCGGAAGATgg                                                                                                                           >  2:572119/1‑141 (MQ=255)
                 aTGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAGATgg                                                                                                                           >  2:572154/1‑141 (MQ=255)
                      aacaCAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACg                                                                                                                      <  1:684174/141‑1 (MQ=255)
                         ccccttAAACATCAATTCAGTCATCGAAAATGCGAAACGCGTTACTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGCCTGAAGATGTCGCAAAATGGTTACGTGa                                                                                                                   <  2:573358/137‑1 (MQ=255)
                         agaGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACTATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGa                                                                                                                   <  2:686430/139‑1 (MQ=255)
                          caGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATTGTTACGTGAt                                                                                                                  >  2:698809/1‑141 (MQ=255)
                            gTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGGTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACATTTGAAGATGTCGCAAAAGAGTTACGTGATGt                                                                                                                >  2:482486/1‑141 (MQ=255)
                              tAAACATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTc                                                                                                              <  1:89743/141‑1 (MQ=255)
                               aaaTATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGttgg                                                                                                             >  2:292631/1‑139 (MQ=255)
                                aaTATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCga                                                                                                            >  2:304540/1‑141 (MQ=255)
                                  tatCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCgaga                                                                                                          <  1:402837/141‑1 (MQ=255)
                                    tCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACAAATTTCGATTTTTGCAGCACGAAATCCATGGGCAGGATTAGTAGCGGGTACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGAtg                                                                                                        >  1:335010/1‑141 (MQ=255)
                                      aaTTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGGtgtg                                                                                                      >  2:498622/1‑141 (MQ=255)
                                       gttCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTgg                                                                                                     <  2:378528/140‑1 (MQ=255)
                                       aTTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTgg                                                                                                     <  2:246379/141‑1 (MQ=255)
                                       aTTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTgg                                                                                                     <  2:263736/141‑1 (MQ=255)
                                       aTTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTgg                                                                                                     <  2:495168/141‑1 (MQ=255)
                                            gtaaTTGAAAATGGGACACATGTTATTACACCATTATCACCAAGTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTGGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATAtt                                                                                                <  2:147165/138‑1 (MQ=255)
                                             tCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATAttt                                                                                               >  1:31908/1‑141 (MQ=255)
                                               aTTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGTATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATCtttcc                                                                                             >  2:448408/1‑139 (MQ=255)
                                               aTTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTc                                                                                             <  1:545254/141‑1 (MQ=255)
                                                ttGAAAATGCGAAACGTGTTATTCCACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGGTGTTCGAGAGGTGGATATTTTcc                                                                                            >  1:272813/1‑141 (MQ=255)
                                                ttGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTcc                                                                                            <  1:317150/141‑1 (MQ=255)
                                                 tGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTccc                                                                                           <  2:676959/141‑1 (MQ=255)
                                                  gAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCTATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCa                                                                                          <  1:158257/141‑1 (MQ=255)
                                                  gAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCa                                                                                          <  1:34859/141‑1 (MQ=255)
                                                    aaaTGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAGATCCATGGGAAGGATTAGAAGCGGATACGTTCGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCaaa                                                                                        >  2:248335/1‑141 (MQ=255)
                                                    aaaTGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAGATCCATGGGAAGGATTAGAAGCGGATACGTTCGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCaaa                                                                                        >  2:247916/1‑141 (MQ=255)
                                                         cGAAACGTGTTATTACACCATTATCACCAATGTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAGTGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAg                                                                                   >  2:433505/1‑141 (MQ=255)
                                                            aaCGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCAt                                                                                <  1:691927/141‑1 (MQ=255)
                                                              cGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTa                                                                              >  1:117968/1‑141 (MQ=255)
                                                                tgtTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAAt                                                                            <  1:458673/141‑1 (MQ=255)
                                                                    attaCACCGTTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAGTGGTTACGTGATGTTCGAGATGTAGATATTTTCCCAAATAAAGCATTAATAGaa                                                                        >  2:569724/1‑141 (MQ=255)
                                                                    attaCACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAAAAGTCACAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAgca                                                                        <  2:389939/141‑3 (MQ=255)
                                                                        cacCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAgag                                                                    >  2:315801/1‑139 (MQ=255)
                                                                         acCAGTATCACCAACTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAgcgc                                                                   <  1:317556/141‑1 (MQ=255)
                                                                             ttATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCtgtg                                                               >  1:108261/1‑141 (MQ=255)
                                                                             ttATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAAAAGAAAGCGCtgtg                                                               >  2:300947/1‑141 (MQ=255)
                                                                             ttATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGAGGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCtgtg                                                               >  1:607465/1‑141 (MQ=255)
                                                                              tATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTgg                                                              >  1:316354/1‑141 (MQ=255)
                                                                              tATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTgg                                                              >  2:494917/1‑141 (MQ=255)
                                                                              tATCACCAATTTCGATTTTTGCAGCACGAAATCCATGAGAAGGGTTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTgg                                                              >  1:29416/1‑141 (MQ=255)
                                                                                      aTTTCGATTTTTGCAGCACGAAATCCCTGGGAAGGGTTAGGAGCGGCTCTATTGGAAGATGTCGCAAAATGGTTACGTGATGTTCGGGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACgtggt                                                      <  1:758235/141‑1 (MQ=255)
                                                                                      aTTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATGAGAAGCGGGTTTTCGTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACgtggt                                                      <  1:758267/141‑1 (MQ=255)
                                                                                         tCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGaa                                                   >  1:279715/1‑141 (MQ=255)
                                                                                           gATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGAGGTGAAtt                                                 >  1:228351/1‑141 (MQ=255)
                                                                                            aTTTTTGCAGAACGAAATCCATGGGAAGGATTAGAAGCGGATAGGTTTGAAGATGTGGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCAGGTGATGAATTa                                                >  2:300504/1‑141 (MQ=255)
                                                                                               tttGCAGCACGAGATCCATGGGAAGGATTAGAAGCGGATACGTTCGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGAt                                             <  1:247916/141‑1 (MQ=255)
                                                                                               tttGCAGCACGAGATCCATGGGAAGGATTAGAAGCGGATACGTTCGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGAt                                             <  1:248335/141‑1 (MQ=255)
                                                                                               tttGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGAt                                             <  1:485593/141‑1 (MQ=255)
                                                                                                 tGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGa                                           <  2:232873/141‑1 (MQ=255)
                                                                                                 tGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGa                                           >  1:755246/1‑141 (MQ=255)
                                                                                                 tGCAGCACGAAATCCATGAGAAGGCTTAAAAACGGATCCAATAGCAGATGTCACAAAATGATTATGAGATACTCGAGAAATGTATAGTTGCCCAAGCAAAGCCTTAAGAGAAAGCGCGTTGGCACGTGGTGAATTAGATGa                                           >  2:181358/1‑141 (MQ=255)
                                                                                                    agaaCGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGGTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCTCTGTGGCACGTGGTGAATTAGATGAAAg                                        <  2:335010/138‑1 (MQ=255)
                                                                                                      caCGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAgtg                                      <  1:745197/141‑1 (MQ=255)
                                                                                                               caggggAAGGATTAGGAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATc                             <  2:301155/138‑1 (MQ=255)
                                                                                                               cATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGCAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATc                             <  1:744336/141‑1 (MQ=255)
                                                                                                                      agtaTTAGAAGCGGATACGTTTGAAGATGTCGTAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGt                      <  1:330938/138‑1 (MQ=255)
                                                                                                                       ggATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGtt                     >  2:719053/1‑141 (MQ=255)
                                                                                                                           tAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTg                 <  2:633197/141‑1 (MQ=255)
                                                                                                                           tAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTg                 <  1:114250/141‑1 (MQ=255)
                                                                                                                              gaGCGGATACGTTTGAAGATGTCGAAAAATGGTTACGTGTTGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGata              <  1:131148/140‑1 (MQ=255)
                                                                                                                                 gggATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGCTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGTACGGTGTGAATTAGATGAAAGTGCATTTAATCAACTTGTTACTGATATGt           <  2:23068/140‑1 (MQ=255)
                                                                                                                                 cGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGATATGt           <  1:609660/141‑1 (MQ=255)
                                                                                                                                 cGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGATATGt           <  1:421143/141‑1 (MQ=255)
                                                                                                                                  ggATACGTTTGAAGGTGTCGCAAAATGGTGACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGATATGtt          <  2:39068/141‑1 (MQ=255)
                                                                                                                                   gATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGATATGTTa         <  2:206159/141‑1 (MQ=255)
                                                                                                                                     tACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGATATGTTaat       >  2:725948/1‑139 (MQ=255)
                                                                                                                                     tACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGATATGTTACt       >  1:304581/1‑141 (MQ=255)
                                                                                                                                       cGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGATATGTTACTTg     >  2:736363/1‑141 (MQ=255)
                                                                                                                                        gTTCGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGATATGTTACTTGa    <  1:45913/141‑1 (MQ=255)
                                                                                                                                          ttGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGATATGTTACTTGAAc  <  2:607465/141‑1 (MQ=255)
                                                                                                                                          ttGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGATATGTTACTTGAAc  <  2:108261/141‑1 (MQ=255)
                                                                                                                                          ttGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTACCAAAGAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGCTACTGATATGTTGCTTGaaa  >  2:600294/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                           
TATAGGAGGTGAAAGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGATATGTTACTTGAAC  >  USA300TCH1516_ALE/478957‑479235

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: