Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 933,033 T→C intergenic (‑74/‑27) nfuA ← / → USA300HOU_RS04645 Fe/S biogenesis protein NfuA/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE933,0330TC98.9% 282.8 / NA 87intergenic (‑74/‑27)nfuA/USA300HOU_RS04645Fe/S biogenesis protein NfuA/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (37/49);  total (37/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGACGAAGACGTTCAATAACTTCTGCTACTTGATCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATATTTATGACTGGCTACAGCCGCTATTAAAAAG  >  USA300TCH1516_ALE/932897‑933169
                                                                                                                                        |                                                                                                                                        
tGGACGAAGACGTTCAATAACTTCTGCTACTTGATCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTa                                                                                                                                       <  2:438534/140‑1 (MQ=255)
    cgcaGACGTTCAATAACTTCTGCTACTTGATCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCaa                                                                                                                                   <  2:820095/137‑1 (MQ=255)
       aGACGTTCAATAACTTCTGCTACTTGCTCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATACAAACCTATCAAAcc                                                                                                                                >  2:152278/1‑140 (MQ=255)
         aCGATCAATCACTTCTGCTACTTGATCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCaa                                                                                                                              <  2:613394/140‑1 (MQ=255)
         aCGATCAATAACTTCTGCTACTGGATCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGACAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCaa                                                                                                                              <  2:613389/140‑1 (MQ=255)
          cGTTCAATGACTTCTGCTACTTGATCAAACATCGTTGTATCTTCAGTAGGAATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATGATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAAt                                                                                                                             >  2:741182/1‑140 (MQ=255)
           gTTCAATAACTTCTGCTACTTGATCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATg                                                                                                                            <  2:402173/140‑1 (MQ=255)
             tCAATAACTTCTGCTACTTGATCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATggg                                                                                                                          <  2:591246/140‑1 (MQ=255)
             tCAATAACTTCTGCTACTTGATCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATAAAACCAATggg                                                                                                                          >  2:810807/1‑140 (MQ=255)
                  aaCTTCTGCTACTTGATCAAACATCGCTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAACCTTTATTAGTATCATTTAATATATAGTTTAATAAAAACCTATCAAACCAATAGGGGATa                                                                                                                     >  2:601385/1‑140 (MQ=255)
                    cTTCTGCTACTTGATCAAACATCGTAGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATaaa                                                                                                                   <  1:606802/140‑1 (MQ=255)
                       cTGCTTCTTGATCAAACCTCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTCAATGATACCCATCAAACCAATGTGGGATAAatat                                                                                                                <  2:664931/140‑1 (MQ=255)
                       cTGATACTTGATCAAAAATCGTTGTATCTGCAGTAGGAATCAGGCAAGTCTCCGTACGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAatat                                                                                                                <  1:20536/140‑1 (MQ=255)
                       atGCTACTTGATCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAatat                                                                                                                <  1:71686/139‑1 (MQ=255)
                        tGCTACTTGATCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCACCTTTCGAAAGATGAATTTGATAATCTTTATTAGTATATAATAATAGATAGTGTAATAAAAACCTATCAAACCAAAGGGGGATAAataaa                                                                                                               >  2:76406/1‑138 (MQ=255)
                               tGATCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAatatat                                                                                                              >  1:30318/1‑134 (MQ=255)
                               tGATCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAatatat                                                                                                              <  2:30318/134‑1 (MQ=255)
                                      aCATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAgt                                                                                                 <  2:347559/140‑1 (MQ=255)
                                      aCATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAgt                                                                                                 <  1:4276/140‑1 (MQ=255)
                                           gTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGt                                                                                            >  2:414478/1‑140 (MQ=255)
                                            ttGTTTCTTGAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATTTATGGAGCACGTGAGTGTGGTa                                                                                           <  2:125470/140‑1 (MQ=255)
                                               tATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATATTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAAGATATGGAGAACGTGAGTGTCGTAGtt                                                                                        >  2:654419/1‑140 (MQ=255)
                                                   agctgTAGGCATCAGGCAATTCTCCTTTGGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAGAAACCTTTCAAACCAATGGGGGAGAAATATATGGGGCAGGTGGGTGTGGTAGTTTATg                                                                                    <  2:76820/136‑1 (MQ=255)
                                                     cAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATggg                                                                                  <  2:326744/140‑1 (MQ=255)
                                                     cAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATggg                                                                                  >  2:855451/1‑140 (MQ=255)
                                                      aGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATCATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATgggg                                                                                 <  1:77653/140‑1 (MQ=255)
                                                         aGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAg                                                                              <  1:586173/140‑1 (MQ=255)
                                                          ggCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGa                                                                             >  2:595736/1‑140 (MQ=255)
                                                          ggCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGGGGTAGTTTATGGGGCAGa                                                                             >  2:201066/1‑140 (MQ=255)
                                                                 ggCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGCGGGGGGTAGTTTATGGGGCAGATGTtgta                                                                      >  1:427909/1‑140 (MQ=255)
                                                                 ggCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTGGTTTATGGGGCAAAAGTtata                                                                      >  1:428014/1‑140 (MQ=255)
                                                                 ggCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTtata                                                                      >  2:596654/1‑140 (MQ=255)
                                                                 ggCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGAGAAATATATGGAGCACGTGAGGGTGGTAGTTTATGGGGCAGATGTtata                                                                      >  2:236459/1‑140 (MQ=255)
                                                                   cAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATAtg                                                                    <  1:452590/140‑1 (MQ=255)
                                                                     aGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGCTACATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATAtgtg                                                                  >  1:7855/1‑140 (MQ=255)
                                                                      gTCTCCTTTCGAGAGATTAATTTTATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGc                                                                 <  2:94421/140‑1 (MQ=255)
                                                                       tctcCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCa                                                                <  1:1006078/140‑1 (MQ=255)
                                                                          ccTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGt                                                             <  1:758273/140‑1 (MQ=255)
                                                                          ccTTTCGAGAAATGAATTTGAAAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGt                                                             <  1:50540/140‑1 (MQ=255)
                                                                               cGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGt                                                        <  1:780099/140‑1 (MQ=255)
                                                                               cGAGAGATGAATTTGATAATCTTTATTAGCACGTTAAAAGAGGTAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGt                                                        <  1:774522/140‑1 (MQ=255)
                                                                                  gagaTGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGAAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTaa                                                     >  2:415270/1‑140 (MQ=255)
                                                                                         aTTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCa                                              <  2:820734/140‑1 (MQ=255)
                                                                                           ttGATAATTTTTATTAGCAGATTATAATAGATAGGTTAATAAAAACCTATCAAACCAATGGGGGATAAATATAGGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAACGCGCCAAc                                            <  2:33711/140‑1 (MQ=255)
                                                                                           ttGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAAc                                            <  2:350059/140‑1 (MQ=255)
                                                                                            ttaTAATCTTTATTAGGATTTTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACa                                           <  1:92233/138‑1 (MQ=255)
                                                                                              aTAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATc                                         >  2:208141/1‑140 (MQ=255)
                                                                                                actCTTTATTAATATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTGGCGTTAATGCGCCAACATCGa                                       <  2:85035/138‑1 (MQ=255)
                                                                                                aaTCTTTATTAGTATATTATAATAGATTGTGTAATAAAAACCTACCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGa                                       <  1:639867/140‑1 (MQ=255)
                                                                                                aaTCTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGa                                       <  1:713172/140‑1 (MQ=255)
                                                                                                aaTCTTGATTAGTATATTATAATAGATAGTGTAATAAAATCCTATCAAACCAATGGGGGATAAAGAGATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGa                                       <  1:575077/140‑1 (MQ=255)
                                                                                                   cTTTATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGAGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACGTCGCAAg                                    >  1:818507/1‑140 (MQ=255)
                                                                                                    tttATTAGTATATTATAATAGTTAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGATTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGa                                   >  1:522013/1‑140 (MQ=255)
                                                                                                    tttATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGa                                   <  2:429140/140‑1 (MQ=255)
                                                                                                    tttATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGa                                   >  2:94854/1‑140 (MQ=255)
                                                                                                    tttATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGa                                   >  1:395285/1‑140 (MQ=255)
                                                                                                    tttATTAGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGAGCCAACATCGAAAGa                                   >  1:395202/1‑140 (MQ=255)
                                                                                                    tctattaGCATATTATAATAGATAGTTTAATAAAAACCTATCAAACCAATGGGGGATAACTATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGa                                   <  1:519096/138‑1 (MQ=255)
                                                                                                         tggTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAGTGGGGGATAAATATATGGAGCACGTGCGTGCGGTAGTTTATGGGGCAGATGTTATGTGTGCAAGTTGGGTTAATGCGCCAACATCGAAAGATAttt                              >  2:623227/3‑140 (MQ=255)
                                                                                                         taGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATAttt                              >  2:393875/1‑140 (MQ=255)
                                                                                                          aGTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATATTTa                             <  1:178467/140‑1 (MQ=255)
                                                                                                           gTATATTATAATAGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATATTTAt                            >  1:872953/1‑140 (MQ=255)
                                                                                                                 aatattcgatagTGAAATAAAAACCTATCAAGCCATGGGGGGGTAAATATAGGGTGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCGAACATCGAAAGCTATTTATGACTgg                      <  2:174547/133‑1 (MQ=255)
                                                                                                                  ataataGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGGGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATATTTATGACTGGc                     <  2:927685/140‑1 (MQ=255)
                                                                                                                  ataataGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATATTTATGACTGGc                     <  1:187496/140‑1 (MQ=255)
                                                                                                                  ataataGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATATTTATGACTGGc                     <  2:835768/140‑1 (MQ=255)
                                                                                                                  ataataGATAGTGCAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATATTTATGACTGGc                     <  1:244514/140‑1 (MQ=255)
                                                                                                                    aataGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATATTTATGACTGGCTa                   >  2:379170/1‑140 (MQ=255)
                                                                                                                    aataGATAGTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGAGTGGTAGTGTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATATTTATGACTGGCTa                   >  2:534163/1‑140 (MQ=255)
                                                                                                                    aataGATAGTGTAATAAAAACCTATCAAACCAATGAGGGATAAATATATAGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGAGTTAATGCGCCAACATCGAAAGAAATTTATGAATGGCTa                   >  2:39438/1‑140 (MQ=255)
                                                                                                                     atagatagTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATATTTATGACTGGCTAc                  >  1:658587/1‑140 (MQ=255)
                                                                                                                      tagatagTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATATTTATGACTGGCTACa                 <  1:347663/140‑1 (MQ=255)
                                                                                                                      aagatagTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATATTTATGACTGGCTACa                 <  1:604877/139‑1 (MQ=255)
                                                                                                                         atagTGTAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTCTGGGGCAGATGTTATCTGTGCAAGTTGCGTTAATGCGCCAACATCGAACGATATTTATGACTGGCTACAGcc              <  2:586073/140‑1 (MQ=255)
                                                                                                                          ttgtggAATAAAAATCTATCAACCCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATATTTATGACTGGCTACAGCCg             <  2:86115/138‑1 (MQ=255)
                                                                                                                          tagTGTAATAAAAACCTATCAAACCAAAGGGAGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATATTTATGACTGGCTACAGCCg             <  1:379264/140‑1 (MQ=255)
                                                                                                                          tagTGTAACTAAAACCGTTCAAACCAATGGTGGATAAATAGATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATATTTATGACTGGCTACAGCCg             <  1:76406/140‑1 (MQ=255)
                                                                                                                            gtgtAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATATTTATGACTGGCTACAGCCGCt           >  1:663536/1‑140 (MQ=255)
                                                                                                                            gtgtAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATATTTATGACTGGCTACAGCCGCt           >  2:891703/1‑140 (MQ=255)
                                                                                                                            gtgtAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACAGCGAAAGATATTTATGACTAGCTACAGCCGCt           >  1:104122/1‑140 (MQ=255)
                                                                                                                             tgtAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATATTTATGACTGGCTACAGCCGCTa          >  1:144135/1‑140 (MQ=255)
                                                                                                                             tgtAATAAAAACCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTCTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATATTTATGACTGGCTACAGCCGCtc          >  1:675186/1‑139 (MQ=255)
                                                                                                                                    aaaaCCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATATTTATGACTGGCTACAGCCGCTATTaaaaa   >  2:286027/1‑140 (MQ=255)
                                                                                                                                    aaaaCCTATCAAACAAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATATTTATGACTGGCTACAGCCGCTATTaaaaa   <  2:146852/140‑1 (MQ=255)
                                                                                                                                     aaaCCTATCAAACCAATGGGGGATAAATGTATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATATTTATGACTGGCTGACGCCGCTTTTAAagag  >  2:867088/1‑140 (MQ=255)
                                                                                                                                     aaaCCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATATTTATGACTGGCTACAGCCGCTATTAAAAAg  <  2:671028/140‑1 (MQ=255)
                                                                                                                                     aaaCCTATCAAACCAATGGGGGATAAAGATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTCGCGTTAATGCGCCAACATCGAAAGATATTTATGACTGGCTACAGCCGCTATTAAAAAg  <  1:602719/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                        
TGGACGAAGACGTTCAATAACTTCTGCTACTTGATCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAGATGTTATATGTGCAAGTTGCGTTAATGCGCCAACATCGAAAGATATTTATGACTGGCTACAGCCGCTATTAAAAAG  >  USA300TCH1516_ALE/932897‑933169

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: ‑