Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 1,971,067 A→C K9Q (AAG→CAG)  cspC → cold‑shock protein CspC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH11,971,0670AC98.7% 1520.8 / 4.2 393K9Q (AAG→CAG) cspCcold‑shock protein CspC
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base C (220/166);  total (223/170)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01

Only 100 of 393 total aligned reads displayed.

AAACACGTTGACGTTGAATACCGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTAAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGAAATTCAGGACGGCC  >  DH1/1970920‑1971214
                                                                                                                                                   |                                                                                                                                                   
aaaCACGTTGACGTTGAATACCGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGt                                                                                                                                                  <  2:355075/151‑1 (MQ=255)
aaaCACGTTGACGTTGAATACCGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGt                                                                                                                                                  <  7:1404997/151‑1 (MQ=255)
 aaCACGTTGACGTTGAATACCGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTg                                                                                                                                                 <  4:250419/151‑1 (MQ=255)
     cgttgacgttgaATACCGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTc                                                                                                                                             <  1:1362296/151‑1 (MQ=255)
        tgacgttgaATACCGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATTCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGAAAAGATTAAAGGTCAGGTTCAGTGGTTCAAc                                                                                                                                          >  2:1455363/1‑151 (MQ=255)
        tgacgttgaATACCGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCaa                                                                                                                                           <  2:430090/150‑1 (MQ=255)
         gacgttgaATACCGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTGAGTTTTATCTTTGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGACAGGTTCAGTGGTTCAAcc                                                                                                                                         >  1:1083587/1‑150 (MQ=255)
         gacgttgaATACCGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACg                                                                                                                                         <  6:1266510/151‑1 (MQ=255)
           cgttgaATACCGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACg                                                                                                                                         >  4:1395449/1‑149 (MQ=255)
            gttgaATACCGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGt                                                                                                                                      >  3:1168730/1‑151 (MQ=255)
                aaTACCGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTa                                                                                                                                   >  8:172500/1‑150 (MQ=255)
                  tACCGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCGGTGGTTCAACGAGTCTaa                                                                                                                                  >  8:31978/1‑149 (MQ=255)
                     cGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGtt                                                                                                                             >  5:302602/1‑151 (MQ=255)
                       cTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGtttt                                                                                                                           >  8:624646/1‑151 (MQ=255)
                         tctcTGCCGAGCATTTTATTGGTGCGTCTTGCAGTTTTGTGTCAGTTTTATCTTTGTTATGCCTGCGGGCGATGAAAACAATCTAAGGATTTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTg                                                                                                                          <  3:219571/150‑1 (MQ=255)
                           tctGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCt                                                                                                                       >  4:722009/1‑151 (MQ=255)
                            ctGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCtt                                                                                                                      >  3:327485/1‑151 (MQ=255)
                              gCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTc                                                                                                                     >  5:674548/1‑150 (MQ=255)
                               ccGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCAt                                                                                                                   <  8:902869/151‑1 (MQ=255)
                                 gatccTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAATAAAACAATCTAAGGAATTTTTCAAATGTCAAAGATTAAATGTCAGGTTGAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTa                                                                                                                 <  3:55898/148‑1 (MQ=255)
                                       tataTTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCgg                                                                                                           >  6:1115776/1‑151 (MQ=255)
                                         tatTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCt                                                                                                         >  4:457861/1‑151 (MQ=255)
                                         tatTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCt                                                                                                         >  7:634963/1‑151 (MQ=255)
                                          atTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCt                                                                                                         <  8:373241/150‑1 (MQ=255)
                                          atTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTg                                                                                                        >  6:674610/1‑151 (MQ=255)
                                           ttGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGTAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTGTAAAGGTTTTGGCTTCATTACTCCGGCTg                                                                                                        >  7:825325/1‑150 (MQ=255)
                                             ggTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGAt                                                                                                      <  6:225356/150‑1 (MQ=255)
                                             ggTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGAt                                                                                                      <  2:252023/150‑1 (MQ=255)
                                             ggTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATg                                                                                                     >  8:416178/1‑151 (MQ=255)
                                              gTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATgg                                                                                                    >  3:442343/1‑151 (MQ=255)
                                                    cATGCAGTAATGTGTCTGTTTTATCTATGTTATGCCTGCGGGCGAAGATATCAATCTAAGGAATTTTTCAATTGGCTAAGATTTAAGGTCTGGTTGAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTGATTACTCCGGCTGATGgcagca                                                                                               <  6:471202/150‑1 (MQ=255)
                                                    cATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGgcagc                                                                                                <  1:1011901/149‑1 (MQ=255)
                                                     aTGCATTATTGTGTGAGTTGTATCTATGTTATGGCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATTGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCaaa                                                                                             <  7:958664/151‑1 (MQ=255)
                                                       gCAGTAATGTGTCAGTTTTATCTATGTTATGCGTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGa                                                                                           <  3:193183/151‑1 (MQ=255)
                                                         aGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGAtg                                                                                         >  5:558189/1‑151 (MQ=255)
                                                          gTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGGTTCATTACTCCGGCTGATGGCAGCAAAGAtgt                                                                                        >  4:990494/1‑151 (MQ=255)
                                                          gTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGAtgt                                                                                        >  7:499626/1‑151 (MQ=255)
                                                            aaTGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGAtgtg                                                                                       <  1:1506123/150‑1 (MQ=255)
                                                             aTGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGtt                                                                                     >  6:78911/1‑151 (MQ=255)
                                                             aTGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGtt                                                                                     <  6:132560/151‑1 (MQ=255)
                                                              tgtgTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTc                                                                                    >  5:1484250/1‑151 (MQ=255)
                                                               gtgtCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCg                                                                                   <  6:538920/151‑1 (MQ=255)
                                                                    aGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTacac                                                                              <  3:92421/151‑1 (MQ=255)
                                                                    aGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTacac                                                                              >  3:1049238/1‑151 (MQ=255)
                                                                    aGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTac                                                                                >  2:1542214/1‑149 (MQ=255)
                                                                        ttttCTATGTTATGCCTGCGGGAGAATAAAACAATCTAAGGAATTTTTCAAATGGGAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGATTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTtc                                                                           <  2:1367896/147‑1 (MQ=255)
                                                                         tATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCGAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTtct                                                                          <  3:755474/150‑1 (MQ=255)
                                                                           tCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCg                                                                       >  5:917733/1‑151 (MQ=255)
                                                                              aTGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCg                                                                       <  1:129530/148‑1 (MQ=255)
                                                                              aTGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTa                                                                    <  5:1276709/151‑1 (MQ=255)
                                                                               tGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTAt                                                                   >  7:1247214/1‑151 (MQ=255)
                                                                               tGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTAt                                                                   >  3:1262683/1‑151 (MQ=255)
                                                                                gTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGGCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTAt                                                                   <  6:1179242/150‑1 (MQ=255)
                                                                                    tGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAgg                                                              >  7:1084468/1‑151 (MQ=255)
                                                                                     gCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAgg                                                              <  2:690367/150‑1 (MQ=255)
                                                                                        tGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCACTGGTTCAACGAGTCTAAAGGTTTTGGCATCATTACTCCGGCTGATGGCAGCAAAGATTTGTTCCTACACTTCTCCGCTAACCatcgta                                                           >  1:568048/1‑145 (MQ=255)
                                                                                         gCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGATCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTaa                                                          <  1:434537/150‑1 (MQ=255)
                                                                                          cGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATg                                                        <  5:1389378/151‑1 (MQ=255)
                                                                                           gggCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATgg                                                       >  5:832133/1‑151 (MQ=255)
                                                                                            ggCGAAGAAAACAATCTAATGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTGAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGc                                                      <  4:1475533/151‑1 (MQ=255)
                                                                                               gaagaaAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTc                                                   >  8:1350898/1‑151 (MQ=255)
                                                                                               gaagaaAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTc                                                   >  7:279116/1‑151 (MQ=255)
                                                                                               gaagaaAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTc                                                   >  6:33087/1‑151 (MQ=255)
                                                                                               gaagaaAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGATTCATTACTCCGGCTCATGCCAGCAAAGATGTGTTCGTACACTTCTCCGCAATCCAGGGTAATGGCtt                                                    >  8:1429767/1‑150 (MQ=255)
                                                                                                 agaaAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCa                                                  >  5:1136234/1‑150 (MQ=255)
                                                                                                  gaaAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCaa                                                 >  6:1018120/1‑150 (MQ=255)
                                                                                                   aaaaCAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCaaa                                                >  6:1478355/1‑150 (MQ=255)
                                                                                                   aaaaCAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTAAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCaaa                                                >  3:839035/1‑150 (MQ=255)
                                                                                                      actagcTAAGGAATTTATCAAATGGCAAAGCTTCAAGGTCAGGTTCAGTGGTTCAACGAGGCTACAGCTTTTGGCTTCATTCCTCCGGCTGATGGCAGCAAAGATGTGTTCGTAGACTTCACCGCTCTCCTGGGTAATGGCTTCAAAActg                                            <  8:1073901/146‑2 (MQ=255)
                                                                                                      aCAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCCGGTTCAGTGGTTCAACGCGTCTAAAGGTTTTGGCTGCAGTACTCCGGCTGCGGGCGGCAACGATGTGTTCGTAGCCTGCTCCGCTATCCAGGGTACTGGCGTCCACActc                                            >  8:151537/1‑151 (MQ=255)
                                                                                                      aCAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAActc                                            <  1:732834/151‑1 (MQ=255)
                                                                                                        aaTCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTg                                          >  2:569651/1‑151 (MQ=255)
                                                                                                          tCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGc                                        >  8:1285309/1‑151 (MQ=255)
                                                                                                           cTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCt                                       >  1:1216108/1‑151 (MQ=255)
                                                                                                             aaGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTg                                      >  7:347394/1‑150 (MQ=255)
                                                                                                                gAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAgg                                  <  7:312480/151‑1 (MQ=255)
                                                                                                                 aaTTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGt                                 >  7:1517540/1‑151 (MQ=255)
                                                                                                                  aTTTTTCAAATGGCAAAGATTAAAGGGCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCGGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTc                                <  7:565243/151‑1 (MQ=255)
                                                                                                                    ttttCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCa                               <  3:1480073/150‑1 (MQ=255)
                                                                                                                    ttttCAAATGGCAAAGATTAAAGGTCAGGTTAAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTc                                >  1:171796/1‑149 (MQ=255)
                                                                                                                       tCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAAc                           >  5:892400/1‑151 (MQ=255)
                                                                                                                        cAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACg                          >  5:4270/1‑151 (MQ=255)
                                                                                                                        cAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACg                          >  6:84475/1‑151 (MQ=255)
                                                                                                                         aaaTGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGt                         >  8:515583/1‑151 (MQ=255)
                                                                                                                         aaaTGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGt                         <  4:1545649/151‑1 (MQ=255)
                                                                                                                         aaaTGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGt                         >  4:593167/1‑151 (MQ=255)
                                                                                                                            tGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTg                       >  5:233120/1‑150 (MQ=255)
                                                                                                                              gCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGt                    <  5:1319917/151‑1 (MQ=255)
                                                                                                                                aaaGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTc                  >  4:1125486/1‑151 (MQ=255)
                                                                                                                                aaaGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTc                  <  2:1267563/151‑1 (MQ=255)
                                                                                                                                 aaGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCg                 >  4:1243924/1‑151 (MQ=255)
                                                                                                                                   gATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGaa               <  8:764337/151‑1 (MQ=255)
                                                                                                                                     ttAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGAAAt             >  2:821783/1‑151 (MQ=255)
                                                                                                                                       aaaGGGCAGGTGCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGAAAtt            <  3:1428958/150‑1 (MQ=255)
                                                                                                                                         aGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGAAATTc           >  7:131261/1‑149 (MQ=255)
                                                                                                                                          ggTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACTCTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGAAATTCAgg        <  1:268313/151‑1 (MQ=255)
                                                                                                                                            tCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGAAATTCAGGAc      >  4:82106/1‑151 (MQ=255)
                                                                                                                                              aGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGAAATTCAGGACgg    <  1:999530/151‑1 (MQ=255)
                                                                                                                                               ggTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGAAATTCAGGACGGc   >  8:1214993/1‑151 (MQ=255)
                                                                                                                                                gTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGAAATTCAGGACGGcc  >  1:822944/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAACACGTTGACGTTGAATACCGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTAAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGAAATTCAGGACGGCC  >  DH1/1970920‑1971214

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