Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,424,608 A→G 6.2% intergenic (‑166/+62) 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,6080AG 6.2% 728.4 / 23.2 242intergenic (‑166/+62)USA300HOU_RS12490/USA300HOU_RS12495hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (84/137);  new base G (7/8);  total (97/145)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 246 total aligned reads displayed.

CCTGTCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAGCA  >  USA300TCH1516_ALE/2424518‑2424698
                                                                                          |                                                                                           
 cTGTCGATGATTGATTCACAATAGCACGTGATTTTTTTTTTTTTCTATATTTTTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTAc                                                                                            >  2:880078/1‑91 (MQ=255)
  tGTCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACa                                                                                           <  2:605825/91‑1 (MQ=255)
  tGTCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACa                                                                                           >  1:283711/1‑91 (MQ=255)
   gTCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaa                                                                                          <  1:3455683/91‑1 (MQ=255)
    tCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaa                                                                                         <  1:317285/91‑1 (MQ=255)
    tCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaa                                                                                         >  1:1912643/1‑91 (MQ=255)
    tCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaa                                                                                         >  1:2654619/1‑91 (MQ=255)
      gatgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                       <  1:3798695/91‑1 (MQ=255)
      gatgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                       <  2:1765843/91‑1 (MQ=255)
       atgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAAt                                                                                      <  1:1579969/91‑1 (MQ=255)
       atgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAAt                                                                                      >  1:1829841/1‑91 (MQ=255)
        tgattgatTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCTAAAACACCCTTTTTAATTTTTACAAAAAtt                                                                                     >  2:123943/1‑91 (MQ=255)
        tgattgatTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTCATTTTTACAGAAAtt                                                                                     >  2:3651516/1‑91 (MQ=255)
        tgattgatTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAAtt                                                                                     <  1:3625611/91‑1 (MQ=255)
        tgattgatTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAAtt                                                                                     >  2:2534053/1‑91 (MQ=255)
        tgattgatTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAAtt                                                                                     <  2:2159978/91‑1 (MQ=255)
        tgattgatTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAAtt                                                                                     >  2:1992892/1‑91 (MQ=255)
        tgattgagTCACAGTAGCACGTGATTTGTTTGTTTTTCAATATTATTCATTTAGTCCATCAAAGACACCCTTTTTAATTTTTAAAAAAAtt                                                                                     <  2:1529319/91‑1 (MQ=255)
         gattgattCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTa                                                                                   <  1:3899268/91‑1 (MQ=255)
          attgattCACAATAGCTCGTGATTTTTTTGTTTTTCAATATTATTCAGTTATTCCATCAAAAACACCCTTTTTAATTTTTGTAAAAATTaa                                                                                  <  2:174445/91‑1 (MQ=255)
          attgattCACAATAGCACGTGATTTTTTTGTTTTTCTTTTTTTTTCATTTATTCCATCAAAAACACTCTTTTTATTTTTTACAAAAATTaa                                                                                  >  2:3554028/1‑91 (MQ=255)
           ttgattCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTaaa                                                                                 <  1:1326208/91‑1 (MQ=255)
             gattCGCAATAGCACGTGATTTTTTTGTTTTTCAATATTCTTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTaaaaa                                                                               <  1:3672534/91‑1 (MQ=255)
               ttCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTGTAAAAATTaaaaagg                                                                             <  2:730235/91‑3 (MQ=38)
               ttCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTaaaaaaa                                                                             >  1:802072/1‑91 (MQ=255)
                tCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAg                                                                            >  1:409756/1‑91 (MQ=255)
                 cacaATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGt                                                                           >  1:2423659/1‑91 (MQ=255)
                 cacaATAGCACGTGATTTTATTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGt                                                                           <  2:1452256/91‑1 (MQ=255)
                  acaATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTCTTTACTTTTTACAAAAATTAAAACAAGTg                                                                          >  1:3814604/1‑91 (MQ=255)
                   caATAGCACGTGATTTTTTTTTTTTTTTTTATTATTCATTTATTCCATCAAAAACACCCTTTTTATTTTTTTCAAAAATTAAAAAAAGTGc                                                                         >  2:1366174/1‑91 (MQ=255)
                   caATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGc                                                                         <  1:2057545/91‑1 (MQ=255)
                    aaTAGCACGTGATTTTTTTGTTTTTCTATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCt                                                                        >  2:1718015/1‑91 (MQ=255)
                       aGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCt                                                                     <  1:83150/91‑1 (MQ=255)
                        gCACTTTATTTTTTTTCTTATCAATATTATTCAATTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTa                                                                    <  2:523898/91‑1 (MQ=255)
                        gCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTa                                                                    >  2:2988199/1‑91 (MQ=255)
                         cACGTGATTTTTTTGTTTTTCTATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTac                                                                   >  2:1143780/1‑91 (MQ=255)
                         cACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTAAAAAAAAAAAAAAAAAGTGCTCCTac                                                                   <  2:389119/91‑1 (MQ=255)
                           cGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTacac                                                                 >  2:2791403/1‑91 (MQ=255)
                           cGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTacac                                                                 >  1:1544422/1‑91 (MQ=255)
                            gTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACt                                                                <  2:3154623/91‑1 (MQ=255)
                             tGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTg                                                               <  2:3213788/91‑1 (MQ=255)
                             tGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTg                                                               <  1:2998238/91‑1 (MQ=255)
                              gATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTGACAAAAATTAAAAAAAGTGCTCCTACACTGc                                                              <  1:3738518/91‑1 (MQ=255)
                                  tttttGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGc                                                          >  1:880341/1‑91 (MQ=255)
                                    tttGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCAt                                                        >  1:2832693/1‑91 (MQ=255)
                                     ttGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATg                                                       >  1:3481787/1‑91 (MQ=255)
                                      tGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGt                                                      <  2:1512691/91‑1 (MQ=255)
                                       gTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTa                                                     <  2:3452775/91‑1 (MQ=255)
                                         ttttCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGa                                                   <  1:23355/91‑1 (MQ=255)
                                         ttttCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGa                                                   <  1:1360238/91‑1 (MQ=255)
                                          tttCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGaa                                                  <  1:2912233/91‑1 (MQ=255)
                                           ttCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGaaa                                                 >  2:3351776/1‑91 (MQ=255)
                                            tCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAac                                                <  1:1048676/91‑1 (MQ=255)
                                              aaTATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAacac                                              <  2:1819603/91‑1 (MQ=255)
                                              aaTATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAacac                                              <  2:2394023/91‑1 (MQ=255)
                                              aaTATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAacac                                              >  1:3126458/1‑91 (MQ=255)
                                                 attattCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACttt                                           <  2:15615/91‑1 (MQ=255)
                                                   tattCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACttttt                                         <  1:1253617/91‑1 (MQ=255)
                                                     ttCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCa                                       <  1:488276/91‑1 (MQ=255)
                                                     ttCATTTATTCCATCAAAAACACCCTTGTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCa                                       <  2:3714811/91‑1 (MQ=255)
                                                      tCATTTATTCCATATAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCATCACTGCTTGCATGTAGAAACACTTTTTCAt                                      <  1:3566114/91‑1 (MQ=255)
                                                        aTTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTg                                    >  1:2158504/1‑91 (MQ=255)
                                                          ttATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTa                                  <  1:911015/91‑1 (MQ=255)
                                                           tATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTaa                                 >  2:2763598/1‑91 (MQ=255)
                                                             ttCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATg                               >  2:3849067/1‑91 (MQ=255)
                                                              tCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGt                              <  2:701127/91‑1 (MQ=255)
                                                              tCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGt                              <  1:3581457/91‑1 (MQ=255)
                                                              tCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGt                              <  1:658066/91‑1 (MQ=255)
                                                               ccATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGtt                             <  1:1616647/91‑1 (MQ=255)
                                                                cATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTa                            >  2:1870809/1‑91 (MQ=255)
                                                                cATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTa                            >  1:2011073/1‑91 (MQ=255)
                                                                 aTCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAt                           <  2:41854/91‑1 (MQ=255)
                                                                 aTCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAt                           <  2:2424512/91‑1 (MQ=255)
                                                                 aTCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAt                           <  1:1654003/91‑1 (MQ=255)
                                                                 aTCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAt                           <  1:2591445/91‑1 (MQ=255)
                                                                   cAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttc                         <  2:1632603/91‑1 (MQ=255)
                                                                   cAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttc                         >  1:2310187/1‑91 (MQ=255)
                                                                      aaaCACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttcttc                      <  2:1995384/91‑1 (MQ=255)
                                                                      aaaCACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttcttc                      >  1:2107779/1‑91 (MQ=255)
                                                                       aaCACCCTTTTTAATTTTTACAAAAATTAAAATAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTCTAATGTTATtcttct                     <  1:3366666/91‑1 (MQ=255)
                                                                          acCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCga                  >  1:324148/1‑91 (MQ=255)
                                                                           ctttttttAATTTTTGTAAAAATTAAAAAGGGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCgag                 >  1:2966162/4‑91 (MQ=25)
                                                                              tttttAATTTTTGTAAAAATTAAAAAGGGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACa              >  1:3923106/1‑91 (MQ=37)
                                                                              tttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACa              <  2:3103851/91‑1 (MQ=255)
                                                                              tttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACa              <  2:2610925/91‑1 (MQ=255)
                                                                              tttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACa              >  1:3759096/1‑91 (MQ=255)
                                                                               ttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACAt             <  2:605477/91‑1 (MQ=255)
                                                                               ttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACAt             <  1:2439981/91‑1 (MQ=255)
                                                                                tttAATTTTTACAAAAATTAAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACAt             >  2:3398061/1‑91 (MQ=255)
                                                                                 ttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATAc           <  1:725986/91‑1 (MQ=255)
                                                                                  tAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACCGCTTGCATGTAGCAACACTTTTACATTGTAATGTTATTCTTCTCGAGACCTAcc          >  1:785541/1‑91 (MQ=255)
                                                                                   ttttttttTAAATATTAAAAAAATTGCTCCTACACTGCTTTCATGTTGATATACTTTTTCTTTGTACTGTTATTCTTCTCGTGACATACCt         >  2:2916976/3‑91 (MQ=255)
                                                                                   aaTTTTTGTAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTTTTGAGACATACCt         >  2:3312298/1‑91 (MQ=255)
                                                                                   aaTTTTTGTAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCt         >  1:1727493/1‑91 (MQ=255)
                                                                                   aaTTTTTACAAACATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCt         >  2:2871406/1‑91 (MQ=255)
                                                                                   aaTTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCt         >  2:1966745/1‑91 (MQ=255)
                                                                                   aaTTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCt         <  1:365025/91‑1 (MQ=255)
                                                                                    aTTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCtt        <  2:830405/91‑1 (MQ=255)
                                                                                      ttttACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCtttt      >  2:3014933/1‑91 (MQ=255)
                                                                                         tACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAGc   >  1:1929527/1‑91 (MQ=255)
                                                                                          |                                                                                           
CCTGTCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAGCA  >  USA300TCH1516_ALE/2424518‑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: