Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 216,236 G→C 6.9% P191A (CCA→GCA)  rocD2_1 ← Ornithine aminotransferase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE216,2360GC 6.9% 342.2 / 8.1 119P191A (CCA→GCA) rocD2_1Ornithine aminotransferase 2
Reads supporting (aligned to +/- strand):  ref base G (48/60);  new base C (5/4);  total (54/65)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 125 total aligned reads displayed.

CATTTCCCAGTTCTACCAAGACCAACTTGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTATATGCGTCGTGGTTAGATAGTGATAATGAGC  >  USA300TCH1516_ALE/216103‑216376
                                                                                                                                     |                                                                                                                                             
cATTTCCCAGTTCTACCAAGACCAACTTGAATTTCATCTGCAATCAATAATATTTGAGGTTTATCACATAGTTGACGCCCAGCTTGAATATATCTTTTCGGTAGTCTCTTGACGCAACATTCACCTTGAATTGCTTcccca                                                                                                                                        >  2:762266/1‑138 (MQ=255)
  tttCCCAGTTCTACCAATTACAACTTGAATTTCATCAGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAaat                                                                                                                                     <  2:586320/141‑1 (MQ=255)
  tttCCCAGTTCTACCAAGACCAACTTGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAaat                                                                                                                                     <  1:898200/141‑1 (MQ=255)
   ttCCCAGTTTTACTAAAACTAACTTGAATTTCATCTGCAATCAATAATCTATGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCAAAATTCACCTTGAATTGGTTCCAAaata                                                                                                                                    <  2:899292/141‑1 (MQ=255)
    tCCCTGCTCTACCAAGACCAACTTGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAaataa                                                                                                                                   <  1:619013/141‑1 (MQ=255)
    tCCCATTTCTACCAAGACCAACTTGAATTTCATCTGCAATCAATATTATTTGATGTTTATCACATAGTTGACGCACATCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAaataa                                                                                                                                   <  1:913076/141‑1 (MQ=255)
    tCCCAGTTCTACCAAGACCAACTTGAATTTCATCTGCAATCAATGATATTTGATGTTTATCACATAGTTGTCGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCGAaataa                                                                                                                                   <  2:180773/141‑1 (MQ=255)
     cccAGTTCTACCAAGACCAACTTGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGTTATATTAACGCCACCTTCACCTTGAATTGTTTCCAAaataat                                                                                                                                  >  2:1025819/1‑141 (MQ=255)
     cccAGTTCTACCAAGACCAACTTGAATTTCATCTGCAAGCAATAACATTTGATGCTTATCACATGGTTGACCCACAGCTTGAGTATATCCTTTCGGTGCTATATTCACGCCCCATTCACCTCGCATACGTTACAAcataaa                                                                                                                                  >  2:450713/1‑140 (MQ=255)
       cAGTTCTACCAAGACCAACTTGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGCGGTATATTAACGCCACCTTCACCTTGACTTCCTTCCAATATAAttt                                                                                                                                >  2:937193/1‑140 (MQ=255)
          ttCTACCAAGACCAACTTGAATTTCATCTGCAGTGTATCATACTGGATGTTTATCACATAGTTGACGCACAGCTAGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATtgctg                                                                                                                             <  1:33280/141‑1 (MQ=255)
           tCTACCAAGACCAACTTGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGCATATTAACGCCACCTTCACCTTGAATTGGTCCCAAAATAATggctgc                                                                                                                            >  1:490261/1‑141 (MQ=255)
             caCCAAGACCAACTTAAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGctgctg                                                                                                                          <  2:1009338/140‑1 (MQ=255)
                cAAGACCAACTTGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTAt                                                                                                                       >  1:949812/1‑141 (MQ=255)
                cAAGACCAACTTGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTAt                                                                                                                       >  1:408886/1‑141 (MQ=255)
                         cTTGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGaaa                                                                                                              >  1:565412/1‑141 (MQ=255)
                          ttGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAAt                                                                                                             <  1:118951/141‑1 (MQ=255)
                           tGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATa                                                                                                            <  1:426225/141‑1 (MQ=255)
                           tGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATa                                                                                                            <  1:426117/141‑1 (MQ=255)
                           tGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTACACAATACTTTCTGCTGTATTCGGTGAAATa                                                                                                            <  1:1029635/141‑1 (MQ=255)
                             aaTTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAACATAATTGCTGCTGTATTCGGTGAAATAGc                                                                                                          >  1:265812/1‑141 (MQ=255)
                             aaTTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGTTGAAATAGc                                                                                                          >  1:950934/1‑141 (MQ=255)
                               tttCAGCTGCAATTAATAATATTTGTTGTTGAGCAGATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCtt                                                                                                        <  2:768534/141‑1 (MQ=255)
                                 tCATCTGTAATCATTAATTTTTGTTGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGTTGCTGTATTCGGTGAAATAGCTtgt                                                                                                      <  2:1024123/141‑1 (MQ=255)
                                 tCACCTGCAATCAATAAAAGTTCCTGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTCTATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTtgt                                                                                                      <  1:355570/141‑1 (MQ=255)
                                   aTCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGCTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTgtgt                                                                                                    >  1:409302/1‑141 (MQ=255)
                                   aTCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGCTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTgtgt                                                                                                    >  1:637476/1‑141 (MQ=255)
                                   aTCTGCAATCAATAATATTATATGTTTATCACATAGTTCACGCACAGCTTGAATATATGCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTgtgt                                                                                                    <  2:621897/141‑1 (MQ=255)
                                     cTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTa                                                                                                  >  1:863561/1‑141 (MQ=255)
                                     cTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTa                                                                                                  >  1:266941/1‑141 (MQ=255)
                                      tgacatcaatAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGATTGTGTTaa                                                                                                 <  2:687755/137‑1 (MQ=255)
                                      tGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACTGCTTGAATATATCCATTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAACATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTaa                                                                                                 <  2:73245/141‑1 (MQ=255)
                                       gCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAAt                                                                                                >  2:136574/1‑141 (MQ=255)
                                        caatcaatAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAAtt                                                                                               >  1:732794/1‑141 (MQ=255)
                                          atcactAAATTTTAATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCTGTGAAATAGCTTGTGTTAATTGt                                                                                             <  2:880141/141‑1 (MQ=255)
                                             aataatATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCa                                                                                          >  1:279775/1‑141 (MQ=255)
                                                  tatTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGtc                                                                                     >  2:55377/1‑141 (MQ=255)
                                                     ttGATGTTGCTCACATAGTGATAGCACAGCTTGAATATATCCTTTCGGCGGTATATTAACGCCACCTTCACCTTGACTTGGTTCCAAAATAATCGCTGCTGTATTCAATGAAATAGCTTGTGTTAATTGTTCAATGTCTcc                                                                                  <  2:60688/141‑1 (MQ=255)
                                                      tGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAACAGATTGGGTTAATTGTTCAATGTCTCCa                                                                                 >  1:435498/1‑141 (MQ=255)
                                                      tGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATCGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTccc                                                                                 >  2:367052/1‑140 (MQ=255)
                                                      tGATGTTTAGCACACAGTCGACGCACAGCTCGAATATATCCTTTCGGTGGTATATTAACGCCACCCTCACCTTGAATTGGGTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCa                                                                                 <  1:595225/141‑1 (MQ=255)
                                                             tATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTa                                                                          >  2:281668/1‑141 (MQ=255)
                                                             tATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTa                                                                          >  2:332751/1‑141 (MQ=255)
                                                             tATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTAAATCTCTCGAAAATCTa                                                                          >  2:711616/1‑141 (MQ=255)
                                                             tATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTATTAATTGTTCAATGTCTCCAAActctt                                                                          >  2:104771/1‑140 (MQ=255)
                                                              aTCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAAGAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAACATATAc                                                                         >  2:890982/1‑141 (MQ=255)
                                                                cacaTAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACtg                                                                       <  1:99409/141‑1 (MQ=255)
                                                                 acacAGTTGACGCACAGCTTGAATATATCCGTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGAATGAGTTAAGTTTTCAATGTCTCCAAAATCTACtgt                                                                      <  1:43312/141‑1 (MQ=255)
                                                                 acaTAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACtgt                                                                      <  1:133999/141‑1 (MQ=255)
                                                                   agaGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACtgttg                                                                    <  2:781803/139‑1 (MQ=255)
                                                                   aTAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAGATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACtgttg                                                                    <  1:382970/141‑1 (MQ=255)
                                                                      gTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAg                                                                 <  1:754463/141‑1 (MQ=255)
                                                                       ttGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGTTGCTGTAATCGGACAAACACCCGTTGTTAATTGTCCAATGTCTCCAAAACCTACcgtggtcgc                                                                >  1:286141/1‑138 (MQ=255)
                                                                           cccccaGCTTGAATATATCCTTTCGGTCGTATATTAACTCCACCTCCACCGTGAATTGGTTCCAAAACAATTGCTGCGGTATTCGGTGAAATAGCTGGTGTTAATTGTCCAATGTCTCCAAAATCTACTGTTGTAGTGCCt                                                            <  1:56103/137‑1 (MQ=255)
                                                                           cGCACAGCTTGAAGATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCt                                                            <  1:792259/141‑1 (MQ=255)
                                                                           aaaacaTCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCt                                                            <  1:1010059/138‑1 (MQ=255)
                                                                             ctcatgttGAATTTCTCCTTTCGGCTGTCTATTAACGCCACCTTCACCTCGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATATCTTGTGTTAATTGTTCAATGTCTCCATAATCTACTGTTGTAGTGCCTTg                                                          <  2:554880/135‑1 (MQ=255)
                                                                             cacaGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTGCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTg                                                          <  1:590378/141‑1 (MQ=255)
                                                                               caGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGaa                                                        >  2:816884/1‑141 (MQ=255)
                                                                                aGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAg                                                       >  2:919351/1‑141 (MQ=255)
                                                                                aGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAg                                                       >  2:185590/1‑141 (MQ=255)
                                                                                     gAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAggg                                                  >  1:951289/1‑141 (MQ=255)
                                                                                     gAATATATCCTTTCAGTGGTATATTAACGCCACCTTACCCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGTAAGAGCTAGTGTTAATTGTTCAATGGCTCAAAAATTTACTGTTGTAGTGCCTTGAAGTAggg                                                  <  2:735049/141‑1 (MQ=255)
                                                                                     cactctctCCTTTCGTTGGTATATTATGGCCGCGTTAACCTTGAATTGCTTCTAAAATAATTGCTTCTTTATTCGGTGAAATAGCTTGTGTTAATTGTTAAATGGCTCCAAAATCCACTGTTGTAGTGCCTTGAAGTAggg                                                  <  2:689288/134‑1 (MQ=255)
                                                                                       atatatCGTTGCGGGTGGATATTAAGGTCACCTTCGCCTTGAATTGGTTCCAAAATAAGTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTACAAAATCTACTGTTGTAGTGCCTTGAAGTAggggg                                                <  2:517026/141‑1 (MQ=255)
                                                                                       atatatCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAggggg                                                <  1:739067/141‑1 (MQ=255)
                                                                                           atcgcttCGTTGTTATATTAACGCCACCTTCACCTTGAATTTGTTCCAAAATAATTGATGCTGTATTCGGTGAAATAGCTTATGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGaa                                            <  2:643283/136‑1 (MQ=255)
                                                                                           atCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGaa                                            >  2:391874/1‑141 (MQ=255)
                                                                                           atCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTAATGTTGTAGTGCCTTGAAGTAGGGGCTGaa                                            >  2:862832/1‑141 (MQ=255)
                                                                                              cTTTCGGTTGTACAGTAACGCCACCTTCTCCTGTAATTGGTTGCAAACTTACTTATATATGTTGCATTAAAATCGCCTGGGTTTATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATc                                         <  2:467266/141‑1 (MQ=255)
                                                                                              cTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATc                                         >  1:55442/1‑141 (MQ=255)
                                                                                              cTTTCGGTGGTATATTAACGCAACATTCACCTTGAATTGGTTCCAAAATAATTGCTGATGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGaaacc                                         <  2:768950/141‑3 (MQ=255)
                                                                                                ttCGGTGGTATATTAACGCCACCTTCACCTTGAATTGCTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCt                                       <  2:637476/141‑1 (MQ=255)
                                                                                                ttCGGTGGTATATTAACGCCACCTTCACCTTGAATTGCTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCt                                       <  2:409302/141‑1 (MQ=255)
                                                                                                tgcGGTGGTATATTAACGCCACCTTCACCTTGAATTGCTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAACTCTACTGTTGTAGTGTCTTGAAGTAGGGGGTGAAATCCt                                       <  2:19909/139‑1 (MQ=255)
                                                                                                  cagtggtATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGc                                     <  1:687732/139‑1 (MQ=255)
                                                                                                    gtggtATATTAACGCCACCTTCACCTTGAATTTGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTCGTGCCTTGAAGTAGGGGGTGAAATCCTGCtt                                   >  2:702054/1‑141 (MQ=255)
                                                                                                     tggtATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCttt                                  >  1:908012/1‑141 (MQ=255)
                                                                                                        tataTTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATTTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTtata                               <  2:794636/141‑1 (MQ=255)
                                                                                                        tataTTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTtata                               <  1:730311/141‑1 (MQ=255)
                                                                                                        tataTTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTtata                               <  1:136574/141‑1 (MQ=255)
                                                                                                            ttACCGCTACCTTTACGTTGAATTGGTTCTGAAAGAATTGCGTCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCGCCAAACTCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTATATGcg                           <  2:9092/141‑1 (MQ=255)
                                                                                                            ttAACGCCACCTTCACCTTGAATTGGTCCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTATATGcg                           <  2:490261/141‑1 (MQ=255)
                                                                                                              aaCGCCACCTTCATCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTATATGcgtc                         <  1:74782/141‑1 (MQ=255)
                                                                                                                   caccttcaccttGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTATATGCGTCGTGGt                    >  1:493684/1‑141 (MQ=255)
                                                                                                                   caccttcaccttGAATTGGTTACGAAAGGGTTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTATATGCGTggtggt                    >  1:494170/1‑141 (MQ=255)
                                                                                                                    accttcaccttGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTATATGCGTCGTGGtt                   >  2:518182/1‑141 (MQ=255)
                                                                                                                     ccttcaccttGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTATATGCGTCGTGGTTa                  >  1:114263/1‑141 (MQ=255)
                                                                                                                       ttcaccttGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTATATGCGTCGTGGTTAGa                >  1:733354/1‑141 (MQ=255)
                                                                                                                        tcaccttGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGTTTTATATTCGTTGTGGTTAGAt               >  1:497558/1‑141 (MQ=255)
                                                                                                                        tcaccttGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTCTATGCGTCGTGGTTAGAt               >  1:185686/1‑141 (MQ=255)
                                                                                                                        tcaccttGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTATATGCGTCGTGGTTAgct               >  1:569375/1‑139 (MQ=255)
                                                                                                                          tccttTAATGGTTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTATATGCGTCGTGGTTAGATAg             <  1:133248/140‑1 (MQ=255)
                                                                                                                          accttGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTATATGCGTCGTGGTTAGATAg             <  1:417156/141‑1 (MQ=255)
                                                                                                                            cttGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTATATGCGTCGTGGTTAGATAGTg           <  2:443009/141‑1 (MQ=255)
                                                                                                                            cttGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTATATGCGTCGTGGTTAGATAGTg           <  2:258132/141‑1 (MQ=255)
                                                                                                                            cgtGAATTGGGTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTATATGCGTCGTGGTTAGATAGTg           <  1:419881/139‑1 (MQ=255)
                                                                                                                                    ggTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTATATGCGTCGTGGTTAGATAGTGATAATGAg   <  2:493684/141‑1 (MQ=255)
                                                                                                                                     gTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTATATGCGTCGTGGTTCGATAGTGATAATGAGc  >  2:484653/1‑141 (MQ=255)
                                                                                                                                     gTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTAAATGCGTCGTGGTTAGATAGTGGTAATGAGc  >  1:356342/1‑141 (MQ=255)
                                                                                                                                     |                                                                                                                                             
CATTTCCCAGTTCTACCAAGACCAACTTGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTATATGCGTCGTGGTTAGATAGTGATAATGAGC  >  USA300TCH1516_ALE/216103‑216376

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