Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 162,116 G→A 8.5% A205V (GCA→GTA)  phnE_2 ← Phosphate‑import permease protein PhnE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE162,1160GA 8.5% 222.4 / 2.9 80A205V (GCA→GTA) phnE_2Phosphate‑import permease protein PhnE
Reads supporting (aligned to +/- strand):  ref base G (35/36);  new base A (4/4);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CGACGATAACTAAAGTAAATAAAATAATCGTTGCTGTTTTTGGATATTGAAATAAACCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCTTCCATTGGACCTGGATCTATCGTTTCCAATGATTCATATAATA  >  USA300TCH1516_ALE/161976‑162254
                                                                                                                                            |                                                                                                                                          
cgacgaTAACTAAAGTAAATAAAATAATCGTTGCTGTTTTTGGATATTGAAATAAACCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATg                                                                                                                                            >  1:895476/1‑141 (MQ=255)
cgacgaTAACTAAAGTAAATAAAATAATCGTTGCTGTTTTTGGATATTGAAATAAACCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATg                                                                                                                                            >  1:839668/1‑141 (MQ=255)
cgacgaTAACTAAAGTAAATAAAATAATCGTTGCTGTTTTTGGATATTGAAATAAACCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCGGAAGCCCGTATATTTACTTCAACTg                                                                                                                                            >  1:764318/1‑141 (MQ=255)
 aacgaTAACTAAAGTAAATAAAATAATCGTTGCTGTTTTTGGATATTGAAATAAACCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGc                                                                                                                                           <  1:165965/140‑1 (MQ=255)
     aTAACTAAAGTAAATAAAATAATCGTTGCTGTTTTTGGATATTGAAATAAACCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCatat                                                                                                                                       <  1:554762/141‑1 (MQ=255)
       aaCTAAAGTAAATAAAATAATCGTTGCTGTTTTTGGATATTGAAATAAACCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATataa                                                                                                                                     <  1:174508/141‑1 (MQ=255)
       aaCTAAAGTAAATAAAATAACCGTTGCTGTTTTTGGATATTGAAATAAACCACGGCTGTAAGCATAAAACAATCCAATACCGCGAGCGCCGACTAAACCAAGAACACCTTAAGCACGTCTATTTACTTCCAATCCATctaa                                                                                                                                     >  2:402275/1‑141 (MQ=255)
         cTAAAGTAAATAAAATAATCGTTGCTGTTTTTGGATATTGAAATAAACCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCCCGTAACTTAACTTCAAATGAATatcata                                                                                                                                   >  2:939863/1‑141 (MQ=255)
            aaGTAAATAAAATAATCGTTGCTGTTTTTGGATATTGAAATAAACCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCGAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGt                                                                                                                                <  2:327462/141‑1 (MQ=255)
                   taaaagaaTTGTTGCTGTTTTTGGAAACTGAAATAATCCAAATGTTGTATCACAAAACAATCCACTACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACAt                                                                                                                         <  1:75913/141‑1 (MQ=255)
                     aaataaTCGTTGCTGTTTTTGGATATTGAAATAAACCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATAAAATGCGTATGACATaa                                                                                                                       >  2:838570/1‑141 (MQ=255)
                      aataatCGTTGCTGTTTTTGGATATTGAAATAAACCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATaaa                                                                                                                      <  1:700471/141‑1 (MQ=255)
                       ataatCGTTGCTGTTTTTGGATATTGAAATAAACCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTCTATTTACTTCAAATGCATATAATACGTATGACATAAAt                                                                                                                     <  1:681420/141‑1 (MQ=255)
                       ataatCGTTGCTGTTTTTGGATATTGAAATAAACCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATATACTTCAAATGCATATAATACGTATGACATAAAt                                                                                                                     >  2:846098/1‑141 (MQ=255)
                       ataatCGTTGCTGTTTATGGATATTGAAATAAACCAAGAGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATCTTTACTTCAAATGCATATAATACGTATGACATAAAt                                                                                                                     >  1:590835/1‑141 (MQ=255)
                        taatCGTTGCTGTTTTTGGATATTGAAATAAACCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATACATATAATACGTATTACATAAATg                                                                                                                    >  1:279878/1‑141 (MQ=255)
                        taatCGTTGCTGTTTTTGGATATTGAAATAAACCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATACATATAATACGTATGACATAAAta                                                                                                                    >  1:279670/1‑140 (MQ=255)
                        taatCGTTGCTGTTTTTGGATATTGAAATAAACCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATACATATAATACGTATGACATAAATg                                                                                                                    >  1:760347/1‑141 (MQ=255)
                          atCGTTGCTGTTTTTGGATATTGAAATAAACCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCGTATATTACGTATGACATAAATGAc                                                                                                                  >  2:589301/1‑141 (MQ=255)
                               tGCTGTTTTTGGATATTGAAATAAACCAAGTGTTTGATCATAAAACAATCCAATCCCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGatat                                                                                                             <  2:596641/141‑1 (MQ=255)
                               tGCTGTTTTTGGATATTGAAATAAACCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGatat                                                                                                             <  1:11661/141‑1 (MQ=255)
                               tGCTGTTTTTGGAAATTGAAATAAACCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGatat                                                                                                             <  1:1013781/141‑1 (MQ=255)
                                    tttgtGGATATTGAAATAAACCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCtt                                                                                                        <  1:355717/141‑1 (MQ=255)
                                        tGGATATTGAAATAAACCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTgg                                                                                                    <  2:930891/141‑1 (MQ=255)
                                            tatTGAAATAAACCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATACATATAATACGTATGACATAAATGGCGATATGGCTTGTgggaaa                                                                                                >  2:631301/1‑137 (MQ=255)
                                                  aaTAAACCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCg                                                                                          <  2:853367/141‑1 (MQ=255)
                                                   aTAAACCGAGTGTTTGATCATAACACAATCCAACCCCGCCAGCGCCGACTATTCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCAGATAATACGTAGGACATAAATGACGATATGGCTTGTGGTACAACACCGa                                                                                         <  2:876819/141‑1 (MQ=255)
                                                   aTAAACCGAGTGTTTGATCATAAAACAATGCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACCTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGa                                                                                         <  1:638911/141‑1 (MQ=255)
                                                     aagccAAGTGTATGATCATAACACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGaaa                                                                                       <  1:1032566/138‑1 (MQ=255)
                                                      aaCCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGaaaa                                                                                      >  1:332435/1‑141 (MQ=255)
                                                      aaCCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGaaaa                                                                                      >  1:912990/1‑141 (MQ=255)
                                                       aCCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAAc                                                                                     >  1:688213/1‑141 (MQ=255)
                                                       aCCAAGTGTTTGATCATAAAACAAACCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGAAGATATAGCTTGGGGTACAACACCAaaaaa                                                                                     >  2:751461/1‑140 (MQ=255)
                                                              gtTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCat                                                                              >  1:997579/1‑141 (MQ=255)
                                                               tttGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCatt                                                                             <  2:839668/141‑1 (MQ=255)
                                                                    tCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCattttatt                                                                        >  2:98575/1‑141 (MQ=255)
                                                                     cATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCattttattt                                                                       >  1:1009597/1‑141 (MQ=255)
                                                                      aTAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCattttatttt                                                                      >  1:800951/1‑141 (MQ=255)
                                                                             aaTCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGGTATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAgcc                                                               >  1:83947/1‑140 (MQ=255)
                                                                               tCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAACTGCATATAATACGTATGCACTCTCTGCGGATAGGGCTTGTGGTACAACACCGAAAACAATCACTTTTATTTTATTAGCGcc                                                             >  1:244255/1‑141 (MQ=255)
                                                                                ccAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCa                                                            <  1:1050166/141‑1 (MQ=255)
                                                                                 cAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCaa                                                           >  1:386854/1‑141 (MQ=255)
                                                                                   attccGCCAGCGCCGACGAATCCAAGCACAGCTGAAGCACGTATATTTACTGCAAATGCAGATAATACGTAGGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACa                                                         <  1:1011642/138‑1 (MQ=255)
                                                                                   aTACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACa                                                         <  1:155867/141‑1 (MQ=255)
                                                                                   aTACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACa                                                         <  2:627652/141‑1 (MQ=255)
                                                                                   aTACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACa                                                         <  2:57073/141‑1 (MQ=255)
                                                                                   aTACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACa                                                         <  2:556025/141‑1 (MQ=255)
                                                                                           aGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTTTGGTACAACACCGAAAACAATCCATTTTATTTTATAAGCGCCAACAGCCgccgt                                                 >  2:424115/1‑137 (MQ=255)
                                                                                           aGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCAt                                                 <  2:875329/141‑1 (MQ=255)
                                                                                            gcgcCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCAtt                                                >  2:360644/1‑141 (MQ=255)
                                                                                            gcgcCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCAtt                                                <  1:875589/141‑1 (MQ=255)
                                                                                            gcgcCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCAtt                                                <  1:875481/141‑1 (MQ=255)
                                                                                             cgcCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTg                                               >  1:852852/1‑141 (MQ=255)
                                                                                             cgaccaCTAATCCAAGAACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTg                                               <  1:306704/136‑1 (MQ=255)
                                                                                               ccGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCt                                             <  1:679209/141‑1 (MQ=255)
                                                                                               ccGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGAATGTAATGCTGATGACATCAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCt                                             <  1:994133/141‑1 (MQ=255)
                                                                                                    tAAGCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACCATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCTTCCAt                                        <  2:25976/141‑1 (MQ=255)
                                                                                                    aaactcAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCTTCCAt                                        <  2:789135/136‑1 (MQ=255)
                                                                                                    aaaTCCAAGCACAGCTGAAGCTCGTATATTTACTTCAAATACATATCATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCTTCCAt                                        <  2:628824/140‑1 (MQ=255)
                                                                                                      gtCCAAGCACAGCTGAAGCACGTATATTTACTTCAAGTGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCTTCCATTg                                      <  2:869979/140‑1 (MQ=255)
                                                                                                      aTCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGGCGACAGCCGTCATTGCTTACATTg                                      >  1:753241/1‑141 (MQ=255)
                                                                                                       tCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATAACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCTTCCCTTgg                                     >  2:307942/1‑141 (MQ=255)
                                                                                                        ccAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCTTCCATTGGa                                    >  1:390324/1‑141 (MQ=255)
                                                                                                             cacaGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCTTCCATTGGACCTgg                               >  2:392497/1‑141 (MQ=255)
                                                                                                             cacaGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTACTTAGCGCCAACAGCCCTCAATGCTTCCATTGTACCTgg                               >  2:466529/1‑141 (MQ=255)
                                                                                                               caGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCTTCCATTGGACCTGGAt                             >  2:379836/1‑141 (MQ=255)
                                                                                                               caGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGAGCCAACAGACGTCATTGCTTCCATTGGAACTGGAt                             >  1:741157/1‑141 (MQ=255)
                                                                                                                      aGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCTTCCATTGGACCTGGATCTATCGt                      <  1:878311/141‑1 (MQ=255)
                                                                                                                      aGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCTTCCATTGGACCTGGATCTATCGt                      >  1:767972/1‑141 (MQ=255)
                                                                                                                      aGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCTTCCATTGGACCTGGATCTATCGt                      <  1:87417/141‑1 (MQ=255)
                                                                                                                      aGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCTTCCATTGGACCTGGATCTATCGt                      <  1:877482/141‑1 (MQ=255)
                                                                                                                      aGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTCTATTAGCGCCAACAGCCGTCATTGCTTCCATAGGACCTGGATCTATCGt                      >  1:1011696/1‑141 (MQ=255)
                                                                                                                         aagTATATATACTACATACGCACACAATACGTATCACATAAATGACGATATGTCTTGTTGTACAACACCGAATACAATCCATTTTATTTAATTAGCGCCAACAGCCGTCATTGCTTCCATTGGACCTGGATCTATCGTTTc                   <  2:822011/139‑1 (MQ=255)
                                                                                                                          cGTATATTTACTTCAAATACATATAATACGTATGACATAAATGACGATATGGCTTGTGGTGCAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCTTCCATTGGACCTGGATCTATCGTTTcc                  <  2:760347/141‑1 (MQ=255)
                                                                                                                          cGTATATTTACTTCAAATACATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACATTCCATGTTATTTTATTAGCGCCAACAGCCGTCATTGCTTCCATTGGACCTGGATCTATCGTTTcc                  <  2:279670/141‑1 (MQ=255)
                                                                                                                          cGTATATTTACTTCAAATACATATAATACGTATGACATAAATGACCATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCTTCCATTGGACCTGGATCTATCGTTTcc                  <  2:279878/141‑1 (MQ=255)
                                                                                                                           gTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCTTCCATTGGACCTGGATCTATCGTTTCCa                 >  2:662074/1‑141 (MQ=255)
                                                                                                                                tttACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCTTCCATTGGACCTGGATCTGTCGTTTCCAATggt            >  1:383208/1‑139 (MQ=255)
                                                                                                                                tttACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCTTCCATTGGACCTGGATCTATCGTTTCCACTGAt            >  1:383187/1‑141 (MQ=255)
                                                                                                                                     ttCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCTTCCATTGGACCTGGATCTATCGTTTCCAATGATTCata       >  1:9726/1‑141 (MQ=255)
                                                                                                                                      tCAAATGCATATAATAGGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCTTCCATTGGACCTGGATCTATCGTTTCCAATGATTCatat      >  1:631271/1‑141 (MQ=255)
                                                                                                                                         aaTGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCTTCCATTGGACCTGGATCTATCGTTTCCAATGATTCATataat   >  1:336695/1‑141 (MQ=255)
                                                                                                                                          aTGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCTTCCATTGGACCTGGATCTATCGTTTCCAATGATTCATataata  >  1:820937/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                          
CGACGATAACTAAAGTAAATAAAATAATCGTTGCTGTTTTTGGATATTGAAATAAACCAAGTGTTTGATCATAAAACAATCCAATACCGCCAGCGCCGACTAATCCAAGCACAGCTGAAGCACGTATATTTACTTCAAATGCATATAATACGTATGACATAAATGACGATATGGCTTGTGGTACAACACCGAAAACAATCCATTTTATTTTATTAGCGCCAACAGCCGTCATTGCTTCCATTGGACCTGGATCTATCGTTTCCAATGATTCATATAATA  >  USA300TCH1516_ALE/161976‑162254

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: