Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 308,084 C→T 5.8% A111V (GCA→GTA)  lytR_1 → Sensory transduction protein LytR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE308,0840CT 5.8% 1057.5 / 36.2 307A111V (GCA→GTA) lytR_1Sensory transduction protein LytR
Reads supporting (aligned to +/- strand):  ref base C (153/136);  new base T (8/10);  total (161/146)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 310 total aligned reads displayed.

CTGCACATGACCAATACGCAGTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCGGCGAATTTTGATCAAAGCTT  >  USA300TCH1516_ALE/307994‑308174
                                                                                          |                                                                                          
 tGCACATGACCAATACGCAGTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCa                                                                                           <  1:4290316/91‑1 (MQ=255)
   cacaTGACCAATACGCAGTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAagtagt                                                                                         <  1:3248286/91‑1 (MQ=255)
     caTGACCAATACGCAGTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCa                                                                                       <  1:2339382/91‑1 (MQ=255)
       tGACCAATACGCAGTACAGGCATTTGAATTAAATGCGACAGAGTATATTTTGAAACCGTTTGGTCAAAAACGTTTTGAACAAGCAGTCAAt                                                                                     <  1:5688428/91‑1 (MQ=255)
       tGACCAATACGCAGTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAAt                                                                                     >  1:3513555/1‑91 (MQ=255)
         aCCAATACGCAGTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATaa                                                                                   <  1:1693242/91‑1 (MQ=255)
           cAATACGCAGTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGTAGTCAATAAAg                                                                                 >  1:933882/1‑91 (MQ=255)
           cAATACGCAGTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAg                                                                                 <  1:67772/91‑1 (MQ=255)
           cAATACGCAGTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAg                                                                                 <  1:5625200/91‑1 (MQ=255)
           cAATACGCAGTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAg                                                                                 <  2:1737638/91‑1 (MQ=255)
           cAATACGCAGTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAg                                                                                 <  1:1675461/91‑1 (MQ=255)
           cAATACGCAGTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAg                                                                                 >  2:3541284/1‑91 (MQ=255)
           cAATACGCAGTACAGGCATATGAATTAAATGCGACAGACTATATTTTGAAAAAAAATAAAAAAAAACGTATTGAACAAGCAGTCAATAAAg                                                                                 <  2:3845703/91‑1 (MQ=255)
            aaTACGCAGTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAgt                                                                                >  1:1446238/1‑91 (MQ=255)
            aaTACGCAGTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAgt                                                                                >  1:3390063/1‑91 (MQ=255)
               aCGCAGTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGATACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAgtgcg                                                                             <  2:5669699/91‑1 (MQ=255)
               aCGCAGTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAgtgcg                                                                             >  2:5867135/1‑91 (MQ=255)
                cGCAGTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAgtgcgt                                                                            >  2:1426847/1‑91 (MQ=255)
                cGCAGTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAgtgcgt                                                                            <  1:46568/91‑1 (MQ=255)
                 gCAGTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAgtgcgtg                                                                           >  1:259716/1‑91 (MQ=255)
                 gCAGTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAgtgcgtg                                                                           >  1:1367963/1‑91 (MQ=255)
                 gCAGTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAgtgcgtg                                                                           >  2:4001307/1‑91 (MQ=255)
                  cAGTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGTAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAgtgcgtgc                                                                          >  2:3960133/1‑91 (MQ=255)
                  cAGTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAgtgcgtgc                                                                          <  1:2189999/91‑1 (MQ=255)
                    gTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGa                                                                        >  1:5416326/1‑91 (MQ=255)
                      aCAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACt                                                                      <  1:889180/91‑1 (MQ=255)
                       cAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTa                                                                     <  2:2100103/91‑1 (MQ=255)
                         ggCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGTAGTCAATAAAGTGCGTGCGACTaaa                                                                   >  2:4558562/1‑91 (MQ=255)
                          gCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAg                                                                  >  1:3305822/1‑91 (MQ=255)
                          gCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAg                                                                  <  2:4875695/91‑1 (MQ=255)
                           cATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGc                                                                 <  2:3674040/91‑1 (MQ=255)
                           cATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGc                                                                 <  1:1307078/91‑1 (MQ=255)
                            aTTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGcc                                                                <  1:3861321/91‑1 (MQ=255)
                            aTTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGcc                                                                <  2:1533271/91‑1 (MQ=255)
                             tttGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCa                                                               <  1:2223015/91‑1 (MQ=255)
                                 aaTTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAga                                                           <  2:2138294/91‑1 (MQ=255)
                                  aTTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAgat                                                          <  2:4381240/91‑1 (MQ=255)
                                   ttAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAgatg                                                         >  2:1660171/1‑91 (MQ=255)
                                   ttAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAgatg                                                         <  1:783072/91‑1 (MQ=255)
                                    tAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAgatga                                                        <  2:5631407/91‑1 (MQ=255)
                                      aaTGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGata                                                      <  2:1610902/91‑1 (MQ=255)
                                       aTGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGataa                                                     <  1:18803/91‑1 (MQ=255)
                                        tGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGataat                                                    >  2:671325/1‑91 (MQ=255)
                                         gCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGataata                                                   >  1:5192185/1‑91 (MQ=255)
                                           gACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAAc                                                 >  1:4738028/1‑91 (MQ=255)
                                            aCAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAGTAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACg                                                >  1:4851580/1‑91 (MQ=255)
                                            aCAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGTGTGCGACTAAAGCCAAAGATGATAATAACg                                                <  2:4301080/91‑1 (MQ=255)
                                            aCAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACg                                                >  2:2248627/1‑91 (MQ=255)
                                                aCTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAATGATGATAATAACGCAAg                                            <  1:4439307/91‑1 (MQ=255)
                                                aCTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAg                                            <  1:4956384/91‑1 (MQ=255)
                                                aCTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAg                                            <  2:5745350/91‑1 (MQ=255)
                                                aCTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAg                                            <  1:4067081/91‑1 (MQ=255)
                                                   atatTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGc                                         <  1:4359018/91‑1 (MQ=255)
                                                    tatTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCa                                        >  1:4473845/1‑91 (MQ=255)
                                                     atTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCaa                                       >  1:1248106/1‑91 (MQ=255)
                                                       tttGAAACCGTTTGGTCAAAAACGTATTGAATAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAAtt                                     <  1:5454412/91‑1 (MQ=255)
                                                       tttGAAACCGTTTGGTCAAAAACGTATTGAACAAGTAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAAtt                                     <  2:3093861/91‑1 (MQ=255)
                                                       tttGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAAtt                                     >  1:3096980/1‑91 (MQ=255)
                                                       tttGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAAtt                                     <  2:1152817/91‑1 (MQ=255)
                                                       tttGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAAtt                                     <  1:700172/91‑1 (MQ=255)
                                                       tttGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAAtt                                     >  2:1804142/1‑91 (MQ=255)
                                                        ttGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTg                                    >  2:1085811/1‑91 (MQ=255)
                                                         tGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGc                                   >  2:2522472/1‑91 (MQ=255)
                                                         tGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGc                                   >  2:5381179/1‑91 (MQ=255)
                                                          gAAACCGTTTGGTCAAAAACGTATTGAACAAGTAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCg                                  >  1:4822395/1‑91 (MQ=255)
                                                           aaaCCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAGTAACGCAAGTGCAATTGCGa                                 >  1:2932247/1‑91 (MQ=255)
                                                           aaaCCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGa                                 >  2:1234029/1‑91 (MQ=255)
                                                           aaaCCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGa                                 <  1:2372870/91‑1 (MQ=255)
                                                           aaaCCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGa                                 >  1:1573198/1‑91 (MQ=255)
                                                            aaCCGTTTGGTCAAAAACGTATTGAACAAGTAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGaa                                >  2:4792525/1‑91 (MQ=255)
                                                             aCCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAAt                               >  2:301513/1‑91 (MQ=255)
                                                                gTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCGATTGCGAATGat                            >  2:3749162/1‑91 (MQ=255)
                                                                gTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGat                            >  1:5796611/1‑91 (MQ=255)
                                                                 tttGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGata                           >  1:2837872/1‑91 (MQ=255)
                                                                     gTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTc                       >  2:1901739/1‑91 (MQ=255)
                                                                      tCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGTATGATATGTCg                      <  2:5822483/91‑1 (MQ=255)
                                                                       cAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCgg                     <  1:2248627/91‑1 (MQ=255)
                                                                         aaaaCGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCGGCg                   >  1:3921919/1‑91 (MQ=255)
                                                                          aaaCGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCGGCGa                  >  2:50125/1‑91 (MQ=255)
                                                                          aaaCGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCGGCGa                  >  1:193766/1‑91 (MQ=255)
                                                                          aaaCGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCGGCGa                  >  2:904110/1‑91 (MQ=255)
                                                                              gTATTGAACAAGTAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCGGCGAAttt              >  1:5299417/1‑91 (MQ=255)
                                                                              gTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCGGCGAAttt              <  1:4143672/91‑1 (MQ=255)
                                                                               tATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCGGCGAAtttt             >  2:3610932/1‑91 (MQ=255)
                                                                                 ttGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCGGCGAATTTTGa           >  1:369170/1‑91 (MQ=255)
                                                                                 ttGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCGGCGAATTTTGa           >  2:2746275/1‑91 (MQ=255)
                                                                                 ttGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCGGCGAATTTTGa           <  1:5381179/91‑1 (MQ=255)
                                                                                  tGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCGGCGAATTTTGAt          <  2:1107994/91‑1 (MQ=255)
                                                                                   gAACAAGTAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCGGCGAATTTTGATc         >  2:3336010/1‑91 (MQ=255)
                                                                                   gAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCGGCGAATTTTGATc         <  1:2747792/91‑1 (MQ=255)
                                                                                   gAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCGGCGAATTTTGATc         <  1:5098602/91‑1 (MQ=255)
                                                                                     aCAAGCAGTCGATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCGGCGAATTTTGATCaa       <  2:4671506/91‑1 (MQ=255)
                                                                                     aCAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCGGCGAATTTTGATCaa       <  1:421058/91‑1 (MQ=255)
                                                                                       aaGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCGGCGAATTTTGATCAAAg     <  1:1097262/91‑1 (MQ=255)
                                                                                        aGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCGGCGAATTTTGATCAAAGc    <  2:4479752/91‑1 (MQ=255)
                                                                                         gCAGTCAATAAAGTGCGTGCGAGTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCCGCGAATTTTGATCAAAGCt   >  1:4507074/1‑91 (MQ=255)
                                                                                         gCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCGGCGAATTTTGATCAAAGCt   >  2:3663101/1‑91 (MQ=255)
                                                                                         gCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCGGCGAATTTTGATCAAAGCt   >  2:3259538/1‑91 (MQ=255)
                                                                                         gCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCGGCGAATTTTGATCAAAGCt   >  2:2975011/1‑91 (MQ=255)
                                                                                          cAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCGGCGAATTTTGATCAAAGCtt  <  1:3731519/91‑1 (MQ=255)
                                                                                          |                                                                                          
CTGCACATGACCAATACGCAGTACAGGCATTTGAATTAAATGCGACAGACTATATTTTGAAACCGTTTGGTCAAAAACGTATTGAACAAGCAGTCAATAAAGTGCGTGCGACTAAAGCCAAAGATGATAATAACGCAAGTGCAATTGCGAATGATATGTCGGCGAATTTTGATCAAAGCTT  >  USA300TCH1516_ALE/307994‑308174

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