Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,607,818 A→T 10.2% H381Q (CAT→CAA)  sad ← succinate semialdehyde dehydrogenase, NAD(P)+‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,607,8180AT10.2% 557.9 / 45.1 187H381Q (CAT→CAA) sadsuccinate semialdehyde dehydrogenase, NAD(P)+‑dependent
Reads supporting (aligned to +/- strand):  ref base A (97/71);  new base T (4/15);  total (101/86)
Fisher's exact test for biased strand distribution p-value = 3.00e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 189 total aligned reads displayed.

GGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAGCCATCTTTTCCCCGCCCAGTAACAAA  >  CP009273/1607671‑1607965
                                                                                                                                                   |                                                                                                                                                   
ggCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCa                                                                                                                                                     >  1:796253/1‑148 (MQ=255)
   cACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGt                                                                                                                                                  <  2:3154301/148‑1 (MQ=255)
    aCTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTAACTATCATTAGCCAGTTCCAGTGCATGtt                                                                                                                                                 >  2:1836226/1‑148 (MQ=255)
     cTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCTTGTTc                                                                                                                                                <  1:3380/148‑1 (MQ=255)
     cTCGCGCGTCGCTGACACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTc                                                                                                                                                >  1:3099359/1‑148 (MQ=255)
       cgcgcgTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTg                                                                                                                                              >  1:181615/1‑148 (MQ=255)
           cgTCGCTGGCACAATAACCCTTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATTTTCTGCATc                                                                                                                                          >  1:210281/1‑148 (MQ=255)
              cGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCttt                                                                                                                                       >  2:2139059/1‑148 (MQ=255)
                 tGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGc                                                                                                                                    >  2:719211/1‑148 (MQ=255)
                  ggCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCa                                                                                                                                   <  1:3029401/148‑1 (MQ=255)
                    cacaATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAAt                                                                                                                                 >  2:2303584/1‑148 (MQ=255)
                     acaATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCTTGTTCTGCATCTTTCGCAATg                                                                                                                                <  1:2312963/148‑1 (MQ=255)
                      caATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCTTGTTCTGCATCTTTCGCAATgg                                                                                                                               <  2:1491464/148‑1 (MQ=255)
                         tAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTAACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGa                                                                                                                            >  2:2375941/1‑148 (MQ=255)
                          aaCCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGAt                                                                                                                           >  1:2919649/1‑148 (MQ=255)
                            ccATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTAACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTg                                                                                                                         >  2:575490/1‑148 (MQ=255)
                               ttGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTAACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCg                                                                                                                      <  1:208726/148‑1 (MQ=255)
                                   tAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAAc                                                                                                                  >  2:1080143/1‑148 (MQ=255)
                                    aaaCACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATCAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACg                                                                                                                 <  1:1339491/148‑1 (MQ=255)
                                      acacCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACggg                                                                                                               <  1:3215328/148‑1 (MQ=255)
                                      acacCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACggg                                                                                                               <  1:3213540/148‑1 (MQ=255)
                                      acacCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACggg                                                                                                               >  2:745409/1‑148 (MQ=255)
                                            cACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCaaa                                                                                                         >  2:89862/1‑148 (MQ=255)
                                            cACCGCATTCCAGACGAGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTAACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCaaa                                                                                                         <  2:3123439/148‑1 (MQ=255)
                                              ccGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCTTCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACa                                                                                                       <  2:20976/148‑1 (MQ=255)
                                              ccGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACa                                                                                                       >  1:69194/1‑148 (MQ=255)
                                              ccGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACa                                                                                                       >  1:69230/1‑148 (MQ=255)
                                                 cATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACAttt                                                                                                    >  2:1064963/1‑148 (MQ=255)
                                                  aTTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCTTGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATttc                                                                                                   <  2:16393/148‑1 (MQ=255)
                                                    tCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCTTGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATttctt                                                                                                 <  2:813148/148‑1 (MQ=255)
                                                     ccAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATttcttc                                                                                                <  2:1928481/148‑1 (MQ=255)
                                                      cAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTcc                                                                                               >  2:2408675/1‑148 (MQ=255)
                                                        gACGTGCAGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACAATCATTAGCCAGATCCAGTGCATGTTCTGAATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTccct                                                                                             >  2:929148/1‑147 (MQ=255)
                                                          cGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGaa                                                                                           <  1:1428033/148‑1 (MQ=255)
                                                          cGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGaa                                                                                           <  1:2055719/148‑1 (MQ=255)
                                                           gTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGaaa                                                                                          <  2:2433123/148‑1 (MQ=255)
                                                            tGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAAc                                                                                         >  2:477609/1‑148 (MQ=255)
                                                             gccgccATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAAcg                                                                                        <  1:645123/148‑1 (MQ=255)
                                                              ccgccATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAAcgc                                                                                       <  2:1893754/148‑1 (MQ=255)
                                                                 ccATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGt                                                                                    >  2:1539031/1‑148 (MQ=255)
                                                                 ccATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTAACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGt                                                                                    <  2:421006/148‑1 (MQ=255)
                                                                   aTCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCTTGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCa                                                                                  <  1:1532293/148‑1 (MQ=255)
                                                                   aTCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCTTGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCa                                                                                  <  1:1532427/148‑1 (MQ=255)
                                                                   aTCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCa                                                                                  >  1:395734/1‑148 (MQ=255)
                                                                    tctgtctgGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCAt                                                                                 >  1:456037/1‑148 (MQ=255)
                                                                     ctgtctgGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCAtt                                                                                <  1:2588385/148‑1 (MQ=255)
                                                                      tgtctgGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCAttt                                                                               <  2:1590184/148‑1 (MQ=255)
                                                                      tgtctgGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTAACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCAttt                                                                               <  2:2602237/148‑1 (MQ=255)
                                                                         ctgGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCTTGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTg                                                                            >  2:2837367/1‑148 (MQ=255)
                                                                         ctgGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTg                                                                            >  2:1379382/1‑148 (MQ=255)
                                                                          tgGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTgg                                                                           >  2:1236016/1‑148 (MQ=255)
                                                                          tgGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTgg                                                                           >  2:2886164/1‑148 (MQ=255)
                                                                            gCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTTCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTgggg                                                                         >  1:2814428/1‑148 (MQ=255)
                                                                            gCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTgggg                                                                         <  1:2867083/148‑1 (MQ=255)
                                                                               tgtgTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTATGGGGTaa                                                                      <  1:944565/148‑1 (MQ=255)
                                                                                  gtTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACAt                                                                   >  1:2922529/1‑148 (MQ=255)
                                                                                   tttCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACAtt                                                                  >  1:100690/1‑148 (MQ=255)
                                                                                   tttCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACAtt                                                                  >  1:100160/1‑148 (MQ=255)
                                                                                    ttCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTc                                                                 >  2:1173452/1‑148 (MQ=255)
                                                                                         cAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTAACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAg                                                            >  1:1388423/1‑148 (MQ=255)
                                                                                          aGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGa                                                           <  2:194442/148‑1 (MQ=255)
                                                                                          aGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGa                                                           <  2:194579/148‑1 (MQ=255)
                                                                                            tGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAAc                                                         <  1:345492/148‑1 (MQ=255)
                                                                                               tAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGt                                                      <  1:332598/148‑1 (MQ=255)
                                                                                                 aaaaTGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTg                                                    >  1:552618/1‑148 (MQ=255)
                                                                                                    aTGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCTAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCg                                                 >  2:1647341/1‑148 (MQ=255)
                                                                                                     tGGTCGCTGAAAGGCCGAACTACCTCTCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACAGTCGCCAGAACCGTTGGCgg                                                <  2:3285826/148‑1 (MQ=255)
                                                                                                       gTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGAt                                              <  1:641632/148‑1 (MQ=255)
                                                                                                        tCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCTTGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGAta                                             >  2:13262/1‑148 (MQ=255)
                                                                                                        tCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGAta                                             <  1:1308222/148‑1 (MQ=255)
                                                                                                        tCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGAGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGAta                                             <  2:64457/148‑1 (MQ=255)
                                                                                                         cGCTGAAAGGCCGAACTAACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGAtag                                            >  1:1641355/1‑148 (MQ=255)
                                                                                                          gCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGAtagt                                           <  1:2761731/148‑1 (MQ=255)
                                                                                                             gAAAGGCCGAACTAACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATagtagt                                        <  1:2915791/148‑1 (MQ=255)
                                                                                                              aaaGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGtt                                       >  2:588253/1‑148 (MQ=255)
                                                                                                              aaaGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGtt                                       >  2:588242/1‑148 (MQ=255)
                                                                                                                 ggCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTAcc                                    >  2:1238093/1‑148 (MQ=255)
                                                                                                                    cGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGc                                 <  2:3099359/148‑1 (MQ=255)
                                                                                                                      aaCTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGccc                               >  1:2631057/1‑148 (MQ=255)
                                                                                                                       aCTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGcccc                              <  2:2975058/148‑1 (MQ=255)
                                                                                                                          cACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAGc                           >  1:1774420/1‑148 (MQ=255)
                                                                                                                          cACTATCATTAGCCAGTTCCAGGGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAGc                           <  2:6537/148‑1 (MQ=255)
                                                                                                                           aCTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCTGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAGcc                          >  1:3236009/1‑148 (MQ=255)
                                                                                                                            cTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAGCCa                         >  2:1863546/1‑148 (MQ=255)
                                                                                                                             tATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAGCCAt                        >  1:1799535/1‑148 (MQ=255)
                                                                                                                                cATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAGCCATCtt                     >  1:844505/1‑148 (MQ=255)
                                                                                                                                 aTTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAGCCATCttt                    >  2:803425/1‑148 (MQ=255)
                                                                                                                                     gCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAGCCATCTTTTccc                >  1:119357/1‑148 (MQ=255)
                                                                                                                                     gCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAGCCATCTTTTccc                <  2:2739629/148‑1 (MQ=255)
                                                                                                                                         gTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAGCCATCTTTTCCCCGcc            <  1:3061839/148‑1 (MQ=255)
                                                                                                                                          ttCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAGCCATCTTTTCCCCGccc           <  2:393297/148‑1 (MQ=255)
                                                                                                                                          ttCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAGCCATCTTTTCCCCGccc           >  1:1626684/1‑148 (MQ=255)
                                                                                                                                          ttCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAGCCATCTTTTCCCCGccc           >  1:2694159/1‑148 (MQ=255)
                                                                                                                                          ttCCAGTGCATGTTCTGCATCTTTCGAAATGGTGATTGCCGCAAAGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACAGTTGGCGGATAGTAGTTACCTGCCCCAGCCATCTTTTCCACGccc           >  1:1626680/1‑148 (MQ=255)
                                                                                                                                           tCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAGCCATCTTTTCCCCGCCCa          <  1:748158/148‑1 (MQ=255)
                                                                                                                                            ccAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAGCCATCTTTTCCCCGCCCAg         >  1:1686008/1‑148 (MQ=255)
                                                                                                                                              aGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAGCCATCTTTTCCCCGCCCAGTa       >  1:3265673/1‑148 (MQ=255)
                                                                                                                                              aGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAGCCATCTTTTCCCCGCCCAGTa       <  2:2790551/148‑1 (MQ=255)
                                                                                                                                               gTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAGCCATCTTTTCCCCGCCCAGTaa      >  2:3272162/1‑148 (MQ=255)
                                                                                                                                                   aTGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAGCCATCTTTTCCCCGCCCAGTAACaaa  <  1:199316/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAGCCATCTTTTCCCCGCCCAGTAACAAA  >  CP009273/1607671‑1607965

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