Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,249,530 G→A 5.7% A956T (GCG→ACG)  smc_2 → Chromosome partition protein Smc

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,249,5300GA 5.7% 202.7 / 3.6 71A956T (GCG→ACG) smc_2Chromosome partition protein Smc
Reads supporting (aligned to +/- strand):  ref base G (35/31);  new base A (2/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.83e-01

CAAGTTTGTCACCAAGATATTTTAGCTATCGAAAATCACTACCAAGATATTAAAGCTGAACAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGTCCTGTAAACTTAAATGCAATTGAACAATTTGAAGAGTTAAATGAACGTTA  >  USA300TCH1516_ALE/1249395‑1249666
                                                                                                                                       |                                                                                                                                        
aaaGTTTGTCAGCAAGATATTTTAGCTAGCGAAAATCACTACCAAGATATTAAAGCTGAACAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGaaa                                                                                                                                     <  1:427070/140‑1 (MQ=255)
       gTCACCAAGATATTTTAGCTATCGAAAATCACTACCAAGATATTAAAGCTGACTAGTCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAat                                                                                                                              >  1:435411/1‑141 (MQ=255)
            cAAGATATTTTAGCTATCGAAAATCACTACCAAGATATTAAAGCTGAACAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATAca                                                                                                                         >  2:168365/1‑140 (MQ=255)
                atatTTTAGCTATCGAAAATCACTACCAAGATATTAAAGCTGAACAATCAAAGCTAGATGTATTAATTCATCATGCGATATATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGtg                                                                                                                     >  2:780983/1‑141 (MQ=255)
                   ttttAGCTATCGAAAATCACTACCAAGATATTAAAGCTGAACAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAGCGTGCGGAATCTGCATATACAATtgata                                                                                                                  >  2:480006/1‑140 (MQ=255)
                   ttttAGCTATCGAAAATCACTACCAAGATATTAAAGCTGAACAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGtgatg                                                                                                                  >  1:292605/1‑141 (MQ=255)
                   ttttAGCTATCGAAAATCACTACCAAGATATTAAAGCTGAACAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGtgatg                                                                                                                  >  1:322231/1‑141 (MQ=255)
                   ttttAGCTATCGAAAATCACTACCAAGATATTAAAGCTGAACAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGtgatg                                                                                                                  >  2:262873/1‑141 (MQ=255)
                     ttAGCTATCGAAAATCACTACCAAGATATTAAAGCTGAACAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGaa                                                                                                                >  2:310739/1‑141 (MQ=255)
                     ttAGCTATCGAAAATCACTACCAAGATATTAAAGCTGAACAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGaa                                                                                                                >  1:604158/1‑141 (MQ=255)
                       aGCTATCGAAAATCACTACCAAGATATTAAAGCTGAACAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATc                                                                                                              <  1:786378/141‑1 (MQ=255)
                        acTATCGAAAATCACTACCAAGATATTAAAGCTGAACAATCAAAGCTAGATGGATTAATTCATCATGCGATAGATCATTTAAATGATGAAGATCAACTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCg                                                                                                             <  1:512094/140‑1 (MQ=255)
                           aTCGAAAATCACTACCAAGATATTAAAGCTGAACAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGAtt                                                                                                          <  1:541186/141‑1 (MQ=255)
                            tCGAAAATCACTACCAAGATATTAAAGCTGAACAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTg                                                                                                         >  1:871927/1‑141 (MQ=255)
                            tCGAAAATCACTACCAAGATATTAAAGCTGAACAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTg                                                                                                         >  2:594949/1‑141 (MQ=255)
                                aaaTCACTACCAAGATATTAAAGCTGAACAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGCCGc                                                                                                     >  2:608081/1‑141 (MQ=255)
                                    caatACCAAAATATTAAAGCTGAACAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGACATCTGAATATACGAGTGATGAATCGATTGACGCATTa                                                                                                 <  2:686327/138‑1 (MQ=255)
                                    cACTACCAAGATAATAAAGCTGAACAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTa                                                                                                 <  1:473138/141‑1 (MQ=255)
                                     aCTACCAAGATATTAAAGCTGAACAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTAc                                                                                                >  2:488182/1‑141 (MQ=255)
                                     aCTACCAAGATATTAAAGCTGAACAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATAGATTGACGCATTaa                                                                                                >  2:372041/1‑140 (MQ=255)
                                       tACCAAGATATTAAAGCTGAACAATCAAAGCTAGATGTATTAATTCATAATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGAGGAAATACga                                                                                              >  2:539878/1‑140 (MQ=255)
                                         ccAAGATATTAAAGCTGAACAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTaa                                                                                            <  1:490114/141‑1 (MQ=255)
                                         ccAAGATATTAAAGCTGAACAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTACGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTaa                                                                                            >  2:158114/1‑141 (MQ=255)
                                         ccAAGATATTAAAGCTGAACAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTACGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTaa                                                                                            >  2:389827/1‑141 (MQ=255)
                                                    aaGCTGAACAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAGTGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATAGATTGATGCATTACGTGAAAacgttaagt                                                                                 >  1:895365/1‑141 (MQ=255)
                                                          atcAATCAAAGTTAGACGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATAGACGAGTGCTGGTACAAATGACGCATTACGTAAAAAAGTTCAGTTAATGa                                                                           <  1:667224/139‑1 (MQ=255)
                                                          aaCAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGa                                                                           <  2:316971/141‑1 (MQ=255)
                                                            cAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAg                                                                         <  1:488593/141‑1 (MQ=255)
                                                              aTCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGAt                                                                       >  2:297997/1‑141 (MQ=255)
                                                               tCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTACTGGAGATg                                                                      >  2:136559/1‑141 (MQ=255)
                                                               tCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTGAATGAAGATg                                                                      >  2:501864/1‑141 (MQ=255)
                                                               tCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGAGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGagg                                                                      >  1:173790/1‑139 (MQ=255)
                                                                  aaGCTATATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCg                                                                   <  2:427024/141‑1 (MQ=255)
                                                                   aGCTAGATGTATTAATACAACATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGa                                                                  <  1:4772/141‑1 (MQ=255)
                                                                    gCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGAt                                                                 <  2:615307/141‑1 (MQ=255)
                                                                     cTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGAtt                                                                >  1:226821/1‑141 (MQ=255)
                                                                        gATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCCATTGACTGTTGAACGTGCGAAATCTGACAATACGAGTGATGAATCGATTGACGCATTACGTAATAAAGTTAAGTTAATGTAAGTGTCGATtgct                                                             >  2:579439/1‑139 (MQ=255)
                                                                        gATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATtgct                                                             >  2:913930/1‑139 (MQ=255)
                                                                        gATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATtgat                                                             >  2:6470/1‑141 (MQ=255)
                                                                        gATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTAAACGTGCGAAATCTCAATATAAGCTTGTTGAACTGATTGACTCATTACGTAGAAAAGTTAAGTTAATCAAAATGTCGATtgat                                                             >  2:78749/1‑141 (MQ=255)
                                                                            tATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAAc                                                         <  2:604158/141‑1 (MQ=255)
                                                                             aTTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTACGTTAATGAAGATGTCGATTGATGAACt                                                        >  2:372458/1‑141 (MQ=255)
                                                                             aTTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACt                                                        <  1:636544/141‑1 (MQ=255)
                                                                             aTTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACt                                                        >  1:841898/1‑141 (MQ=255)
                                                                              ttAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTa                                                       >  1:13332/1‑141 (MQ=255)
                                                                                 aTTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGt                                                    <  1:103290/141‑1 (MQ=255)
                                                                                 aTTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGt                                                    >  1:653731/1‑141 (MQ=255)
                                                                                     atcatGCGATAGATCATTTAAATGCTGAGTATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGTCCTg                                                >  1:578201/1‑141 (MQ=255)
                                                                                     atcatGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGTCCTg                                                <  1:327928/141‑1 (MQ=255)
                                                                                         tGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGGGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGTCCTGTaaa                                            <  1:98721/141‑1 (MQ=255)
                                                                                             aTAGATCATTTAAATGATGAATATCAATTGACTGTTAAACGTGAGAAATCTGAATAAACGGGGGATGAATCGATTGAAGCATTAGGTAAAAACGTTAAGTTAATGAAGATGTCGATTGATGAAGTAGGTCCTGtccacgaa                                        >  2:135780/1‑134 (MQ=255)
                                                                                               aGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCGGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGTCCTGTAAACTTaaa                                      >  2:316957/1‑141 (MQ=255)
                                                                                                aaTCTTTTAAATTATGAATATCAATTCACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGTCCTGTAAACTTAAAt                                     <  1:892230/140‑1 (MQ=255)
                                                                                                aaTCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGTCCTGTAAACTTAAAt                                     <  1:170739/140‑1 (MQ=255)
                                                                                                        aaaTGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGTCCTGTAAACTTAAATGCAATTGa                             <  2:64062/141‑1 (MQ=255)
                                                                                                             atgaATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGTCCTGTAAACTTAAATGCAATTGAACAAt                        >  2:535289/1‑141 (MQ=255)
                                                                                                                 atatCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTGGATTGATTAACGAGGTCTTGTAAACTTAAATGCAATTGAACAATTTGa                    >  2:189355/1‑141 (MQ=255)
                                                                                                                 atatCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGTCCTGTAAACTTAAATGCAATTGAACAATTTGa                    <  2:614954/141‑1 (MQ=255)
                                                                                                                 atatCAATTGACTGTTGAACGCGCGAAATCTGAATATACGATTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGTCCTGTAAACTTAAATGCAATTTAACAATTTGa                    <  2:615065/141‑1 (MQ=255)
                                                                                                                  tatCAATTGACAGGTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGCTTGATGCGTTACGTAAAAAAGTTAAGTTCATGAAGATGTCGATTGATGAACTAGGTGTTTGAAACTTAGATGCAATTGAACAATTTGaa                   <  2:602508/141‑1 (MQ=255)
                                                                                                                      aaTTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGTCCTGTAAACCTAAATGCAATAGAACAATCTGAAGAGt               <  2:440240/141‑1 (MQ=255)
                                                                                                                          gaatggtGAACGTCAGAAAGCTGAATAAACGAGGGATGAATCGATTGAAAAAAAAAAAAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGTCCTGTAAACTTAAATGCAATTGAACAATTTGAAGAGTTaaa           <  2:58658/135‑1 (MQ=255)
                                                                                                                           aCTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGTCCTGTAAACTTAAATGCAATTGAACAATTTGAAGAGTTAAAt          <  1:561710/141‑1 (MQ=255)
                                                                                                                           aCTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGTCCTGTAAACTTAAATGCAATTGAACAATTTGAAGAGTTAAAt          <  1:142916/141‑1 (MQ=255)
                                                                                                                           aCTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGTCCTGTAAACTTAAATGCAATTGAACAATTTGAAGAGTTAAAt          <  1:825881/141‑1 (MQ=255)
                                                                                                                             tGTTGAACGAGAGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGTCCTGTAAACTTAAATGCAATTGAACAATTTGAAGAGTTAAATGa        <  2:723806/141‑1 (MQ=255)
                                                                                                                                  aaCGTACGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGTCCTGTAAACTTAAATGCAATTGAACAATTTGAAGAGTTAAATGAACGtt   <  1:158114/141‑1 (MQ=255)
                                                                                                                                  aaCGTACGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGTCCTGTAAACTTAAATGCAATTGAACAATTTGAAGAGTTAAATGAACGtt   <  1:389827/141‑1 (MQ=255)
                                                                                                                                   aCGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGTCCTGTAAACTTAAATGCAATTGAACAATTTGAAGAGTTAAATGACAGTta  >  2:133040/1‑141 (MQ=255)
                                                                                                                                   aCGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGTCCTGTAAACTTAAATGCAATTGAACAATTTGAAGAGTTAAATGAACGTta  <  2:134364/141‑1 (MQ=255)
                                                                                                                                   aCGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGTCCTGTAAACTTAAATGCAATTGAACAATTTGAAGAGTTAAATGAACGTta  <  2:248610/141‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                        
CAAGTTTGTCACCAAGATATTTTAGCTATCGAAAATCACTACCAAGATATTAAAGCTGAACAATCAAAGCTAGATGTATTAATTCATCATGCGATAGATCATTTAAATGATGAATATCAATTGACTGTTGAACGTGCGAAATCTGAATATACGAGTGATGAATCGATTGACGCATTACGTAAAAAAGTTAAGTTAATGAAGATGTCGATTGATGAACTAGGTCCTGTAAACTTAAATGCAATTGAACAATTTGAAGAGTTAAATGAACGTTA  >  USA300TCH1516_ALE/1249395‑1249666

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: