Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,963,768 C→G 24.0% V365L (GTA→CTA)  hemY ← Protoporphyrinogen oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,963,7680CG24.0% 270.5 / 91.1 137V365L (GTA→CTA) hemYProtoporphyrinogen oxidase
Reads supporting (aligned to +/- strand):  ref base C (57/47);  new base G (18/15);  total (75/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

Only 100 of 139 total aligned reads displayed.

CATACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCTTTT  >  USA300TCH1516_ALE/1963628‑1963907
                                                                                                                                            |                                                                                                                                           
cATACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAc                                                                                                                                             <  1:970016/141‑1 (MQ=255)
 aTACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCACTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACa                                                                                                                                            <  1:823700/141‑1 (MQ=255)
    cTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATAc                                                                                                                                         <  1:1063157/141‑1 (MQ=255)
      tttCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGc                                                                                                                                       >  2:1136823/1‑141 (MQ=255)
      tttCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGc                                                                                                                                       >  2:143284/1‑141 (MQ=255)
      tttCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACATGGCTTACCTACATCCGc                                                                                                                                       >  2:52043/1‑141 (MQ=255)
        tCGTCAAACCATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCAc                                                                                                                                     <  1:305513/141‑1 (MQ=255)
        tCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACGACAGTATCACCTGGTTTACCTACATACGCAc                                                                                                                                     >  2:640393/1‑141 (MQ=255)
        tCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCCCAGTATCACCCGGTTTACCTCCATACGCAc                                                                                                                                     >  2:471016/1‑141 (MQ=255)
        tCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCAc                                                                                                                                     <  2:712900/141‑1 (MQ=255)
        tCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGCTTACCTAGATACCCAc                                                                                                                                     >  1:596038/1‑141 (MQ=255)
        tCGGCAAACGAGTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCCCAGAATCACCTGGTTTACCTACATACGCAc                                                                                                                                     >  2:937911/1‑141 (MQ=255)
         cGGCAAACGATTGACACTTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACg                                                                                                                                    <  2:1016295/141‑1 (MQ=255)
         cGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATActctcg                                                                                                                                    <  2:1037036/141‑6 (MQ=255)
         cGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACg                                                                                                                                    >  2:523643/1‑141 (MQ=255)
         cGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACg                                                                                                                                    <  2:1118846/141‑1 (MQ=255)
         cGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACg                                                                                                                                    <  1:183022/141‑1 (MQ=255)
          ggCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGa                                                                                                                                   <  2:895934/141‑1 (MQ=255)
          ggCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTAAAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGa                                                                                                                                   <  2:576020/141‑1 (MQ=255)
            cAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAAt                                                                                                                                 <  2:1191145/141‑1 (MQ=255)
             aaaCGATTGACAATGGTAAATTCGGGATCACCTTTAAATGTCATCATTTGACTTAAATCGCTACGTACAATCGATACTAATTCATTATCTGTATGATCATGAGCCACAGTATCACCTGGTTTACCTACCTACCACCGAATc                                                                                                                                >  2:640848/1‑141 (MQ=255)
               aCGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCCGAAATCaa                                                                                                                              >  2:99301/1‑141 (MQ=255)
                 gATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCaaaa                                                                                                                            >  1:851531/1‑141 (MQ=255)
                  aTTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAac                                                                                                                           >  1:302014/1‑141 (MQ=255)
                  aTTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGACCGAATCAAAac                                                                                                                           >  2:241371/1‑141 (MQ=255)
                   ttGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACACCAGTATCACCTGGTTTAACTACATACGCCCGAATCAAAacc                                                                                                                          >  2:472516/1‑141 (MQ=255)
                    tGACAATTGTAAATTCAGGATCACTTTTAAATGTCATCATTTGACTTAAAGCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAacct                                                                                                                         <  2:522160/141‑1 (MQ=255)
                         aTTGTAAATTCAGGATCACCTTGAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCAGAGTAGCACCTGGTTTACCTAGATACGCACGAATCAAAaccttacct                                                                                                                    <  1:525365/141‑1 (MQ=255)
                              aaaTTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCgg                                                                                                               >  1:626791/1‑141 (MQ=255)
                              aaaTTCAGGATCACCTTTAAATGTCAGCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCgg                                                                                                               <  2:171079/141‑1 (MQ=255)
                                aTTCAGGATCACCTTAAAATGTCATCATTTGACTTAAATCTCTACGCACAATCGATACTACAACCTTATCTGTATGACCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAACCCTTACCTTCCGgtg                                                                                                             <  1:457711/141‑1 (MQ=255)
                                 ttCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGgtgt                                                                                                            <  2:483468/141‑1 (MQ=255)
                                        tCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAAt                                                                                                     >  2:1097797/1‑141 (MQ=255)
                                        tCACCTGTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAAt                                                                                                     >  2:1066070/1‑141 (MQ=255)
                                           ccTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGc                                                                                                  <  1:259960/141‑1 (MQ=255)
                                             tttAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCa                                                                                                <  2:204175/141‑1 (MQ=255)
                                              ttAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGccct                                                                                               >  2:52392/1‑139 (MQ=255)
                                                   tGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAAACACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTCGTAAATGGCCAtttttt                                                                                          >  2:303442/1‑141 (MQ=255)
                                                      catcatTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGa                                                                                       >  1:1163041/1‑141 (MQ=255)
                                                      catcatTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGa                                                                                       >  1:897530/1‑141 (MQ=255)
                                                       ggcatTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTAACTACATACGCACTAATCAAAACCTTCCATTCCGGATTAGTAAATGGCCATTTTTTCGAt                                                                                      >  2:190896/3‑141 (MQ=255)
                                                        tgatTTAACTTAAATCTCTAAGTACACTCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATg                                                                                     <  2:95289/139‑1 (MQ=255)
                                                        tcatTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGAtt                                                                                     >  2:413723/1‑140 (MQ=255)
                                                              gACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAg                                                                               <  1:948307/141‑1 (MQ=255)
                                                               aCTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAACCCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGt                                                                              <  2:986364/141‑1 (MQ=255)
                                                               aCTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGt                                                                              <  1:143284/141‑1 (MQ=255)
                                                                gttAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGtc                                                                             >  2:735502/2‑140 (MQ=255)
                                                                cTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTa                                                                             >  2:179620/1‑141 (MQ=255)
                                                                cTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTCCCGGTGTTGTAAATGGCCATTTTTTCGTTGTCCAAGtc                                                                             >  2:1001827/1‑140 (MQ=255)
                                                                cTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGCTGTCCAAGTa                                                                             >  1:249577/1‑141 (MQ=255)
                                                                attAAATCACTACGTACAATAGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTCGTAAATGGCCATTTTTTCGATGTCCAAGTa                                                                             >  2:403817/2‑141 (MQ=255)
                                                                 ttAAATCTCTACGTACACTCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTAc                                                                            <  1:153898/141‑1 (MQ=255)
                                                                 ttAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTAc                                                                            <  1:146320/141‑1 (MQ=255)
                                                                  tAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACa                                                                           <  1:1134503/141‑1 (MQ=255)
                                                                   aaaTCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGACTTCCGAGTAAAt                                                                          >  2:997075/1‑141 (MQ=255)
                                                                    aatccgtACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATg                                                                         <  2:156413/135‑1 (MQ=255)
                                                                      tctcTACGTACAATCTATACTAACTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCg                                                                       <  1:411515/141‑1 (MQ=255)
                                                                      tctcTACGTACAATCGATACTAATTCATTAGCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCg                                                                       <  1:472516/141‑1 (MQ=255)
                                                                      tccttgcGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGATTCAAATCCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCg                                                                       <  1:410937/135‑1 (MQ=255)
                                                                        tctACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACGTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACAGTCGGt                                                                     >  1:36154/1‑141 (MQ=255)
                                                                        tctACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCAACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGt                                                                     <  2:891696/141‑1 (MQ=255)
                                                                           aCGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAAt                                                                  >  1:1037095/1‑141 (MQ=255)
                                                                               aCAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTTCAAGTGCATGCGGTaaggact                                                              >  1:1169669/1‑136 (MQ=255)
                                                                                   tCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTAt                                                          <  1:502303/141‑1 (MQ=255)
                                                                                     gATACTCATCCACTACCTGTATTATCATCCACCACAGTCTAACATGGTTTCCCTGCATTCGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTCTCGATGTCCAAGTACATGCGGTAATGTCTGTATCa                                                        <  1:161420/141‑1 (MQ=255)
                                                                                     gATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGGAATGTCTGTATCa                                                        >  2:349085/1‑141 (MQ=255)
                                                                                            aTTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTCCCTACATACGCACGAATCAAAACCTTACCTTCCGGCGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTACCACTCGTtc                                                 >  2:1170876/1‑141 (MQ=255)
                                                                                              tCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGGTTCACGCGTtctc                                               >  1:447575/1‑141 (MQ=255)
                                                                                                ctcaTCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGc                                             <  1:555182/138‑1 (MQ=255)
                                                                                                     cTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTa                                        >  1:787171/1‑141 (MQ=255)
                                                                                                     cTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTa                                        >  1:641077/1‑141 (MQ=255)
                                                                                                      tGTATGATCATCAACCCCAATGTCACCTGGTTTACCTACATACGCACGAATGAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCCATTTCCCAGGACATGCGGTAATGGCGGTATCACTCGTTCTCGcgactag                                       >  2:850595/1‑135 (MQ=255)
                                                                                                          tGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCGCGAATAAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAAGGTCTATGTCACTCGTTCTCGCACTTACGAAg                                   >  1:1127592/1‑141 (MQ=255)
                                                                                                          tGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAg                                   >  1:1127896/1‑141 (MQ=255)
                                                                                                          tGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCACAATTACGAAg                                   >  1:820781/1‑141 (MQ=255)
                                                                                                          tGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTGCTCGCAATTACGAAg                                   >  1:339023/1‑141 (MQ=255)
                                                                                                          tGATCATCAACCAACGTATCACCTGGTTTACCTACATACGCACAAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCGATTACGAAg                                   >  1:235838/1‑141 (MQ=255)
                                                                                                           gATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTGCTTTCTTTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGc                                  >  1:513658/1‑141 (MQ=255)
                                                                                                           gATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGc                                  <  2:1139370/141‑1 (MQ=255)
                                                                                                            atcatcAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGAGGTCCAAGTACATGCGGTAATGTCGGTATCACTCGTTCTCGCAATCACGAAGcc                                 >  2:460768/1‑141 (MQ=255)
                                                                                                             tcatcaACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGcaa                                >  1:925960/1‑139 (MQ=255)
                                                                                                              catcaACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAg                               <  1:478324/141‑1 (MQ=255)
                                                                                                                   aCCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATTGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTAGCAATTACGTAGCCAGTAccc                          >  1:30212/1‑140 (MQ=255)
                                                                                                                    ccACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCCGTGCcaa                         >  2:41881/1‑140 (MQ=255)
                                                                                                                      acaGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTCGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCCCGCAATTACGAAGCCAGTACcatca                       >  2:449359/1‑141 (MQ=255)
                                                                                                                      acaGTACCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcatca                       <  2:258718/141‑1 (MQ=255)
                                                                                                                       caGTATCACCTGGTTTACCTGGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcatcat                      <  2:168432/141‑1 (MQ=255)
                                                                                                                       caGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcatcat                      >  2:861935/1‑141 (MQ=255)
                                                                                                                       caGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTTTTTTTTTTTTTACATGCAGTAATGTCTGTATCACTCGTTCTCGCAATTGCGAAGCCAGTACcatcat                      >  2:282046/1‑141 (MQ=255)
                                                                                                                       caGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCCGTATCACTCTTTCTCGCAATTACGAAGCCAGTACcatcat                      <  1:1066070/141‑1 (MQ=255)
                                                                                                                         gTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCACTCACGAAGCCAGTACCATCATgg                    >  1:268353/1‑141 (MQ=255)
                                                                                                                           ttCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGt                  <  2:941867/140‑1 (MQ=255)
                                                                                                                               ccTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTAttt              >  1:1199567/1‑141 (MQ=255)
                                                                                                                                cTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGACCATTTGTTCGATGTCCAAGTACGTGCGGTAATGTCTGTATCACTCGTTCTCGCAATTCAGAAGCCAGTACCATCATGGGCAgttt             >  2:328331/1‑141 (MQ=255)
                                                                                                                                 tGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTc            <  2:36154/141‑1 (MQ=255)
                                                                                                                                  cgcttACCTACATACACACCAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGACGTCCAAGTACATGCGGTAATGTCCGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCa           <  2:268729/138‑1 (MQ=255)
                                                                                                                                     ttACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTCTCACTCGTTCGCGCAATTACGAAGCCAGTACCATCATGGGTCGTTTCaagg        >  1:682544/1‑139 (MQ=255)
                                                                                                                                        ccTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCt     >  1:445929/1‑141 (MQ=255)
                                                                                                                                         cTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCtt    >  1:303244/1‑141 (MQ=255)
                                                                                                                                           aCATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGTTGTCGAAGTGCATGCGCGAATGTCTGTATCACTGGTTCTCGCAATTACGCAGCCAGTACCATCATGGGTATTTTCGCTGTCtttt  >  2:960416/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                           
CATACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCTTTT  >  USA300TCH1516_ALE/1963628‑1963907

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: