Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,771,843 C→A 8.3% F214L (TTC→TTA)  phoB → Alkaline phosphatase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,771,8430CA 8.3% 262.6 / 3.0 85F214L (TTC→TTA) phoBAlkaline phosphatase 3
Reads supporting (aligned to +/- strand):  ref base C (20/55);  new base A (3/6);  total (24/61)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGACTTCTAGAGATGATAAAAACGAAATCGCTAAACAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAATATGCCACTTGCAATAGATGCTTCTAAAGACGAACCATCATTAGCT  >  USA300TCH1516_ALE/2771703‑2771978
                                                                                                                                            |                                                                                                                                       
gTGACTTCTAGAGATGATAAAAACGAAATCGCTAAACAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTc                                                                                                                                         >  1:7686/1‑141 (MQ=255)
gTGACTTCTAGAGATGATAAAAACGAAATCGCTAAACAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCAGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTc                                                                                                                                         >  2:645449/1‑141 (MQ=255)
      tCTAGAGATGATAAAAACGAAATCGCTAAACAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCaaaaaa                                                                                                                                   <  1:470328/141‑1 (MQ=255)
      tCTAGAGAGGATAAAAACGAAATCGCTAAACAATTTTAGAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCaaaaaa                                                                                                                                   <  1:77263/141‑1 (MQ=255)
      tCTAGAGAGAATAAAAACCAAATCGCTAAACAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCaaaaaa                                                                                                                                   <  1:483346/141‑1 (MQ=255)
        tAGAGATGATAAGAACGAAATCGCTAAAAAAATTTATAAAGATAAAATGGACGGAAAGCATAAAGTGGGTGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAGAAGGAGTGGCTATGTTGAAAAAAAAAAAAAAAAAGa                                                                                                                                 <  2:508066/141‑1 (MQ=255)
           agaTGATAAAAACGAAATCGCTAAACAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGTAATTTAGATAAAAAATTCAAAAAAGATgg                                                                                                                              <  1:654600/141‑1 (MQ=255)
            gatgatAAAAACGAAATCGCTAAACAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGt                                                                                                                             >  2:78964/1‑141 (MQ=255)
            gatgatAAAAACGAAATCGCTAAACAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGt                                                                                                                             >  2:843091/1‑141 (MQ=255)
               gaaaaaaaCGAAATCGATAAACCATTTTATAAAGATAAAAACAACGGGAAGCATTAAGTTGATGTTTTATGAGGCGGCGGCGCTAAATATTTCACTAAAAGTAGTGGCAATTTATATAAAAAATTCAAAAAAGATGGTTAc                                                                                                                          <  2:610910/138‑1 (MQ=255)
                  aaaaaCGAAATCGCTAAACAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGAc                                                                                                                       <  1:239578/141‑1 (MQ=255)
                   aaaaCGAAATCGCTAAACAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGAcc                                                                                                                      >  1:103737/1‑141 (MQ=255)
                     aaCGAAATCGCTAAACAATTTTATAAAGATAAAATCAACGGAAAGCAGAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCtt                                                                                                                    >  1:485846/1‑141 (MQ=255)
                      aCGAAATCGCTAAACAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAAGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGACAATTTAGATAAAAAATGCAAAAAAGATGGTTCAAACCTTg                                                                                                                   >  2:306974/1‑141 (MQ=255)
                       cGAAATCGCTAAACAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGc                                                                                                                  <  1:580604/141‑1 (MQ=255)
                       cGAAATCGCTAAACAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGc                                                                                                                  <  2:764243/141‑1 (MQ=255)
                       cGAAATCGCTAAACAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGc                                                                                                                  <  1:374948/141‑1 (MQ=255)
                       cGAAAGCGCTAAACGATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGc                                                                                                                  <  2:457825/141‑1 (MQ=255)
                        gAAATCTCTATACAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGcc                                                                                                                 <  2:192193/141‑1 (MQ=255)
                        gAAATCGCTAAACAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGATTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAATGGTGAAAAAGATGGACATGACCTTGcc                                                                                                                 >  2:947019/1‑141 (MQ=255)
                         aaaTCGCTAAACAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCa                                                                                                                >  2:938614/1‑141 (MQ=255)
                          aaTCGCTAAACAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCAc                                                                                                               >  2:740711/1‑141 (MQ=255)
                          aaTCGCTAAACAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTAAAAAAAGATGGTTACGACCTTGCCAc                                                                                                               >  2:266339/1‑141 (MQ=255)
                           aTCGCTAAACAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACg                                                                                                              <  1:644633/141‑1 (MQ=255)
                              gCTAAACAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAAc                                                                                                           <  2:520424/141‑1 (MQ=255)
                               cTAAACAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGAAGTTCATCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACa                                                                                                          <  1:686478/141‑1 (MQ=255)
                               cTAAACAATTTTAAAAAGATAAAATCAACGGAAAGCATAAAGGTGATGTTTTATTAGGCGGCGGCGCGAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACa                                                                                                          <  1:686750/141‑1 (MQ=255)
                                   aCAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTaa                                                                                                      <  1:809493/141‑1 (MQ=255)
                                    cAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGTGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAg                                                                                                     >  2:260530/1‑141 (MQ=255)
                                      aTTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGa                                                                                                   <  2:510218/141‑1 (MQ=255)
                                         ttATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGCCCTTGCCCCGAACAGTCAGGAGtt                                                                                                >  1:336097/1‑141 (MQ=255)
                                            taAAGATAAAATCAACGGAAACCATAAAGTTGATGTTTTATTAGGCGCCGGCGCTAAATATTTCGGTAAAAGTAATTGCAATTTAGATAAAAAATTCAAAAAAGATGTTTACTACCTTGCCACGAACAGTAAGGAGTTATc                                                                                             <  1:318252/141‑1 (MQ=255)
                                                 aTAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATACAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAAt                                                                                        <  2:64578/141‑1 (MQ=255)
                                                  tAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATc                                                                                       >  2:796786/1‑141 (MQ=255)
                                                  aaaaaTCAACGGAAATGATAAAGTTGATGTTTTATTAGGAGGCGGAGCTAAATATTTAGGTAAAAGTAATGGCAATTTAGAAAAAAAAATAAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATc                                                                                       <  2:101278/140‑1 (MQ=255)
                                                   aagagaaaCGGAAAGAAAAAAGTGGATGATTTATTAGGCGGCGGCGCAAAATATGTCAGTAAAAGTAATGGCAATTGAGAAAAAAAAAAAAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCt                                                                                      <  2:208833/135‑1 (MQ=255)
                                                   aaaaACAAAGGAAAGCATAAAGTGGATGTCTGATTAGGCGGCGGCGATAAATATTTCGGTAAAAGTAATGGCAATTTAGAAAAAAAAAAAAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCt                                                                                      <  2:246988/141‑1 (MQ=255)
                                                     aaTCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTga                                                                                    >  1:251860/1‑141 (MQ=255)
                                                          aCGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTgataaag                                                                               >  1:715919/1‑141 (MQ=255)
                                                             gAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTgataaagata                                                                            <  1:4770/141‑1 (MQ=255)
                                                             gAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTAAAAAAAGGCGGTTAGGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTgataaagata                                                                            >  1:261768/1‑141 (MQ=255)
                                                             gAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATGTCGGTAAAAGGAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTgataaagata                                                                            >  1:141691/1‑141 (MQ=255)
                                                             gAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAAGATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTgataaagata                                                                            <  1:22298/141‑1 (MQ=255)
                                                              aaaGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTgataaagataa                                                                           <  2:461033/141‑1 (MQ=255)
                                                              aaaGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTgataaagataa                                                                           <  1:297636/141‑1 (MQ=255)
                                                                    tAAGGTTGATGTTTGTTTGGGCAGGGGCGCTGACTATTTAGGTAACGGTGATTGCATTTTAAATAAAAAATTCAAAGAAGATGGTTACTTCCTGTGCACGAACAGGAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCt                                                                     <  2:786065/141‑1 (MQ=255)
                                                                    tAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCTCGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCt                                                                     <  1:648621/141‑1 (MQ=255)
                                                                    aaaaGTTGATGTTTTATTAGGCGGCGGGGCTGAAAATTCCTGTAAAAGTAATAGCAAATTAGATAAAAAATTCAAAAAATATGTGTACGACCTTTCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCt                                                                     <  2:785983/140‑1 (MQ=255)
                                                                       aGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCCGATAAAGATAAAGTTCTTgg                                                                  <  1:662386/141‑1 (MQ=255)
                                                                       aGTTGATGTTTTATTAGGCGGCCGCGCTAAAGAGTTCGGAAAACGATTTGGCGAATGAAATAAAAAATTCAAAAGAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTgg                                                                  <  1:213565/141‑1 (MQ=255)
                                                                        gTTGATGTTTTATTAGGTGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGgt                                                                 <  2:841218/141‑1 (MQ=255)
                                                                         ttGATGTTTTATTAGGCGGCGACGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGgtt                                                                >  1:485014/1‑141 (MQ=255)
                                                                              gTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGgtttgttt                                                           <  1:657812/141‑1 (MQ=255)
                                                                              gTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGgtttgttt                                                           <  1:158655/141‑1 (MQ=255)
                                                                                 ttataaGGCGACGGCGATAAATATTTAGGAAAAAGTAATGGCAATTTAAAAAAAAAAAAAAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCa                                                        <  2:113549/141‑1 (MQ=255)
                                                                                     taGGCGGCGGTGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATa                                                    <  1:260530/141‑1 (MQ=255)
                                                                                       ggcggcggcGCTAAATATTTTGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATaaa                                                  <  2:683077/141‑1 (MQ=255)
                                                                                       ggcggcggcGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGATCAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATaaa                                                  <  1:567620/141‑1 (MQ=255)
                                                                                       ggcggcggcGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATaaa                                                  <  1:536875/141‑1 (MQ=255)
                                                                                        gcggcggcgCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATGTGATCAGGATAAAGTTCTTGGTTTGTTTGCAGTTaata                                                 >  1:745336/1‑139 (MQ=255)
                                                                                           gcggcgCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAata                                              <  1:78964/141‑1 (MQ=255)
                                                                                            cggcgCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAatat                                             <  2:898100/141‑1 (MQ=255)
                                                                                              gcgcTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGCTGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGCTAAAAATATGc                                           <  1:947019/141‑1 (MQ=255)
                                                                                              gcgcTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTGTGTTCGCAGATAAAAATATGc                                           <  1:484164/141‑1 (MQ=255)
                                                                                              gcgcTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTCAAGAGCTTTTCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAATATGc                                           <  1:449089/141‑1 (MQ=255)
                                                                                               cgcTAAATAGTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAATATGcc                                          <  2:715919/141‑1 (MQ=255)
                                                                                                gcTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTAGAAAAAGTCGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAATATGCCa                                         >  1:366242/1‑141 (MQ=255)
                                                                                                  tAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAATATGCCACt                                       <  1:843091/141‑1 (MQ=255)
                                                                                                  aaaaTATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAAGTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAATATGCCACt                                       <  2:772480/140‑1 (MQ=255)
                                                                                                        tttCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAATATGCCACTTGCAAt                                 >  2:436391/1‑141 (MQ=255)
                                                                                                          tCGTTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAATATGGTTACTACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAATATGCCACTTGCAATAg                               <  2:359642/141‑1 (MQ=255)
                                                                                                          tCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAATATGCCACTTGCAATAg                               >  1:168614/1‑141 (MQ=255)
                                                                                                            ggTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAATATGCCACTTGCAATAGAt                             <  1:654521/141‑1 (MQ=255)
                                                                                                            ggTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCAACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAATATGACGCTTGCAATAGAt                             >  2:810748/1‑141 (MQ=255)
                                                                                                            cgTAAAAGAAATGGCAATTTAGATAAAAAATTCAAAAAAGATTGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAATATGCCACTTGCAATAGAt                             <  1:350522/140‑1 (MQ=255)
                                                                                                              tAAATGTAATACCAATATAGATAAAATATGCAAACAAGATTGTTTCCACCTTGACCCCAACAGTAAGGAGTTATCAGAATATGATGAAGATCAATTTCTTGGTTTGTTTGCAGATAAACATATGTCACTTGCAATAGATGc                           <  2:287237/141‑1 (MQ=255)
                                                                                                              tAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAATATGCCACTTGCAATAGATGc                           <  1:167783/141‑1 (MQ=255)
                                                                                                                 aaGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAATATGCCACTTGCAATAGATGCttt                        >  2:845331/1‑140 (MQ=255)
                                                                                                                  aGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAATATGCCACTTGCAATAGATGCTTCt                       <  1:321641/141‑1 (MQ=255)
                                                                                                                       tGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGCCCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAATATGCCACTTGCAATAGATGCTTCTAAAGa                  <  1:727458/141‑1 (MQ=255)
                                                                                                                       tGGCAATTTAGATAAAAAACTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAATATGCCACTTGCAATAGATGCTTCTAAAGa                  <  1:910441/141‑1 (MQ=255)
                                                                                                                        ggCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAATATGCCACTTGCAATAGATGCTTCTAAAGAc                 <  2:7686/141‑1 (MQ=255)
                                                                                                                                 gATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAATATGACACTTGCAATAGATGCTTCTAAAGAAAAACcttca        >  2:783711/1‑141 (MQ=255)
                                                                                                                                   aaaaaaaTTCAAAAAAGATGGTTACGACCTTGCCAAGAACAGTAAGGAGTTAGCAAAATCTGATTAAGATAAATTTCTTTGTTTGTTTGCAGATAAAAATATGCCACTTGCAATAGATGCTTCTAAAGACGAACCATCAtt      <  2:336097/140‑1 (MQ=255)
                                                                                                                                    aaaaaaTTAAAAAAATATAGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAATATGCCACTTGCAATAGATGCTTCTAAAGACGAACCATCATTa     <  2:103737/141‑1 (MQ=255)
                                                                                                                                       aaaTTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAATATGCCACTTGCAATAGATGCTTCTAAAGACGAACCATCATTAGCt  <  1:645449/141‑1 (MQ=255)
                                                                                                                                       aaaTTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAATATGCCACTTGCAATAGATGCTTCTAAAGACGAACCATCATTAGCt  <  1:796786/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                       
GTGACTTCTAGAGATGATAAAAACGAAATCGCTAAACAATTTTATAAAGATAAAATCAACGGAAAGCATAAAGTTGATGTTTTATTAGGCGGCGGCGCTAAATATTTCGGTAAAAGTAATGGCAATTTAGATAAAAAATTCAAAAAAGATGGTTACGACCTTGCCACGAACAGTAAGGAGTTATCAAAATCTGATAAAGATAAAGTTCTTGGTTTGTTTGCAGATAAAAATATGCCACTTGCAATAGATGCTTCTAAAGACGAACCATCATTAGCT  >  USA300TCH1516_ALE/2771703‑2771978

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