Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,135,045 G→A 5.3% R55C (CGT→TGT)  USA300HOU_RS05620 ← Nicotinamide‑nucleotide adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,135,0450GA 5.3% 228.8 / 2.8 76R55C (CGT→TGT) USA300HOU_RS05620Nicotinamide‑nucleotide adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (33/39);  new base A (2/2);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.23e-01

ACTACCAAATGCAATTGTATCGACACTCATATAATCTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCATTTGTAAGTTTTTTAGATTGATTAATATGATATTGATGCCCATTATGAAAGGGATTATAT  >  USA300TCH1516_ALE/1134905‑1135181
                                                                                                                                            |                                                                                                                                        
aCTACCAAATGCAATTGTATCGACACTCATATAATCTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTCATAATGCCATTTTTGCACg                                                                                                                                          <  1:577712/141‑1 (MQ=255)
 cTACCAAATGCAATTGTATCGACACTCATATAATCTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGa                                                                                                                                         >  1:147561/1‑141 (MQ=255)
    ccAAATGCAATTGTATCGACACTCATATAATCTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGGTGATAATGCAATTTTTGCACGGGTa                                                                                                                                      >  2:12368/1‑141 (MQ=255)
      aaaTGCAATTGTATCGACACTCATATAATCTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATCACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTaaa                                                                                                                                    >  2:54034/1‑141 (MQ=255)
      aaaTGCAATTGTATCGACACTCATATAATCTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTaaa                                                                                                                                    >  2:462120/1‑141 (MQ=255)
       aaTGCAATTGTATCGACACTCATATAATCTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAAc                                                                                                                                   >  2:202257/1‑141 (MQ=255)
       aaTGCAATTGTATCGACACTCATATAATCTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAAc                                                                                                                                   >  2:655937/1‑141 (MQ=255)
       aaTGCAATTGTATCGACACTCATATAATCTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACAAGTAAAc                                                                                                                                   >  2:143334/1‑141 (MQ=255)
       aaTGCAATTGTATCGACACTCATATAATCTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACAAGTAAAc                                                                                                                                   >  2:720877/1‑141 (MQ=255)
            aaTTGTATCGACACTCATATAATCTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACttatt                                                                                                                              <  2:289342/141‑1 (MQ=255)
                 tATCGACACTCATATACTCTGCGACTTTAACTGCCCGTTCGGCAAGTTGTTCACCAGCTTATAAACATTCAGTTACCGCTTCTTCGACAGCTAAATCTTTTCGTGATGATGCCATTTTTGCACGAGTAAACTTATTATAGa                                                                                                                         <  2:181488/141‑1 (MQ=255)
                 tATCGACACTCATATAATCTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGa                                                                                                                         >  1:88711/1‑141 (MQ=255)
                 gattgACACTCATATAAGCTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGa                                                                                                                         <  2:101518/137‑1 (MQ=255)
                       cacTCATATAATCTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTg                                                                                                                   <  1:655937/141‑1 (MQ=255)
                         cTCCTATAATCTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGt                                                                                                                 <  1:247831/141‑1 (MQ=255)
                           cATATAATCTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTc                                                                                                               <  1:158004/141‑1 (MQ=255)
                           cATATAATCTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAAATGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTGGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTc                                                                                                               >  2:258564/1‑141 (MQ=255)
                            atatAATCTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCg                                                                                                              >  2:496590/1‑141 (MQ=255)
                             tctaATCTGCGACTTTAACTGCGAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGCTAATGCCATTTTTGCACGAGTAACCTTATTATAGATTGCTGGTTCGc                                                                                                             <  2:623722/139‑1 (MQ=255)
                             tataATCTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGc                                                                                                             <  1:128137/141‑1 (MQ=255)
                              ataATCTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGcc                                                                                                            <  1:517899/141‑1 (MQ=255)
                                   cTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCa                                                                                                       >  1:93177/1‑141 (MQ=255)
                                   cTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCa                                                                                                       >  1:628412/1‑141 (MQ=255)
                                   cTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCa                                                                                                       >  2:723594/1‑141 (MQ=255)
                                   cTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCa                                                                                                       >  1:12552/1‑141 (MQ=255)
                                       caCTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGAc                                                                                                   <  1:12368/140‑1 (MQ=255)
                                         cTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATCATGCTTTTTTTCGACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACaa                                                                                                 <  2:278763/141‑1 (MQ=255)
                                           ttAACTGCTAGTTCGGCAAAATGGTCTCCAGATGATAACCTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAg                                                                                               <  2:122827/141‑1 (MQ=255)
                                             aaCTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAATTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGtt                                                                                             <  1:497978/141‑1 (MQ=255)
                                             aaCTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGtt                                                                                             >  1:715168/1‑141 (MQ=255)
                                              aCTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTTCTCGTTTTTTCTTATTATAGATTGCTGGTTAGCCACGCATGACAAAGTTa                                                                                            >  2:689608/1‑141 (MQ=255)
                                              aCTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTa                                                                                            >  2:689626/1‑141 (MQ=255)
                                                  ttggTTCCGCAATATGATCGCCAGATGATAAACTGGCAGTTGCTGGTTGTTCGATAACTATATCATCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTCGAGATTGCTGGTTCGCCACGCATGACAAAGTTACCAc                                                                                        <  2:47811/138‑1 (MQ=255)
                                                   tAGTTCGGCAAAATGATCGCCAGATGATAAATTGGCAGTTGCTGGTAGTGAGATAACTAAAACAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCAcc                                                                                       >  2:162595/1‑140 (MQ=255)
                                                   tAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCAcc                                                                                       >  2:8992/1‑140 (MQ=255)
                                                   tAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACt                                                                                       <  1:403209/141‑1 (MQ=255)
                                                    aGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATATTTGACAGAGTTAACTTATGATAGATTGCTTGTTCGCCACGCATGACAAAGTTACCACTc                                                                                      <  2:546208/141‑1 (MQ=255)
                                                       tCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCatt                                                                                   <  1:699551/141‑1 (MQ=255)
                                                       tCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCatt                                                                                   <  2:118148/141‑1 (MQ=255)
                                                        ctgCAAAATGATCCCCAGATGAGAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAAATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCatta                                                                                  <  2:692270/139‑1 (MQ=255)
                                                         ggCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAACCTTATTTTAGATTGCTGGTTGGCACCGCATGACAAAGTTACAACTCattat                                                                                 >  1:595113/1‑141 (MQ=255)
                                                            aaaaTGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGc                                                                              <  1:160446/141‑1 (MQ=255)
                                                              aaTGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGATGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCaa                                                                            >  2:766004/1‑141 (MQ=255)
                                                                tGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAGAGTTACCACGCATAATTGCAATa                                                                          >  1:448575/1‑141 (MQ=255)
                                                                  aTCGACAGATGAGAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCCGTTGATAATGCCATTTTTGCACGAGTAAACATATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGt                                                                        <  1:762353/141‑1 (MQ=255)
                                                                    cGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAAc                                                                     <  2:675826/141‑1 (MQ=255)
                                                                     gCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAAc                                                                     >  2:329307/1‑141 (MQ=255)
                                                                             atAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACAAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCAt                                                             <  1:143334/141‑1 (MQ=255)
                                                                             atAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACAAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCAt                                                             <  1:720877/141‑1 (MQ=255)
                                                                              tAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAACCTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCAtt                                                            <  2:620059/141‑1 (MQ=255)
                                                                              tAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCAtt                                                            <  2:521974/141‑1 (MQ=255)
                                                                                  ttGGCAGTTGCGGGTGGTTGGATAACTAAATCAGCTGTTGATAATGCCATTGTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACGCATTATTGCACTAGTAACGTCAGCATTTGTa                                                        <  2:335826/140‑1 (MQ=255)
                                                                                  cTGGCAGTTGCTGTTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCATTTGTa                                                        <  2:624074/141‑1 (MQ=255)
                                                                                  cTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCCCGCATGACAATGTCACCACTCATTATTGCAATAGTAACGTCAGCATTTGTa                                                        <  1:408675/141‑1 (MQ=255)
                                                                                   tGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCATTTGTaa                                                       <  1:258564/141‑1 (MQ=255)
                                                                                         ttGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCATTTGTAAGttttt                                                 <  1:272221/141‑1 (MQ=255)
                                                                                         ttGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCATTGGTAAGttttt                                                 >  2:734019/1‑141 (MQ=255)
                                                                                          ggcgggTAGGTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCATTTGTAAGtttttt                                                <  2:379571/137‑1 (MQ=255)
                                                                                              ggTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCATTTGTAAGTTTTTTAgat                                            >  1:344388/1‑141 (MQ=255)
                                                                                               gTAGTTCGATAAGTAAATCAGCTGTTGATAATTGCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCACGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCCTTTGTAAGTTTTTTAgatt                                           <  2:433140/141‑1 (MQ=255)
                                                                                                   ttCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCATTTGTAAGTTTTTTAgattgatt                                       <  2:173043/141‑1 (MQ=255)
                                                                                                     cGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCATTTGTAAGTTTTTTAGATTGGTTaa                                     >  1:648950/1‑141 (MQ=255)
                                                                                                         accTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCATTTGTAAGTTTTTTAGATTGATTAATATg                                 <  2:637496/139‑1 (MQ=255)
                                                                                                            aaaaGCAGCGTTTGATAATGTCATTTTTGTACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCATTTGTAAGTTTTTTAGATTGATTAATATGata                              <  2:215113/140‑1 (MQ=255)
                                                                                                               aTCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCATTTGTAAGTTTTTTAGATTGATTAATATGATATTg                           >  2:56823/1‑141 (MQ=255)
                                                                                                                 cAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCATTTGTAAGTTTTTTAGATTGATTAATATGATATTGAt                         <  2:528930/141‑1 (MQ=255)
                                                                                                                    cTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCATTTGTAAGTTTTTTAGATTGATTAATATGATATTGATGcc                      <  1:428592/141‑1 (MQ=255)
                                                                                                                      gTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCATTTGTAAGTTTTTTAGATTGATTAATATGATATTGATGCCCa                    >  2:678726/1‑141 (MQ=255)
                                                                                                                           tAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCATTTGTAAGTTTTTTAGATTGATTAATATGATATTGATGCCCATTATg               <  2:250872/141‑1 (MQ=255)
                                                                                                                                   tttttGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCATTTGTAAGTTTTTTAGATTGATTAATATGATATTGATGCCCATTATGAAAGGGAt       <  1:766004/141‑1 (MQ=255)
                                                                                                                                    ttttGCACGCGTACACTTTTCACGGCTTGTTTGTTCGCCCCGCATGTCTCCGCCACGACCGATTATTGCGATTTTAACGTCAACATTTGTAAGTTTTGTAGATTGATTAATATGATATTGATGCCCATTATGAAAGGGAtt      <  2:448575/141‑1 (MQ=255)
                                                                                                                                     tttGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCATTTGTAAGTTTTTTAGATTGATTAATATGATATTGATGCCCATTATGCACGGGAGta     >  2:726563/1‑141 (MQ=255)
                                                                                                                                     tttGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCATTTGTAAGTTTTTTAGATTGATTAATATGATATTGATGCCCATTATGAAAGGGATta     >  1:586582/1‑141 (MQ=255)
                                                                                                                                      ttGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCATTTGTAAGTTTTTTAGATTGATTAATATGATATTGATGCCCATTATGAAAGGGATtat    >  1:360558/1‑141 (MQ=255)
                                                                                                                                        gCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCATTTGTAAGTTTTTTGGAGGGATTAATATGATATTGATGCCCATTATGAAAGGGATTatat  >  1:294654/1‑141 (MQ=255)
                                                                                                                                        gCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCATTTGTAAGTTTTTTAGATTGATTAATATGATATTGATTGCCATTATGAAGAGGATTatat  >  2:434845/1‑141 (MQ=255)
                                                                                                                                        gCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCATTTGTAAGTTTTTTAGATTGATTAATATGATATTGATGCCCATTATGAAAGGGATTatat  >  2:235075/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                        
ACTACCAAATGCAATTGTATCGACACTCATATAATCTGCGACTTTAACTGCTAGTTCGGCAAAATGATCGCCAGATGATAAACTGGCAGTTGCTGGTAGTTCGATAACTAAATCAGCTGTTGATAATGCCATTTTTGCACGAGTAAACTTATTATAGATTGCTGGTTCGCCACGCATGACAAAGTTACCACTCATTATTGCAATAGTAACGTCAGCATTTGTAAGTTTTTTAGATTGATTAATATGATATTGATGCCCATTATGAAAGGGATTATAT  >  USA300TCH1516_ALE/1134905‑1135181

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