Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 4,449 C→A 8.0% Q171K (CAA→AAA)  recF → DNA replication and repair protein RecF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE4,4490CA 8.0% 782.1 / 11.8 240Q171K (CAA→AAA) recFDNA replication and repair protein RecF
Reads supporting (aligned to +/- strand):  ref base C (116/99);  new base A (2/23);  total (118/122)
Fisher's exact test for biased strand distribution p-value = 7.81e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 245 total aligned reads displayed.

TAAAAGGCTCTCCTCAAATAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATTCATGCTGGTATCACAAATGA  >  USA300TCH1516_ALE/4309‑4588
                                                                                                                                            |                                                                                                                                            
    aGGCTCTCCTCAAATAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCaaaa                                                                                                                                         >  2:9057/1‑141 (MQ=255)
      gCTCTCCTCAAATAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCaaaaaa                                                                                                                                       <  1:919444/141‑1 (MQ=255)
         ctcctcAAATAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAgg                                                                                                                                    >  1:444971/1‑141 (MQ=255)
         ctcctcAAATAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAgg                                                                                                                                    <  3:209903/141‑1 (MQ=255)
         ctcctcAAATAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAgg                                                                                                                                    >  3:556444/1‑141 (MQ=255)
         ctcctcAAATAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAgg                                                                                                                                    >  1:1465761/1‑141 (MQ=255)
           cctcAAATAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGAc                                                                                                                                  >  1:1033372/1‑141 (MQ=255)
           cctcAAATAAGACGACGCTTTATAGAAATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTAAAAAAAAAAAAAAAAAAGGAc                                                                                                                                  <  2:226110/141‑1 (MQ=255)
            ctcAAATAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACt                                                                                                                                 <  1:955359/141‑1 (MQ=255)
             tcAAATAAGAGGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAACAGTTACAATTAGGCCAAAAAAAGGACtt                                                                                                                                <  2:1496308/141‑1 (MQ=255)
             tcAAATAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTAAAATAAAAAAAAAAAAAGGACtt                                                                                                                                <  2:425334/141‑1 (MQ=255)
                aaTAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAAc                                                                                                                             <  1:888514/141‑1 (MQ=255)
                 aTAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTAAAAAAAAAAAAAAAAAAGGACTTAACa                                                                                                                            <  2:1215691/141‑1 (MQ=255)
                      acgacgCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGtt                                                                                                                       <  2:774541/141‑1 (MQ=255)
                       cgacgCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCa                                                                                                                                            >  2:1150746/1‑119 (MQ=255)
                       cgacgCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCa                                                                                                                                            <  1:1150746/119‑1 (MQ=255)
                         acgCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGa                                                                                                                    >  3:181879/1‑141 (MQ=255)
                          cgCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGaa                                                                                                                   >  3:161844/1‑141 (MQ=255)
                          cgCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGaa                                                                                                                   <  1:631103/141‑1 (MQ=255)
                            cTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGt                                                                                                                 >  3:121898/1‑141 (MQ=255)
                             tttATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTa                                                                                                                <  1:835993/141‑1 (MQ=255)
                              ttATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTAt                                                                                                               <  3:352735/141‑1 (MQ=255)
                              ttATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGAATTAACAATGTTGGAAGTAt                                                                                                               >  1:1590989/1‑141 (MQ=255)
                                 taGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTaa                                                                                                            >  1:649052/1‑141 (MQ=255)
                                 taGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTaa                                                                                                            >  4:113193/1‑141 (MQ=255)
                                 taGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTaa                                                                                                            <  1:384585/141‑1 (MQ=255)
                                  aGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAATAATAATTACTTAAAGCATTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTaaa                                                                                                           <  4:325287/141‑1 (MQ=255)
                                   aatatGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAAAAAAAAAAAAGGACTTAACAATGTTGGGAGTATTAAAt                                                                                                          <  2:696528/140‑1 (MQ=255)
                                     tatGGAGTTGGGCCCAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATca                                                                                                        <  1:240347/141‑1 (MQ=255)
                                      atGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATcag                                                                                                       >  2:954487/1‑141 (MQ=255)
                                      atGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATcag                                                                                                       >  3:527941/1‑141 (MQ=255)
                                         gAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATcagcag                                                                                                    <  1:1364069/141‑1 (MQ=255)
                                          aGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGt                                                                                                   <  1:180876/141‑1 (MQ=255)
                                          aGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGt                                                                                                   >  4:163777/1‑141 (MQ=255)
                                          aGTTGGGCCAAATTTCGGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGt                                                                                                   <  1:1335941/141‑1 (MQ=255)
                                           gTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGtt                                                                                                  >  3:299623/1‑141 (MQ=255)
                                           gTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGtt                                                                                                  >  4:522685/1‑141 (MQ=255)
                                           gTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGtt                                                                                                  >  1:899107/1‑141 (MQ=255)
                                            ttGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTACAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGttt                                                                                                 <  1:1607070/141‑1 (MQ=255)
                                             tGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTg                                                                                                >  3:306804/1‑141 (MQ=255)
                                              gggCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGc                                                                                               <  1:517658/141‑1 (MQ=255)
                                                ggcAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTg                                                                                             <  2:1071334/139‑1 (MQ=255)
                                                gCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATAAAAAAAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTg                                                                                             <  4:269605/141‑1 (MQ=255)
                                                     aTTTCTGCTGTTTACTTAAATGAATTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCAAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAatat                                                                                        <  4:28962/141‑1 (MQ=255)
                                                        tCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCa                                                                                     >  3:499027/1‑141 (MQ=255)
                                                          tgctgTTTACTTAAATTATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAAt                                                                                   <  1:1256445/141‑1 (MQ=255)
                                                          tgctgTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAAt                                                                                   >  2:459870/1‑141 (MQ=255)
                                                          tgctgTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAAt                                                                                   >  1:1501842/1‑141 (MQ=255)
                                                             tgTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGaa                                                                                <  1:675844/141‑1 (MQ=255)
                                                                ttACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCaaaa                                                                                                                                         <  1:1559075/81‑1 (MQ=255)
                                                                     tAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTg                                                                        >  2:142733/1‑141 (MQ=255)
                                                                       aaTGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGAt                                                                      >  2:254277/1‑141 (MQ=255)
                                                                        aTGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAat                                                                                          >  2:38086/1‑120 (MQ=255)
                                                                         tGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATaa                                                                    <  3:33366/141‑1 (MQ=255)
                                                                            tttAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCAAAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAAc                                                                  <  4:476703/141‑1 (MQ=255)
                                                                             ttAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGt                                                                >  2:1413863/1‑141 (MQ=255)
                                                                             ttAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGt                                                                >  2:262879/1‑141 (MQ=255)
                                                                              tAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTg                                                               <  1:1189614/141‑1 (MQ=255)
                                                                              tAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTg                                                               <  1:459870/141‑1 (MQ=255)
                                                                               aGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCAAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGc                                                              <  4:25987/141‑1 (MQ=255)
                                                                                 gacAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTACCAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGcac                                                            <  1:741306/139‑1 (MQ=255)
                                                                                 cTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGCCTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGcac                                                            >  1:572090/1‑141 (MQ=255)
                                                                                   cAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACAt                                                          <  2:682581/141‑1 (MQ=255)
                                                                                   cAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACAt                                                          >  3:253891/1‑141 (MQ=255)
                                                                                         cAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTAtt                                                    <  2:1322738/141‑1 (MQ=255)
                                                                                            cGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATcagcaaa                                                                                                   <  1:794319/91‑3 (MQ=255)
                                                                                            cGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCaa                                                 >  2:1357631/1‑141 (MQ=255)
                                                                                             gTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAag                                                >  2:597959/1‑141 (MQ=255)
                                                                                             gTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAag                                                >  1:544740/1‑141 (MQ=255)
                                                                                                  ttAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCt                                            <  2:892565/140‑1 (MQ=255)
                                                                                                       gCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTc                                      >  1:154041/1‑141 (MQ=255)
                                                                                                         aaaaGAATAATTACTTCAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGt                                    <  4:628501/141‑1 (MQ=255)
                                                                                                         aaaaGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGt                                    <  4:114186/141‑1 (MQ=255)
                                                                                                           aaGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTCTTCAAGAGCTAGAGTCGTTa                                  >  1:1262406/1‑141 (MQ=255)
                                                                                                             gAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGc                                >  2:1550552/1‑141 (MQ=255)
                                                                                                              tgaaatTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCt                               <  3:401275/138‑1 (MQ=255)
                                                                                                              aataatTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCt                               <  4:173526/141‑1 (MQ=255)
                                                                                                              aacaatTCCTTAAAGAAGTTAAAAAAAAAAAAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCGAGAGCTAGAGTCGTTAGCt                               <  4:595785/138‑1 (MQ=255)
                                                                                                               ataatTACTTAAAGCAGTTACAATAAAAAAAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTa                              <  2:1590989/141‑1 (MQ=255)
                                                                                                                taatTACTTAAAGCAATTAAAAATAGGCAAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTaa                             <  2:1295688/141‑1 (MQ=255)
                                                                                                                 aatTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTATCAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTaaa                            <  1:1106624/141‑1 (MQ=255)
                                                                                                                  atTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAAc                           >  2:1285451/1‑141 (MQ=255)
                                                                                                                   ttACTTAAAGCAGTTACCATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAAcc                          >  2:1658431/1‑141 (MQ=255)
                                                                                                                   ttACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAAcc                          >  4:507228/1‑141 (MQ=255)
                                                                                                                       ttAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGAtt                      <  2:1020760/141‑1 (MQ=255)
                                                                                                                         aaaGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATTCa                    >  2:622945/1‑141 (MQ=255)
                                                                                                                           aGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATTCATg                  >  1:1614177/1‑141 (MQ=255)
                                                                                                                            gCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATTCATGc                 >  1:405151/1‑141 (MQ=255)
                                                                                                                            gCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAACTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATTCATGc                 <  1:1301017/141‑1 (MQ=255)
                                                                                                                              aGTTACAATTAGGCCAANAAAAGGACTTGACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAACCCGATTCATGCTg               >  1:1082634/1‑141 (MQ=255)
                                                                                                                              aGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATTCATGCTg               <  1:582022/141‑1 (MQ=255)
                                                                                                                               gTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATTCATGCTgg              >  3:650913/1‑141 (MQ=255)
                                                                                                                                ttACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATTCATGCTGGt             >  4:347944/1‑141 (MQ=255)
                                                                                                                                ttACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATTCATGCTGGt             >  1:424158/1‑141 (MQ=255)
                                                                                                                                ttACAATTAGGCAAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATTCATGCTGGt             <  2:620288/141‑1 (MQ=255)
                                                                                                                                 tACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATTCACGCTGGTa            <  3:559528/141‑1 (MQ=255)
                                                                                                                                      ttAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATTCATGCTGGTATcaca       <  2:405151/141‑1 (MQ=255)
                                                                                                                                      ttAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATTCATGCTGGTATcaca       >  1:289867/1‑141 (MQ=255)
                                                                                                                                         ggCCAAAAAAAGGACTTAACAATGTGGGAGGTAGTGAATCAGCAGTTTGCTGAATATGCAGTGAAAGTAACTGATAAACGTGCACATTTTAATCAAGAGCTAGAGTCGTTAGCTAAACCGATACATGCTGGTCTCAGAAAt    >  2:825057/1‑141 (MQ=255)
                                                                                                                                           ccAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATTCATGCTGGTATCACAAATGa  >  4:675730/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                            
TAAAAGGCTCTCCTCAAATAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATTCATGCTGGTATCACAAATGA  >  USA300TCH1516_ALE/4309‑4588

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