Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,424,609 C→T 7.9% intergenic (‑167/+61) USA300HOU_RS12490 ← / ← USA300HOU_RS12495 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,424,6090CT 7.9% 721.4 / 28.3 241intergenic (‑167/+61)USA300HOU_RS12490/USA300HOU_RS12495hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (84/135);  new base T (12/8);  total (96/145)
Fisher's exact test for biased strand distribution p-value = 9.30e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 244 total aligned reads displayed.

CTGTCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAGCA  >  USA300TCH1516_ALE/2424519‑2424698
                                                                                          |                                                                                          
cTGTCGATGATTGATTCACAATAGCACGTGATTTTTTTTTTTTTCTATATTTTTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTAc                                                                                            >  2:880078/1‑91 (MQ=255)
 tGTCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACa                                                                                           >  1:283711/1‑91 (MQ=255)
 tGTCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACa                                                                                           <  2:605825/91‑1 (MQ=255)
   tCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaa                                                                                         <  1:317285/91‑1 (MQ=255)
   tCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaa                                                                                         >  1:2654619/1‑91 (MQ=255)
   tCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaa                                                                                         >  1:1912643/1‑91 (MQ=255)
    cGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTAAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaa                                                                                        >  2:3254742/1‑91 (MQ=255)
     gatgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                       <  2:1765843/91‑1 (MQ=255)
     gatgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                       <  2:1554926/91‑1 (MQ=255)
     gatgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                       <  1:2418033/91‑1 (MQ=255)
      atgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAAt                                                                                      <  1:1579969/91‑1 (MQ=255)
      atgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAAt                                                                                      <  1:896247/91‑1 (MQ=255)
       tgattgatTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTCATTTTTACAGAAAtt                                                                                     >  2:3651516/1‑91 (MQ=255)
       tgattgatTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTGTAAAAAtt                                                                                     <  1:48079/91‑1 (MQ=255)
       tgattgatTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAAtt                                                                                     >  2:2534053/1‑91 (MQ=255)
       tgattgatTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAAtt                                                                                     >  2:1992892/1‑91 (MQ=255)
       tgattgagTCACAGTAGCACGTGATTTGTTTGTTTTTCAATATTATTCATTTAGTCCATCAAAGACACCCTTTTTAATTTTTAAAAAAAtt                                                                                     <  2:1529319/91‑1 (MQ=255)
        gattgattCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTa                                                                                   >  1:3450879/1‑91 (MQ=255)
         attgattCACAATAGCTCGTGATTTTTTTGTTTTTCAATATTATTCAGTTATTCCATCAAAAACACCCTTTTTAATTTTTGTAAAAATTaa                                                                                  <  2:174445/91‑1 (MQ=255)
         attgattCACAATAGCACGTGATTTTTTTGTTTTTCTTTTTTTTTCATTTATTCCATCAAAAACACTCTTTTTATTTTTTACAAAAATTaa                                                                                  >  2:3554028/1‑91 (MQ=255)
          ttgattCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTaaa                                                                                 <  1:1326208/91‑1 (MQ=255)
           tgattCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTaaaa                                                                                <  1:3641045/91‑1 (MQ=255)
              ttCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTGTAAAAATTaaaaagg                                                                             <  2:730235/91‑3 (MQ=38)
              ttCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTaaaaaaa                                                                             >  1:802018/1‑91 (MQ=255)
                 acaATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTCTTTACTTTTTACAAAAATTAAAACAAGTg                                                                          >  1:3814604/1‑91 (MQ=255)
                  caATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGc                                                                         <  1:2057545/91‑1 (MQ=255)
                     tAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTcc                                                                      <  1:805720/91‑1 (MQ=255)
                       gCACTTTATTTTTTTTCTTATCAATATTATTCAATTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTa                                                                    <  2:523898/91‑1 (MQ=255)
                          cGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTTAAAAAAGTGCTCCTacac                                                                 >  2:1124023/1‑91 (MQ=255)
                          cGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTacac                                                                 >  1:1544422/1‑91 (MQ=255)
                           gTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACt                                                                <  2:3154623/91‑1 (MQ=255)
                            tGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTg                                                               <  1:2998238/91‑1 (MQ=255)
                            tGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTg                                                               <  1:588365/91‑1 (MQ=255)
                             taTTTTTTTGTTTTTCAATATTATTCATTTATTCTATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGc                                                              <  2:1292967/90‑1 (MQ=255)
                             taTTTGTTTGTTTTACAAATTTATAAAGACATTCGAACATAAACCCCTATGTTCATTTTTACAGAAACGAAAAAAAGTGCTCCTACACTGc                                                              <  2:1509392/90‑1 (MQ=255)
                               tttttttGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCtt                                                            <  1:3567962/91‑1 (MQ=255)
                                ttttttGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTg                                                           >  1:3798602/1‑91 (MQ=255)
                                    ttGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATg                                                       >  1:3481787/1‑91 (MQ=255)
                                      gTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTa                                                     <  2:3452775/91‑1 (MQ=255)
                                        ttttCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGa                                                   <  1:1360238/91‑1 (MQ=255)
                                        ttttCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGa                                                   <  1:23355/91‑1 (MQ=255)
                                        ttttCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGa                                                   <  1:3083816/91‑1 (MQ=255)
                                         tttCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGaa                                                  <  1:2911445/91‑1 (MQ=255)
                                          ttCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGaaa                                                 >  2:3351776/1‑91 (MQ=255)
                                          ttCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGaaa                                                 <  1:875633/91‑1 (MQ=255)
                                           tCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAac                                                <  1:1048676/91‑1 (MQ=255)
                                             aaTATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAacac                                              <  2:1819603/91‑1 (MQ=255)
                                             aaTATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAacac                                              <  2:2394023/91‑1 (MQ=255)
                                                attattCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACttt                                           <  1:727327/91‑1 (MQ=255)
                                                attattCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACttt                                           >  1:367196/1‑91 (MQ=255)
                                                  tattCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACttttt                                         <  1:1253617/91‑1 (MQ=255)
                                                   attCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTc                                        >  1:701698/1‑91 (MQ=255)
                                                       aTTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTg                                    >  1:2158504/1‑91 (MQ=255)
                                                          tATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTaa                                 >  2:2763598/1‑91 (MQ=255)
                                                           aTTCCATCAAAAACACCCTTTTCAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGACACACTTTTTCATTGTAAt                                <  2:1175069/91‑1 (MQ=255)
                                                            ttCCATCAAAAACGCCCTTTTTAATTTTTACAAAAATTAAAAAGAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATg                               <  1:3958734/91‑1 (MQ=255)
                                                            ttCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATg                               >  2:3849067/1‑91 (MQ=255)
                                                             tCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGt                              <  2:701127/91‑1 (MQ=255)
                                                             tCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGt                              <  2:2998557/91‑1 (MQ=255)
                                                              ccATCAAAAACACCCTTTTTAATTTTTACAAAAGTTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGtt                             <  1:3240183/91‑1 (MQ=255)
                                                              ccATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGtt                             <  1:1616647/91‑1 (MQ=255)
                                                               cATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTa                            <  2:2895505/91‑1 (MQ=255)
                                                               cATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTa                            >  1:2011073/1‑91 (MQ=255)
                                                               cATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTa                            >  2:1870809/1‑91 (MQ=255)
                                                                aTCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAt                           <  2:41854/91‑1 (MQ=255)
                                                                aTCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAt                           <  2:3750830/91‑1 (MQ=255)
                                                                aTCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAt                           <  1:2591444/91‑1 (MQ=255)
                                                                aTCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAt                           <  1:1929233/91‑1 (MQ=255)
                                                                aTCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAt                           <  1:1654003/91‑1 (MQ=255)
                                                                  cAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttc                         <  2:1632603/91‑1 (MQ=255)
                                                                  cAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttc                         >  1:2310187/1‑91 (MQ=255)
                                                                  cAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttc                         <  2:3892335/91‑1 (MQ=255)
                                                                   aaaaaCACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttct                        <  1:3874368/91‑1 (MQ=255)
                                                                     aaaCACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttcttc                      >  1:2107779/1‑91 (MQ=255)
                                                                     aaaCACCCTTTGTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttcttc                      <  2:1995370/91‑1 (MQ=255)
                                                                      aaCACCCTTTTTAATTTTTACAAAAATTAAAATAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTCTAATGTTATtcttct                     <  1:3366666/91‑1 (MQ=255)
                                                                          ctttttttAATTTTTGTAAAAATTAAAAAGGGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCgag                 >  1:2966162/4‑91 (MQ=25)
                                                                            cTTTTTTTTTTTTTTTAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCGTTGTAATGTTATTCTTCTATAGAc               >  2:1670582/1‑91 (MQ=255)
                                                                            cTTTTTAATTTTTACAAAAATTAAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCgaga                <  1:3701749/91‑1 (MQ=255)
                                                                            cTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGAc               >  1:3535372/1‑91 (MQ=255)
                                                                             tttttAATTTTTGTAAAAATTAAAAAGGGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACa              >  1:3923079/1‑91 (MQ=38)
                                                                             tttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACa              >  1:3759096/1‑91 (MQ=255)
                                                                             tttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATCGTAATGTTATTCTTCTCGAGACa              >  2:2332445/1‑91 (MQ=255)
                                                                             tttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACGGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACa              <  2:2610966/91‑1 (MQ=255)
                                                                              ttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACAt             <  2:605477/91‑1 (MQ=255)
                                                                              ttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACAt             <  1:2439981/91‑1 (MQ=255)
                                                                               tttAATTTTTACAAAAATTAAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACAt             >  2:3398061/1‑91 (MQ=255)
                                                                                  ttttttttTAAATATTAAAAAAATTGCTCCTACACTGCTTTCATGTTGATATACTTTTTCTTTGTACTGTTATTCTTCTCGTGACATACCt         >  2:2916976/3‑91 (MQ=255)
                                                                                  aaTTTTTGTAAAAATTAAAAAGGGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCt         >  2:1370037/1‑91 (MQ=38)
                                                                                  aaTTTTTGTAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTTTTGAGACATACCt         >  2:3312298/1‑91 (MQ=255)
                                                                                  aaTTTTTGTAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCt         >  1:1727493/1‑91 (MQ=255)
                                                                                  aaTTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCt         >  2:2866140/1‑91 (MQ=255)
                                                                                  aaTTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCt         >  2:1966745/1‑91 (MQ=255)
                                                                                   aTTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCtt        <  2:830405/91‑1 (MQ=255)
                                                                                    tttttACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCttt       >  1:1811596/1‑91 (MQ=255)
                                                                                     ttttACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCtttt      >  2:3014933/1‑91 (MQ=255)
                                                                                       ttACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAg    <  1:282604/91‑1 (MQ=255)
                                                                                        tACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAGc   >  1:3598868/1‑91 (MQ=255)
                                                                                        tACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAGc   <  2:3132142/91‑1 (MQ=255)
                                                                                         aCAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAGCa  >  2:2426935/1‑91 (MQ=255)
                                                                                          |                                                                                          
CTGTCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAGCA  >  USA300TCH1516_ALE/2424519‑2424698

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