Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,963,768 C→G 90.1% V365L (GTA→CTA)  hemY ← Protoporphyrinogen oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,963,7680CG90.1% 331.8 / 26.0 111V365L (GTA→CTA) hemYProtoporphyrinogen oxidase
Reads supporting (aligned to +/- strand):  ref base C (8/3);  new base G (43/57);  total (51/60)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.69e-01

Only 100 of 111 total aligned reads displayed.

ACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCTT  >  USA300TCH1516_ALE/1963631‑1963905
                                                                                                                                         |                                                                                                                                         
ccctttCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATa                                                                                                                                        <  1:317391/138‑1 (MQ=255)
acgtttCAGCAAACTACTGACAAGTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATa                                                                                                                                        <  1:671460/138‑1 (MQ=255)
aCTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATa                                                                                                                                        <  1:650050/141‑1 (MQ=255)
  ttttCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACg                                                                                                                                      >  1:842980/1‑141 (MQ=255)
   tttCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGc                                                                                                                                     >  1:378268/1‑141 (MQ=255)
      cGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACAAAAGTATCACCTGGTTTACCTAGATCCGCACg                                                                                                                                  >  2:430007/1‑141 (MQ=255)
       ggCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATATCATCATTTCACTTGAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACAACAGGATCACCTGGTTTACCTAGATACGCACGa                                                                                                                                 <  2:330998/141‑1 (MQ=255)
         cAAACGCTTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAAt                                                                                                                               <  2:393144/141‑1 (MQ=255)
          aaaCGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATc                                                                                                                              <  1:490355/141‑1 (MQ=255)
              gATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCaaaa                                                                                                                          >  1:610900/1‑141 (MQ=255)
                 tGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGCTCATCAACCACAGTACCACCTGGTTTACCTAGATACGCACGCATCAAAacct                                                                                                                       >  2:257238/1‑141 (MQ=255)
                 tGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCGCCTGGTTTACCTAGATACGCACGAATCAAAacct                                                                                                                       >  1:948022/1‑141 (MQ=255)
                 tGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAacct                                                                                                                       >  1:457681/1‑141 (MQ=255)
                    cAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTCCGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACATGGGTTACCTAGATACGCACGAATCAAAacctccc                                                                                                                    >  1:979058/1‑138 (MQ=255)
                    cAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGTTTTACCTAGATACGCACGAATCAAAaccttac                                                                                                                    >  1:697907/1‑141 (MQ=255)
                    cAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCATGAATCAAAaccttac                                                                                                                    <  1:786395/141‑1 (MQ=255)
                    cAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAaccttac                                                                                                                    <  2:684085/141‑1 (MQ=255)
                    cAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAaccttac                                                                                                                    >  2:713821/1‑141 (MQ=255)
                    cAATTGTAAATTCAGGATCACCTTTAAAGTTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAaccttac                                                                                                                    <  1:67497/141‑1 (MQ=255)
                    aaaTTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAaccttac                                                                                                                    <  2:445740/140‑1 (MQ=255)
                      aTTGTAACATCAGGATCACCTTTAAATGTCATCCTTTGACGTAAATCTCTACGTACATTCGATACGAATTCATTATCTGTATGATCATCATCCACTGTATCACCTGGTTTACCTTGATACGCACGAATTAAAaccttacct                                                                                                                  <  2:335794/141‑1 (MQ=255)
                      aTTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCATCTGGTTTACCTAGATACGTACTAATCAAAaccttacct                                                                                                                  <  1:592511/141‑1 (MQ=255)
                      aTTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAaccttacct                                                                                                                  <  1:522931/141‑1 (MQ=255)
                         gTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTcc                                                                                                               <  2:136015/141‑1 (MQ=255)
                            aaTTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGgt                                                                                                            >  2:563063/1‑141 (MQ=255)
                             aTTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCGCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGgtg                                                                                                           <  1:454384/141‑1 (MQ=255)
                              ttCAGGATCACCTTTAAATGTCCTCATTTGACTTAAACCCCTACGTACAACCGATACTAATTCATTATCTGTATGATCATCAACCACATTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGgtgt                                                                                                          <  1:315134/141‑1 (MQ=255)
                              ttCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATAACCCGAATCCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGgtgt                                                                                                          >  1:771797/1‑141 (MQ=255)
                                 aGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTgtagt                                                                                                       >  2:834181/1‑141 (MQ=255)
                                  ggATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGGgtagta                                                                                                      >  1:887223/1‑141 (MQ=255)
                                   gATCACCTTTAAATGTCATCCTTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAACCAAAACCTTACCTTCCGGTGTAGTaa                                                                                                     >  1:641613/1‑141 (MQ=255)
                                    aTCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATAACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTaaa                                                                                                    >  1:1015828/1‑141 (MQ=255)
                                      cACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTCCCGGTGGATGAAATg                                                                                                  >  1:134501/1‑141 (MQ=255)
                                         cTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGcc                                                                                               >  1:148272/1‑141 (MQ=255)
                                           ttAAATGTCATCCTTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCAt                                                                                             <  1:432619/141‑1 (MQ=255)
                                             aaaTGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCAttt                                                                                           <  1:919095/141‑1 (MQ=255)
                                                  ttatcatTTGACTTAAATCTCTACGTACAATCGATACTACTGCATTATCTGTATGATCCGCAACCACAGTATCACCTGTTTTACCTAGATACGCACGAATCAAAACCTTAGCTTCCGGTGTAGTAAATGGCCATTTTTTCg                                                                                      <  2:377733/139‑1 (MQ=255)
                                                  tcatcaTTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGCATGACCATCAAACACAGTATCACCTGGGTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCg                                                                                      <  2:377905/141‑1 (MQ=255)
                                                  tcatcaTTTGACTTAAATCTCGACGAACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCg                                                                                      <  1:746827/141‑1 (MQ=255)
                                                   tatcatTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGa                                                                                     <  1:844676/140‑1 (MQ=255)
                                                   catcatTTGACTTAAATCTCTCCGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGa                                                                                     <  1:843670/141‑1 (MQ=255)
                                                   catcatTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGa                                                                                     <  1:502026/141‑1 (MQ=255)
                                                   catcatTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGa                                                                                     <  1:731118/141‑1 (MQ=255)
                                                   catcatTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGa                                                                                     <  1:609096/141‑1 (MQ=255)
                                                    atcatTTGACTTAAACCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGAt                                                                                    <  1:656444/141‑1 (MQ=255)
                                                      catTTGACTTAAATCTCTACGTACAATCGATACTAATTCCTACACTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGt                                                                                  >  1:728947/1‑141 (MQ=255)
                                                      catTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGt                                                                                  >  1:107386/1‑141 (MQ=255)
                                                         ttGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCa                                                                               >  2:911336/1‑141 (MQ=255)
                                                         ttGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTACGGTGTAGTAAATGGCCATTTTTTCGATGTCCa                                                                               >  2:592052/1‑141 (MQ=255)
                                                           gACTTAACTCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAg                                                                             <  2:854659/141‑1 (MQ=255)
                                                           gACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAg                                                                             >  1:1008751/1‑141 (MQ=255)
                                                           gACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAg                                                                             <  1:851897/141‑1 (MQ=255)
                                                           gACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGAGGTCCAAg                                                                             >  2:975172/1‑141 (MQ=255)
                                                           gACTCAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAg                                                                             <  2:638379/141‑1 (MQ=255)
                                                            ccTTCAAGCTCCACGTACAATCGAGACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAACTGGCCATTTTTTCGATGTCCAAGt                                                                            <  1:470095/140‑1 (MQ=255)
                                                               tAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTACGATGTCCAAGTACa                                                                         <  1:789878/141‑1 (MQ=255)
                                                                 aaTCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATg                                                                       >  2:4320/1‑141 (MQ=255)
                                                                 aaTCTCTACGTACAATCGATACTAATTCATTATCCTTATGATCATCAACCACAGTATCCCCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATg                                                                       <  2:684075/141‑1 (MQ=255)
                                                                  aTCTCTACGTACAAGCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGc                                                                      <  1:486894/141‑1 (MQ=255)
                                                                    ctctCCGTCTAATCGCTACTAGTTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCgg                                                                    <  2:395024/141‑1 (MQ=255)
                                                                      ctccGTACCATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTCGATACGCACGAATCAACAACTGGCCTGCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTa                                                                  <  1:508773/138‑1 (MQ=255)
                                                                         cGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCATGTACATGCGGTAATg                                                               >  2:285716/1‑141 (MQ=255)
                                                                           tACAATCGATACTCATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTc                                                             <  1:929447/141‑1 (MQ=255)
                                                                           tACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTc                                                             <  2:907633/141‑1 (MQ=255)
                                                                               aTCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTa                                                         <  1:655255/141‑1 (MQ=255)
                                                                                ccGATACTAATTCATTTTCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTCCATGCGGTAATGTCTGTAt                                                        <  2:534052/140‑1 (MQ=255)
                                                                                            cATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCTTTTTTTCGATGCCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCg                                            <  2:435160/141‑1 (MQ=255)
                                                                                            cATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCg                                            >  1:928542/1‑141 (MQ=255)
                                                                                                aTCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCTTTTTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCCCAAt                                        >  2:639011/1‑141 (MQ=255)
                                                                                                aTCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCGAGTGCATGCGGTACTGTCTGTATCACTGGTTCTGCCAAt                                        >  1:261050/1‑141 (MQ=255)
                                                                                                aTCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAAt                                        >  1:616378/1‑141 (MQ=255)
                                                                                                aTCTGTATGATCATCAACCACAGTATCACCTGGGTTACCTAGATACGCACGAATCAAGACCTTCCCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGCTAGTGTCTGTATCCCTCGTTCTCGCAAt                                        >  1:490475/1‑141 (MQ=255)
                                                                                                  cTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTa                                      <  1:901551/141‑1 (MQ=255)
                                                                                                      aTGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACgca                                  >  1:648954/1‑139 (MQ=255)
                                                                                                        gATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGc                                <  1:640286/141‑1 (MQ=255)
                                                                                                         atcatcAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGcc                               >  1:909186/1‑141 (MQ=255)
                                                                                                         atcatcAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGcc                               <  1:1022888/141‑1 (MQ=255)
                                                                                                          tcatcaACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCa                              >  1:974144/1‑141 (MQ=255)
                                                                                                            atcaACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCTATTTCCAACTACATGCGGTAATGTCACTATCACTCGTACTCGCAATTACGAAGCCAGt                            >  2:413037/1‑141 (MQ=255)
                                                                                                                aCCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACca                        <  2:891562/141‑1 (MQ=255)
                                                                                                                 ccACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcac                       >  1:282337/1‑140 (MQ=255)
                                                                                                                  cacgGTATCCCCTGGTTTCGCTAGATACCCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGCTGTCTAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcatc                      <  2:323123/141‑1 (MQ=255)
                                                                                                                  cacaGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACgatc                      >  2:827231/1‑141 (MQ=255)
                                                                                                                   acaGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcatca                     <  1:673636/141‑1 (MQ=255)
                                                                                                                    cggTATCACCTGGTTTACCTAGATACGCCCGAATCAAAACCTTACCTTCCGGTGTAGTAAGTGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcatcat                    <  2:941630/139‑1 (MQ=255)
                                                                                                                     aGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATCACGAAGCCAGTACCATCATg                   <  1:670559/141‑1 (MQ=255)
                                                                                                                     aGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGGAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCACGTACCATCATg                   >  1:380317/1‑141 (MQ=255)
                                                                                                                      gTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATgg                  >  1:1034905/1‑141 (MQ=255)
                                                                                                                      gTATCACCTGGTTTACCTAGATACGCACGAACTAAAACCTTACCTTCCGGTGTAGTAAATGGCCATGTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATgg                  <  2:574586/141‑1 (MQ=255)
                                                                                                                      gTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATgg                  >  1:896084/1‑141 (MQ=255)
                                                                                                                       tATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCTATTTTCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAGGCCAGTACcatccttgg                 >  2:506271/1‑136 (MQ=255)
                                                                                                                        aTCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGt                <  1:726277/141‑1 (MQ=255)
                                                                                                                                  ttACCTAGATACGCACGAATCAAACCCTTACCTTCCGGTGTAGTAAATGGCCATTTTTGCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATg      <  2:300996/141‑1 (MQ=255)
                                                                                                                                  ttACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCTTTTTTTTTGTACATGCGGTAATGTCCGTATCACTCGTTCTCGCAATTGCGAAGCCACTACCGCCCTGGGTATTTTTAAtt      >  2:899690/1‑140 (MQ=255)
                                                                                                                                  ttACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATg      >  1:492035/1‑141 (MQ=255)
                                                                                                                                  ttACCTAGACACGCACGAAACAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCCGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATg      <  1:414323/141‑1 (MQ=255)
                                                                                                                                  ttACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTCCCATCATGGGTAGTTTCAAtt      >  1:11309/1‑140 (MQ=255)
                                                                                                                                   tACCTAGATACGCACGGATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGt     >  1:525861/1‑141 (MQ=255)
                                                                                                                                    aCCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAGGCCAGTACCATCATGGGTCTTTTCACTGTc    >  2:559075/1‑141 (MQ=255)
                                                                                                                                      cTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCtt  >  2:188701/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                         
ACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCTT  >  USA300TCH1516_ALE/1963631‑1963905

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: