Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,963,502 C→T 6.7% T453T (ACG→ACA hemY ← Protoporphyrinogen oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,963,5020CT 6.7% 417.1 / 11.6 134T453T (ACG→ACAhemYProtoporphyrinogen oxidase
Reads supporting (aligned to +/- strand):  ref base C (57/68);  new base T (3/6);  total (60/74)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 135 total aligned reads displayed.

CACATTACTTCCACCATTATATAACTTAAAAATGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGAATATGACCGACATGGTACTGTGGCATACTTTTCGGCA  >  USA300TCH1516_ALE/1963362‑1963641
                                                                                                                                            |                                                                                                                                           
cacaTTACTTCCACCATTATATAACTTAAAAATGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGc                                                                                                                                             >  1:965586/1‑141 (MQ=255)
 acaTTACTTCCACCATTATATAACTTAAAAATGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCg                                                                                                                                            >  1:1030471/1‑141 (MQ=255)
     tACTTCCACCATTATATAACTTAAAAATGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAAt                                                                                                                                        <  1:1033746/141‑1 (MQ=255)
     tACTTCCACCATTATATAACTTAAAAATGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAAt                                                                                                                                        <  2:95219/141‑1 (MQ=255)
      aCTTCCACCATTATATAACTTAAAAATGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATa                                                                                                                                       >  1:793478/1‑141 (MQ=255)
       cTTCCACCATTATATAACTTAAAAATGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATAc                                                                                                                                      <  2:1049067/141‑1 (MQ=255)
       cTTCCACCATTATATAACTTAAAAATGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATAc                                                                                                                                      <  2:774262/141‑1 (MQ=255)
         tCCACCATTATATAACTTAAAAATGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACaa                                                                                                                                    <  2:799811/141‑1 (MQ=255)
         tCCACCATTATATAACTTAAAAATGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACaa                                                                                                                                    <  2:855103/141‑1 (MQ=255)
          ccaccaTTATATAACTTAAAAATGACGATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAAt                                                                                                                                   <  2:861242/141‑1 (MQ=255)
             ccaTTATATAACTTAAAAATGACTATATTTCATCAAACATTATCTAAAGGGGTGGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGAACCTTAACCATGCGAAAAACAACCAg                                                                                                                                >  3:115675/1‑141 (MQ=255)
             ccaTTATATAACTTAAAAATGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGCTTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAg                                                                                                                                >  1:1024986/1‑141 (MQ=255)
                 tatataACTTAAAAATGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAg                                                                                                                            <  1:930744/141‑1 (MQ=255)
                  atataACTTAAAAATGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGt                                                                                                                           >  1:89479/1‑141 (MQ=255)
                   tataACTTAAAAATGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTc                                                                                                                          >  2:742718/1‑141 (MQ=255)
                   tataACTTAAAAATGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGAGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTc                                                                                                                          <  2:130641/141‑1 (MQ=255)
                     taACTTAAAAATGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCa                                                                                                                        >  2:818022/1‑141 (MQ=255)
                       aCTTAAAAATGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAAc                                                                                                                      >  1:60930/1‑141 (MQ=255)
                        cTTAAAAATGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAAcc                                                                                                                     <  2:291390/141‑1 (MQ=255)
                         ttAAAAATGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCg                                                                                                                    >  1:248576/1‑141 (MQ=255)
                            aaaaTGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGcct                                                                                                                 >  1:773851/1‑139 (MQ=255)
                            aaaaTGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCtt                                                                                                                 <  1:404770/141‑1 (MQ=255)
                            aaaaTGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCtt                                                                                                                 <  2:748203/141‑1 (MQ=255)
                            aaaaTGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCtt                                                                                                                 <  3:449985/141‑1 (MQ=255)
                              aaTGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCa                                                                                                               <  2:584825/141‑1 (MQ=255)
                              aaTGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCa                                                                                                               >  1:383050/1‑141 (MQ=255)
                                  aCTATATTTCATCAAACATTATCTCAAGGAGTAGCACCTACACAAACACAATCAAAAAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTACTCCACTCCGGTAGTGCAACCGCTTCAAAAg                                                                                                           >  2:539786/1‑141 (MQ=255)
                                  aCTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAGCAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAg                                                                                                           >  2:285504/1‑141 (MQ=255)
                                  aCTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAg                                                                                                           >  2:370530/1‑141 (MQ=255)
                                  aCTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAg                                                                                                           >  1:724810/1‑141 (MQ=255)
                                   cTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGa                                                                                                          <  2:266386/141‑1 (MQ=255)
                                      tatTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCGTCAGCAGCAGCTTTACCTTGTGTAATACAGTCAGGTAGTCCAACCGCTTCAAAAGATGc                                                                                                       <  2:1058923/141‑1 (MQ=255)
                                           catcaAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAg                                                                                                  <  1:499259/141‑1 (MQ=255)
                                             tcaAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGtt                                                                                                >  1:831473/1‑141 (MQ=255)
                                             tcaAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGAACCAGtt                                                                                                <  1:898824/141‑1 (MQ=255)
                                                aaCATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTAct                                                                                             <  2:806641/141‑1 (MQ=255)
                                                 aCATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTActc                                                                                            <  3:180421/141‑1 (MQ=255)
                                                 aCATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTActc                                                                                            >  2:906384/1‑141 (MQ=255)
                                                  cATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGTGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTActct                                                                                           >  1:396547/1‑141 (MQ=255)
                                                  cATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTActct                                                                                           <  1:473257/141‑1 (MQ=255)
                                                  cATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTActct                                                                                           <  1:765808/141‑1 (MQ=255)
                                                   aTTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTa                                                                                          >  1:627387/1‑141 (MQ=255)
                                                     tATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGTGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAg                                                                                        <  2:40228/141‑1 (MQ=255)
                                                      aTCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGt                                                                                       <  1:663561/141‑1 (MQ=255)
                                                        cTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCg                                                                                     <  1:294073/141‑1 (MQ=255)
                                                         tAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGt                                                                                    >  2:683652/1‑141 (MQ=255)
                                                            aGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGa                                                                                 <  4:376403/141‑1 (MQ=255)
                                                            aGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGa                                                                                 >  1:744252/1‑141 (MQ=255)
                                                            aGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGATTCGAAAGATGCACCAGTTACTCTAAGTCGTGGa                                                                                 <  2:976922/141‑1 (MQ=255)
                                                             ggCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGat                                                                                >  2:1060556/1‑141 (MQ=255)
                                                             ggCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGat                                                                                >  4:212391/1‑141 (MQ=255)
                                                               cGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTCCGCCGTGGatat                                                                              >  2:250654/1‑141 (MQ=255)
                                                                 tCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATAtgt                                                                            >  2:62890/1‑141 (MQ=255)
                                                                  cGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATAtgtt                                                                           >  1:906533/1‑141 (MQ=255)
                                                                  cGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATAtgtt                                                                           >  3:422777/1‑141 (MQ=255)
                                                                   gCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATAtgttt                                                                          >  1:125694/1‑141 (MQ=255)
                                                                      ccTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATAtgtttgtt                                                                       >  1:840193/1‑141 (MQ=255)
                                                                      ccTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATAtgtttgtt                                                                       >  1:513805/1‑141 (MQ=255)
                                                                         acaccaacaccaTCCGACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTGTACCTTGCGTAATACAATCGGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTaa                                                                    <  1:242344/141‑1 (MQ=255)
                                                                         acaccaacaccaTCCAACAATTAACTTACAACGCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGACGCACCAGTTACTCTAAGTCGTGGATATgttggtttca                                                                    >  2:227093/1‑139 (MQ=255)
                                                                          caccaaaaccaTCCAACAATTAACTTACAACTCAGCGATTACTTCTTCAGCAACCACTTTACCTAGCGTAATACAATCAGGTAGTCCAACAGCTTCAAACGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAat                                                                   >  2:544604/1‑141 (MQ=255)
                                                                           accaacaccaTCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAATAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAata                                                                  <  2:747022/141‑1 (MQ=255)
                                                                             caacaccaTCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATAtg                                                                <  3:106365/141‑1 (MQ=255)
                                                                               acaccaTCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTGATAtgtg                                                              <  1:860978/141‑1 (MQ=255)
                                                                                caccaTCCAACAATTAACTTACAACTCTGCGATTAGTTCTTCTGCAGCAACTTTACCTTGCGTAATAAAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGc                                                             >  1:934266/1‑141 (MQ=255)
                                                                                  ccaTCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCtt                                                           >  1:608241/1‑141 (MQ=255)
                                                                                   caTCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTg                                                          <  2:202190/141‑1 (MQ=255)
                                                                                        aacaaTTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCt                                                     <  4:293304/141‑1 (MQ=255)
                                                                                        aacaaTTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCt                                                     <  2:947954/141‑1 (MQ=255)
                                                                                          caaTTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATCCAATCAGGTAGTCCAACCGCTTCAAACGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGt                                                   <  1:648553/141‑1 (MQ=255)
                                                                                               aaCTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAAt                                              >  3:403970/1‑141 (MQ=255)
                                                                                                     cAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATttgttg                                        <  3:277291/141‑1 (MQ=255)
                                                                                                     cAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATttgttg                                        <  2:510368/141‑1 (MQ=255)
                                                                                                     cAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATttgttg                                        <  2:506305/141‑1 (MQ=255)
                                                                                                     cAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATttgttg                                        <  2:418583/141‑1 (MQ=255)
                                                                                                       aCTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGaa                                      <  1:329911/141‑1 (MQ=255)
                                                                                                        ctctGCGATTACTTCTTCAGCAGCAACTTTAGCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGAat                                     <  2:487364/141‑1 (MQ=255)
                                                                                                           tGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGAATATg                                  <  1:430959/141‑1 (MQ=255)
                                                                                                           tGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGAATATg                                  <  4:298004/141‑1 (MQ=255)
                                                                                                            gCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGAATATGa                                 >  4:124709/1‑141 (MQ=255)
                                                                                                            gCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGAATATGa                                 <  3:45033/141‑1 (MQ=255)
                                                                                                              gATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGAATATGAcc                               <  1:353738/141‑1 (MQ=255)
                                                                                                                 tACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGAATATGACCGAc                            <  1:226980/141‑1 (MQ=255)
                                                                                                                   cttcttCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGAATATGACCGACAt                          >  3:98812/1‑141 (MQ=255)
                                                                                                                             agcaACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGAATATGACCGACATGGTACTGTgg                <  4:4415/141‑1 (MQ=255)
                                                                                                                              gcaACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGAATATGACCGACATGGTACTGTGGc               <  2:613990/141‑1 (MQ=255)
                                                                                                                                aaCTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGAATATGACCGACATGGTACTGTGGCAt             <  3:285274/141‑1 (MQ=255)
                                                                                                                                 aCTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGAATATGACCGACATGGTACTGTGGCATa            >  3:80612/1‑141 (MQ=255)
                                                                                                                                   tttACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGAATATGACCGACATGGTACTGTGGCATACt          >  3:382634/1‑141 (MQ=255)
                                                                                                                                   tttACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGAATATGACCGACATGGTACTGTGGCATACt          >  1:56899/1‑141 (MQ=255)
                                                                                                                                    ttACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGGTGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGAATATGACCGACATGGTACTGTGGCATACtt         <  1:762683/141‑1 (MQ=255)
                                                                                                                                     tACCTTGCGTAATACAATCAGGTAGTTCAACCGCTTCAAAATATGCACCAGTTACTCGAATGGTTTTTTATTTGTTTTTAATTTGGGTTTGATTCTGTCTAATTTGTTGTATATGATTGACATTGTACTTTGGCATAtttt        >  3:224220/1‑141 (MQ=255)
                                                                                                                                     tACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTGTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGAATATGACCGACATGGTACTGTGGCATACttt        <  1:1053680/141‑1 (MQ=255)
                                                                                                                                       ccTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTACTATGTTCTTGAATCTGTCTAATTTGTTGAATATGACCGACATGGTACTGTGGCATGCTTTTc      >  2:60327/1‑141 (MQ=255)
                                                                                                                                       ccTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGAATATGACCGACATGGTACTGTGGCATACttttt      >  4:449574/1‑140 (MQ=255)
                                                                                                                                       ccTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGAATATGACCGACATGGTACTGTGGCATACTTTTc      <  1:746884/141‑1 (MQ=255)
                                                                                                                                        cTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGAATATGACCGACATGGTACTGTGGCATACTTTTCg     >  1:997740/1‑141 (MQ=255)
                                                                                                                                         ttGCTNAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGAATATGACCGACATGGTACTGTGGCATACTTTTCgg    <  3:141385/141‑1 (MQ=255)
                                                                                                                                          tGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGAATATGACCGACATGGTACTGTGGCATACTTTTCGGc   >  3:319927/1‑141 (MQ=255)
                                                                                                                                           gCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGAATATGACCGACATGGTACTGTGGCATACTTTTCGGCa  <  1:170652/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                           
CACATTACTTCCACCATTATATAACTTAAAAATGACTATATTTCATCAAACATTATCTAAAGGCGTCGCACCTACACCAACACCATCCAACAATTAACTTACAACTCTGCGATTACTTCTTCAGCAGCAACTTTACCTTGCGTAATACAATCAGGTAGTCCAACCGCTTCAAAAGATGCACCAGTTACTCTAAGTCGTGGATATGTTTGTTTAATATGTGCTTGAATCTGTCTAATTTGTTGAATATGACCGACATGGTACTGTGGCATACTTTTCGGCA  >  USA300TCH1516_ALE/1963362‑1963641

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