Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,607,793 C→A 5.7% E390* (GAG→TAG)  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,7930CA 5.7% 797.9 / 15.7 208E390* (GAG→TAG) sadsuccinate semialdehyde dehydrogenase, NAD(P)+‑dependent
Reads supporting (aligned to +/- strand):  ref base C (98/98);  new base A (7/5);  total (105/103)
Fisher's exact test for biased strand distribution p-value = 7.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 211 total aligned reads displayed.

AGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAGC  >  CP009273/1607646‑1607940
                                                                                                                                                   |                                                                                                                                                   
aGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTc                                                                                                                                                     >  1:3113156/1‑148 (MQ=255)
 gCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCa                                                                                                                                                    <  1:3267068/148‑1 (MQ=255)
 gCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCa                                                                                                                                                    <  1:1660500/148‑1 (MQ=255)
  ccACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCAc                                                                                                                                                   >  2:1041848/1‑148 (MQ=255)
   cACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACt                                                                                                                                                  <  2:3093420/148‑1 (MQ=255)
        tttttCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCAt                                                                                                                                             <  1:2698077/148‑1 (MQ=255)
         ttttCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGGCGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCAtt                                                                                                                                            >  1:2251190/1‑148 (MQ=255)
         ttttCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCAtt                                                                                                                                            >  1:2075314/1‑148 (MQ=255)
              aCGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCa                                                                                                                                       >  1:2596540/1‑148 (MQ=255)
               cGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAg                                                                                                                                      <  2:165178/148‑1 (MQ=255)
                 ccacccAAGGCCACTCGCGCGTCGCTGGCACAATCACCAGTGAGAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTAACTATCATTAGCCAGtt                                                                                                                                    <  2:3003088/148‑1 (MQ=255)
                  caccaAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTc                                                                                                                                   >  2:2462711/1‑148 (MQ=255)
                     caAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAg                                                                                                                                <  2:1893766/148‑1 (MQ=255)
                     caAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAg                                                                                                                                >  1:2498608/1‑148 (MQ=255)
                        aGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGc                                                                                                                             >  1:1017886/1‑148 (MQ=255)
                        aGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGc                                                                                                                             >  1:3393889/1‑148 (MQ=255)
                        aGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGc                                                                                                                             >  1:1491235/1‑148 (MQ=255)
                             aCTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGtt                                                                                                                        <  1:2656184/148‑1 (MQ=255)
                              cTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTc                                                                                                                       <  2:3633451/148‑1 (MQ=255)
                               tCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCt                                                                                                                      <  1:3427159/148‑1 (MQ=255)
                                  cgcgTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCa                                                                                                                   >  1:1727147/1‑148 (MQ=255)
                                   gcgTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCAt                                                                                                                  <  2:2908426/148‑1 (MQ=255)
                                   gcgTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCAt                                                                                                                  <  2:3591714/148‑1 (MQ=255)
                                       cGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCttt                                                                                                              >  2:1183113/1‑148 (MQ=255)
                                       cGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCttt                                                                                                              >  1:3203001/1‑148 (MQ=255)
                                         cTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTAACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCg                                                                                                            >  1:898604/1‑148 (MQ=255)
                                            gCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCaa                                                                                                         >  2:3256377/1‑148 (MQ=255)
                                              acaATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATg                                                                                                       <  1:621231/148‑1 (MQ=255)
                                                aaTAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGt                                                                                                     >  2:3673320/1‑148 (MQ=255)
                                                  tAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGa                                                                                                   <  2:3758508/148‑1 (MQ=255)
                                                  tAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGa                                                                                                   <  1:2147508/148‑1 (MQ=255)
                                                  tAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGa                                                                                                   <  2:1500502/148‑1 (MQ=255)
                                                  tAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGa                                                                                                   >  2:1253729/1‑148 (MQ=255)
                                                  tAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTAACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGa                                                                                                   >  1:1492253/1‑148 (MQ=255)
                                                       aTTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGcc                                                                                              >  2:289217/1‑148 (MQ=255)
                                                         tGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCCATGGTGATTGCCGc                                                                                            >  2:24114/1‑148 (MQ=255)
                                                          gATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCa                                                                                           <  1:2301300/148‑1 (MQ=255)
                                                          gATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCa                                                                                           <  1:319435/148‑1 (MQ=255)
                                                            tAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAAc                                                                                         >  1:1747877/1‑148 (MQ=255)
                                                             aaaCACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACg                                                                                        >  1:3653730/1‑148 (MQ=255)
                                                             aaaCACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACg                                                                                        <  1:1915073/148‑1 (MQ=255)
                                                               acacCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACggg                                                                                      >  1:2964475/1‑148 (MQ=255)
                                                                cacCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACgggg                                                                                     <  1:3499779/148‑1 (MQ=255)
                                                                cacCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACgggg                                                                                     >  1:3360097/1‑148 (MQ=255)
                                                                        cGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACAt                                                                             <  1:3786013/148‑1 (MQ=255)
                                                                        cGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACAt                                                                             >  1:833729/1‑148 (MQ=255)
                                                                          cATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGCTTGCCGCCACGGGGCCAAACAttt                                                                           >  2:2712032/1‑148 (MQ=255)
                                                                          cATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACAttt                                                                           >  2:2711372/1‑148 (MQ=255)
                                                                           aTTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATttc                                                                          <  2:1388346/148‑1 (MQ=255)
                                                                             tCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATttctt                                                                        <  2:2257143/148‑1 (MQ=255)
                                                                             tCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATttctt                                                                        >  2:1624124/1‑148 (MQ=255)
                                                                             tCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATttctt                                                                        <  2:2116790/148‑1 (MQ=255)
                                                                               cAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTcc                                                                      <  2:185475/148‑1 (MQ=255)
                                                                               cAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTcc                                                                      <  1:369552/148‑1 (MQ=255)
                                                                                 gACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCg                                                                    >  2:2440253/1‑148 (MQ=255)
                                                                                   cGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGaa                                                                  <  2:822525/148‑1 (MQ=255)
                                                                                   cGTGCCGCCATCTGTCTGGCCTGTGTTTAGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGaa                                                                  >  2:1277494/1‑148 (MQ=255)
                                                                                    gTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGaaa                                                                 <  1:988784/148‑1 (MQ=255)
                                                                                    gTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGaaa                                                                 <  2:3139828/148‑1 (MQ=255)
                                                                                      gccgccATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTAACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAAcg                                                               <  1:2073575/148‑1 (MQ=255)
                                                                                           cATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTc                                                          >  1:2137536/1‑148 (MQ=255)
                                                                                            aTCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCa                                                         <  1:3689242/148‑1 (MQ=255)
                                                                                            aTCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCa                                                         <  1:2374020/148‑1 (MQ=255)
                                                                                             tctgtctgGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCAt                                                        <  2:944906/148‑1 (MQ=255)
                                                                                                 tctgTCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGAGGTCATTTCt                                                    >  2:2573356/1‑148 (MQ=255)
                                                                                                 tctgGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCt                                                    >  2:1809893/1‑148 (MQ=255)
                                                                                                 tctgGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCt                                                    >  2:3700429/1‑148 (MQ=255)
                                                                                                    ggCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTggg                                                 <  2:3264307/148‑1 (MQ=255)
                                                                                                     gCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTgggg                                                <  2:2881620/148‑1 (MQ=255)
                                                                                                        tgtgTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTaa                                             <  2:2547303/148‑1 (MQ=255)
                                                                                                        tgtgTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTaa                                             >  2:2651584/1‑148 (MQ=255)
                                                                                                         gtgtTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAAc                                            >  1:1144356/1‑148 (MQ=255)
                                                                                                         gtgtTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAAc                                            <  1:2197484/148‑1 (MQ=255)
                                                                                                         gtgtTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAAc                                            >  1:81734/1‑148 (MQ=255)
                                                                                                           gtTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACAt                                          >  2:2817602/1‑148 (MQ=255)
                                                                                                            tttCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCACTGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTCACAtt                                         >  2:2322032/1‑148 (MQ=255)
                                                                                                               cGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGc                                      <  1:1624124/148‑1 (MQ=255)
                                                                                                                 tCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCa                                    >  2:2995762/1‑148 (MQ=255)
                                                                                                                  cAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAg                                   >  1:1291104/1‑148 (MQ=255)
                                                                                                                  cAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAg                                   >  1:2257177/1‑148 (MQ=255)
                                                                                                                   aGTGGTAAAAATGGTCGCTGAAAGGCCGAACTAACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGa                                  <  1:1886607/148‑1 (MQ=255)
                                                                                                                    gTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGaa                                 <  2:3203001/148‑1 (MQ=255)
                                                                                                                     tGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAAc                                >  2:1883893/1‑148 (MQ=255)
                                                                                                                       gTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCg                              <  1:1183113/148‑1 (MQ=255)
                                                                                                                       gTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCg                              >  1:2736021/1‑148 (MQ=255)
                                                                                                                        tAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGt                             <  2:3165113/148‑1 (MQ=255)
                                                                                                                         aaaaaTGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGtt                            <  1:2884018/148‑1 (MQ=255)
                                                                                                                          aaaaTGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTg                           >  1:304952/1‑148 (MQ=255)
                                                                                                                           aaaTGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTgg                          >  2:792411/1‑148 (MQ=255)
                                                                                                                                gTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGAt                     >  1:3563864/1‑148 (MQ=255)
                                                                                                                                   gCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGAtagt                  <  2:60467/148‑1 (MQ=255)
                                                                                                                                       aaaGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGtt              >  2:44896/1‑148 (MQ=255)
                                                                                                                                       aaaGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGtt              >  1:454592/1‑148 (MQ=255)
                                                                                                                                        aaGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTa             >  1:2181243/1‑148 (MQ=255)
                                                                                                                                        aaGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTa             >  2:2747337/1‑148 (MQ=255)
                                                                                                                                         aGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTAc            <  2:1733391/148‑1 (MQ=255)
                                                                                                                                          ggCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTAcc           <  2:3185945/148‑1 (MQ=255)
                                                                                                                                              gAACTAACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGcc       >  2:1127951/1‑148 (MQ=255)
                                                                                                                                               aaCTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGccc      <  1:2440253/148‑1 (MQ=255)
                                                                                                                                               aaCTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGccc      <  2:557002/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAGC  >  CP009273/1607646‑1607940

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: