Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 2,030,815 G→T A90S (GCA→TCA)  USA300HOU_RS10250 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,030,8150GT98.0% 351.8 / 0.6 98A90S (GCA→TCA) USA300HOU_RS10250hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (48/48);  total (48/50)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

ACAATTGAACCAACTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTGCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGTCTCGAACATTGACAGCA  >  USA300TCH1516_ALE/2030680‑2030952
                                                                                                                                        |                                                                                                                                         
acaATTGAACCAACTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACg                                                                                                                                      <  4:96090/141‑1 (MQ=255)
 caATTGAACCAACTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGt                                                                                                                                     >  1:614259/1‑141 (MQ=255)
 caATTGAACCAACTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGt                                                                                                                                     >  4:234228/1‑141 (MQ=255)
   aTTGAACCAACTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTc                                                                                                                                   >  2:196771/1‑141 (MQ=255)
    ttGAACCAACTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCa                                                                                                                                  >  1:97898/1‑141 (MQ=255)
       aaCCAACTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGa                                                                                                                               <  1:116670/141‑1 (MQ=255)
         ccAACTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGAtt                                                                                                                             <  1:438726/141‑1 (MQ=255)
         ccAACTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGAtt                                                                                                                             >  1:485291/1‑141 (MQ=255)
         ccAACTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGAtt                                                                                                                             <  1:483673/141‑1 (MQ=255)
         ccAACTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGAtt                                                                                                                             >  2:551511/1‑141 (MQ=255)
            aCTGCCACGGCATCACTTATTCAAGGTTATCGTAGATTACCTTCTACGGTAATTTGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTCTTACATATGCATTGAGTTCAACGTTCACGATTTTa                                                                                                                          >  4:234594/1‑141 (MQ=255)
            aCTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTa                                                                                                                          >  4:234625/1‑141 (MQ=255)
            aCTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTa                                                                                                                          >  1:137392/1‑141 (MQ=255)
               gCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTc                                                                                                                       <  1:266970/141‑1 (MQ=255)
                     gCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACa                                                                                                                 >  1:513598/1‑141 (MQ=255)
                     gCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACa                                                                                                                 >  4:157691/1‑141 (MQ=255)
                       aTCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACaaa                                                                                                               <  2:631015/141‑1 (MQ=255)
                         cACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAAc                                                                                                             >  1:5/1‑141 (MQ=255)
                            ttattaAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTc                                                                                                           <  1:451568/141‑1 (MQ=255)
                               ttaAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAAt                                                                                                       <  3:270290/141‑1 (MQ=255)
                                taAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAAtt                                                                                                      <  2:338027/141‑1 (MQ=255)
                                 aaaGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGa                                                                                                                               <  2:613098/115‑1 (MQ=255)
                                 aaaGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGa                                                                                                                               >  1:613098/1‑115 (MQ=255)
                                 aaaGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTg                                                                                                     >  1:462976/1‑141 (MQ=255)
                                  aaGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGc                                                                                                    <  2:429991/141‑1 (MQ=255)
                                         tCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGTTTTTATTCTTTACAAAACTCAAATTGCAAGTTgg                                                                                             <  3:7656/141‑1 (MQ=255)
                                         tCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTgg                                                                                             <  4:204435/141‑1 (MQ=255)
                                          cGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGt                                                                                            >  4:225788/1‑141 (MQ=255)
                                          cGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGt                                                                                            >  1:77703/1‑141 (MQ=255)
                                          cGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGt                                                                                            >  1:387145/1‑141 (MQ=255)
                                           gTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTa                                                                                           <  4:156622/141‑1 (MQ=255)
                                             aGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACg                                                                                         <  1:445498/141‑1 (MQ=255)
                                              gATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGa                                                                                        <  2:177438/141‑1 (MQ=255)
                                              gATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGa                                                                                        >  1:125047/1‑141 (MQ=255)
                                               aTTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGaa                                                                                       >  4:81677/1‑141 (MQ=255)
                                                  aCCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATTTTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTAGTAcagcgac                                                                                    <  3:213017/141‑7 (MQ=255)
                                                     tGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGc                                                                                 >  1:449007/1‑141 (MQ=255)
                                                     tGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTGATTCTGCGATCAAGCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGc                                                                                 <  1:383833/141‑1 (MQ=255)
                                                      gCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCt                                                                                <  2:338057/141‑1 (MQ=255)
                                                      gCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCt                                                                                >  2:312930/1‑141 (MQ=255)
                                                          cGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTAc                                                                            <  1:278184/141‑1 (MQ=255)
                                                          cGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTAc                                                                            <  4:15833/141‑1 (MQ=255)
                                                           ggTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACt                                                                           >  3:162873/1‑141 (MQ=255)
                                                             tAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTaa                                                                         <  2:133876/141‑1 (MQ=255)
                                                             tAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTaa                                                                         >  3:10416/1‑141 (MQ=255)
                                                               aTTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAAca                                                                       >  4:183434/1‑141 (MQ=255)
                                                               aTTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAAca                                                                       >  2:241280/1‑141 (MQ=255)
                                                                ttGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAAcat                                                                      <  1:96232/141‑1 (MQ=255)
                                                                   gTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAAcatcat                                                                   <  1:287274/141‑1 (MQ=255)
                                                                     gCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACACAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTa                                                                 >  4:242093/1‑141 (MQ=255)
                                                                     gCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTa                                                                 >  4:241855/1‑141 (MQ=255)
                                                                          cTTTGCAGTATTATTTACATATTTATTCGGCGAACAATCACCATTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTAt                                                            <  1:125432/141‑1 (MQ=255)
                                                                             ttGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGa                                                          >  2:402826/1‑141 (MQ=255)
                                                                              tGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGAt                                                         >  3:227915/1‑141 (MQ=255)
                                                                                cAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTc                                                       <  2:355769/141‑1 (MQ=255)
                                                                                cAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTc                                                       <  1:141990/141‑1 (MQ=255)
                                                                                  gTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCa                                                     >  2:519088/1‑141 (MQ=255)
                                                                                       attTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAGCATCATTAGCTATGATTCCAGGTAt                                                <  1:59634/141‑1 (MQ=255)
                                                                                        tttACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTAtt                                               >  1:344464/1‑141 (MQ=255)
                                                                                           aCATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTCAAATCATTAGCTATGATTCCAGGTATTCat                                            <  2:485291/141‑1 (MQ=255)
                                                                                            cATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCatg                                           >  2:502817/1‑141 (MQ=255)
                                                                                              tatttattCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCatgat                                         >  3:94455/1‑141 (MQ=255)
                                                                                                 ttattCGGCGATCAATCACCTTTTACATATTCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGcc                                      <  3:75746/141‑1 (MQ=255)
                                                                                                 ttattCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTTTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGcc                                      >  1:579734/1‑141 (MQ=255)
                                                                                                 ttattCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGcc                                      <  2:467794/141‑1 (MQ=255)
                                                                                                 ttattCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGcc                                      >  2:307917/1‑141 (MQ=255)
                                                                                                 ttattCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGcc                                      >  1:48789/1‑141 (MQ=255)
                                                                                                     tCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACa                                  <  4:162873/141‑1 (MQ=255)
                                                                                                        gCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACatat                               <  4:3853/141‑1 (MQ=255)
                                                                                                        gCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACatat                               >  2:462245/1‑141 (MQ=255)
                                                                                                         cGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATAtt                              <  1:402826/141‑1 (MQ=255)
                                                                                                         cGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATAtt                              <  3:81677/141‑1 (MQ=255)
                                                                                                           atcaatcaCCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCa                            >  1:322180/1‑141 (MQ=255)
                                                                                                            tcaatcaCCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCaa                           >  2:281894/1‑141 (MQ=255)
                                                                                                            tcaatcaCCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCAGTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCaa                           >  4:246656/1‑141 (MQ=255)
                                                                                                             caatcaCCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAAc                          >  1:516926/1‑141 (MQ=255)
                                                                                                              aatcaCCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACt                         <  1:295314/141‑1 (MQ=255)
                                                                                                              aatcaCCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACt                         <  2:5/141‑1 (MQ=255)
                                                                                                               atcaCCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACtt                        >  3:226026/1‑141 (MQ=255)
                                                                                                                 caCCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGt                                                                                    >  4:236187/1‑79 (MQ=255)
                                                                                                                 caCCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGt                                                                                    <  3:236187/79‑1 (MQ=255)
                                                                                                                     ttttACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGtc                  <  3:229957/141‑1 (MQ=255)
                                                                                                                     ttttACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGtc                  <  1:241280/141‑1 (MQ=255)
                                                                                                                     ttttACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGtc                  <  4:41562/141‑1 (MQ=255)
                                                                                                                      tttACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGtct                 <  2:619413/141‑1 (MQ=255)
                                                                                                                         aCATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGTCTCGa              <  1:298260/141‑1 (MQ=255)
                                                                                                                         aCATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGTCTCGa              >  4:33537/1‑141 (MQ=255)
                                                                                                                         aCATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGTCTCGa              <  2:125047/141‑1 (MQ=255)
                                                                                                                         aCATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGTCTCGa              >  2:169359/1‑141 (MQ=255)
                                                                                                                          cATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGTCTCGaa             >  4:135920/1‑141 (MQ=255)
                                                                                                                          cATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCGGTACTAACATCATTAGCTATGATACCAGGTATTCATGATGCCTACATATTCAACTTTTTGTCTCGaa             >  2:383523/1‑141 (MQ=255)
                                                                                                                             atGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGTCTCGAACAt          <  1:635465/141‑1 (MQ=255)
                                                                                                                             atGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGTCTCGAACAt          <  3:135192/141‑1 (MQ=255)
                                                                                                                                 aTTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGTCTCGAACATTGAc      <  1:312930/141‑1 (MQ=255)
                                                                                                                                  ttGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGTCTCGAACATTGACa     <  3:22151/141‑1 (MQ=255)
                                                                                                                                    gAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGTCTCGAACATTGACAGc   >  3:11893/1‑141 (MQ=255)
                                                                                                                                    gAGTGCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACTTATTCAACTTTTTGTCTCGAACATTGACAGc   <  2:282661/141‑1 (MQ=255)
                                                                                                                                     aGTGCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGTCTCGAACATTGACAGCa  <  2:273075/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                         
ACAATTGAACCAACTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTGCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGTCTCGAACATTGACAGCA  >  USA300TCH1516_ALE/2030680‑2030952

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