Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,607,793 C→A 6.3% 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 6.3% 773.0 / 17.7 204E390* (GAG→TAG) sadsuccinate semialdehyde dehydrogenase, NAD(P)+‑dependent
Reads supporting (aligned to +/- strand):  ref base C (95/96);  new base A (4/9);  total (99/105)
Fisher's exact test for biased strand distribution p-value = 2.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 206 total aligned reads displayed.

AGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAG  >  CP009273/1607646‑1607939
                                                                                                                                                   |                                                                                                                                                  
aGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTc                                                                                                                                                    >  2:194766/1‑148 (MQ=255)
 gCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCa                                                                                                                                                   <  1:1925870/148‑1 (MQ=255)
  ccACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCAc                                                                                                                                                  >  1:934776/1‑148 (MQ=255)
  ccACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCAc                                                                                                                                                  >  1:11496/1‑148 (MQ=255)
   cACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACt                                                                                                                                                 <  2:361808/148‑1 (MQ=255)
      tctTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATc                                                                                                                                              >  2:1372420/1‑148 (MQ=255)
        tttttCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCAt                                                                                                                                            <  2:314749/148‑1 (MQ=255)
         ttttCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCAtt                                                                                                                                           <  2:491075/148‑1 (MQ=255)
            tCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTAACTATCATTAGc                                                                                                                                        <  2:2325511/148‑1 (MQ=255)
             cACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGcc                                                                                                                                       <  1:2569474/148‑1 (MQ=255)
                  caccaAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTc                                                                                                                                  <  2:1225404/148‑1 (MQ=255)
                   accaAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTcc                                                                                                                                 >  2:1601082/1‑148 (MQ=255)
                    ccaAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCa                                                                                                                                >  2:1388591/1‑148 (MQ=255)
                    ccaAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCa                                                                                                                                <  1:2732442/148‑1 (MQ=255)
                        aGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGc                                                                                                                            >  1:3503783/1‑148 (MQ=255)
                        aGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTAACTATCATTAGCCAGTTCCAGTGc                                                                                                                            >  1:3566829/1‑148 (MQ=255)
                         ggCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCa                                                                                                                           >  1:2668756/1‑148 (MQ=255)
                          gCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCAt                                                                                                                          <  1:1320896/148‑1 (MQ=255)
                          gCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCAt                                                                                                                          >  2:1130047/1‑148 (MQ=255)
                                  cgcgTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCa                                                                                                                  >  1:1344575/1‑148 (MQ=255)
                                  cgcgTCGCTGGCACAATAAACATTTATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTTTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACACACTATCATTAGCCACTTCCAGTGCATGTTCTGCa                                                                                                                  <  2:2134143/148‑1 (MQ=255)
                                      tCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCtt                                                                                                              <  2:2157772/148‑1 (MQ=255)
                                      tCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCtt                                                                                                              <  1:2478415/148‑1 (MQ=255)
                                      tCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCtt                                                                                                              >  1:1809613/1‑148 (MQ=255)
                                       cGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCttt                                                                                                             >  2:2384389/1‑148 (MQ=255)
                                        gCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTc                                                                                                            >  2:3186388/1‑148 (MQ=255)
                                           ggCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCa                                                                                                         <  1:1053172/148‑1 (MQ=255)
                                              acaATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATg                                                                                                      <  1:2293772/148‑1 (MQ=255)
                                              acaATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATg                                                                                                      <  1:2252328/148‑1 (MQ=255)
                                                 aTAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTg                                                                                                   <  1:678978/148‑1 (MQ=255)
                                                  tAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGa                                                                                                  <  1:2458204/148‑1 (MQ=255)
                                                  tAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGa                                                                                                  >  2:2571441/1‑148 (MQ=255)
                                                  tAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTAACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGa                                                                                                  <  2:2013759/148‑1 (MQ=255)
                                                     ccATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTg                                                                                               >  2:2824468/1‑148 (MQ=255)
                                                     ccATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTg                                                                                               <  1:1472318/148‑1 (MQ=255)
                                                      cATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGc                                                                                              >  2:2118706/1‑148 (MQ=255)
                                                      cATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGc                                                                                              >  1:2182940/1‑148 (MQ=255)
                                                      cATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGc                                                                                              <  2:834629/148‑1 (MQ=255)
                                                       aTTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACAATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGcc                                                                                             >  1:590636/1‑148 (MQ=255)
                                                        ttGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCg                                                                                            >  2:3344518/1‑148 (MQ=255)
                                                           aTAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCAAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCaa                                                                                         <  1:2910727/148‑1 (MQ=255)
                                                             aaaCACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACg                                                                                       >  2:3489353/1‑148 (MQ=255)
                                                             aaaCACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACg                                                                                       >  1:2497676/1‑148 (MQ=255)
                                                               acacCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACggg                                                                                     >  1:3094637/1‑148 (MQ=255)
                                                               acacCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACggg                                                                                     >  2:1478368/1‑148 (MQ=255)
                                                               acacCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACggg                                                                                     >  2:828854/1‑148 (MQ=255)
                                                                cacCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAATGCCGAACTCACTATCATTAGCCAGGTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACgggg                                                                                    <  2:3011245/148‑1 (MQ=255)
                                                                  ccccACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGcc                                                                                  >  1:1299549/1‑148 (MQ=255)
                                                                    cccccGCCTTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCaa                                                                                <  1:2022312/145‑1 (MQ=255)
                                                                    ccACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCTTGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGCCCaa                                                                                >  2:2135683/1‑148 (MQ=255)
                                                                    ccACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCaa                                                                                >  1:501622/1‑148 (MQ=255)
                                                                      aCCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAAc                                                                              >  1:1319736/1‑148 (MQ=255)
                                                                       ccGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACa                                                                             >  1:982596/1‑148 (MQ=255)
                                                                       ccGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTAACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACa                                                                             <  1:491010/148‑1 (MQ=255)
                                                                           aTTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATttc                                                                         <  1:1466723/148‑1 (MQ=255)
                                                                           aTTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATttc                                                                         <  1:1465687/148‑1 (MQ=255)
                                                                            ttCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATttct                                                                        >  1:1909929/1‑148 (MQ=255)
                                                                            ttCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATttct                                                                        <  1:432284/148‑1 (MQ=255)
                                                                            ttCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATttct                                                                        <  1:3538620/148‑1 (MQ=255)
                                                                            ttCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATttct                                                                        >  1:109158/1‑148 (MQ=255)
                                                                               cAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTcc                                                                     >  1:3426237/1‑148 (MQ=255)
                                                                                aGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTccc                                                                    >  1:1480130/1‑148 (MQ=255)
                                                                                      gccgccATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAAcg                                                              <  1:3452300/148‑1 (MQ=255)
                                                                                          ccATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGt                                                          >  2:3057396/1‑148 (MQ=255)
                                                                                          ccATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGt                                                          <  1:564489/148‑1 (MQ=255)
                                                                                           cATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTc                                                         >  1:3491927/1‑148 (MQ=255)
                                                                                            aTCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCa                                                        >  1:1006183/1‑148 (MQ=255)
                                                                                             tctgtctgGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCAt                                                       <  1:1372420/148‑1 (MQ=255)
                                                                                              ctgtctgGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCAtt                                                      >  2:2597271/1‑148 (MQ=255)
                                                                                                gtctgGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTc                                                    >  1:822887/1‑148 (MQ=255)
                                                                                                  ctgGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTg                                                  <  2:1809613/148‑1 (MQ=255)
                                                                                                  ctgGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTg                                                  >  1:3029804/1‑148 (MQ=255)
                                                                                                    ggCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTggg                                                >  2:600136/1‑148 (MQ=255)
                                                                                                    ggCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTggg                                                <  1:2312863/148‑1 (MQ=255)
                                                                                                       cTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTa                                             <  1:1933573/148‑1 (MQ=255)
                                                                                                        tgtgTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTaa                                            >  1:392846/1‑148 (MQ=255)
                                                                                                            tttCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCGTCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACAtt                                        >  1:2102308/1‑148 (MQ=255)
                                                                                                            tttCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACAtt                                        >  2:2721899/1‑148 (MQ=255)
                                                                                                                gTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGcc                                    <  1:752642/148‑1 (MQ=255)
                                                                                                                 tCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCa                                   >  1:1534648/1‑148 (MQ=255)
                                                                                                                  cAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAg                                  >  2:2756293/1‑148 (MQ=255)
                                                                                                                     tGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAAc                               >  2:2677597/1‑148 (MQ=255)
                                                                                                                       gTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCg                             <  2:3470855/148‑1 (MQ=255)
                                                                                                                       gTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCg                             <  2:151870/148‑1 (MQ=255)
                                                                                                                        tAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGt                            >  1:2268916/1‑148 (MQ=255)
                                                                                                                         aaaaaTGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGtt                           <  1:1677454/148‑1 (MQ=255)
                                                                                                                           aaaTGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTgg                         <  2:872087/148‑1 (MQ=255)
                                                                                                                           aaaTGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTgg                         >  1:3066067/1‑148 (MQ=255)
                                                                                                                              tGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGACCCGTTGGCgg                      >  2:1037633/1‑148 (MQ=255)
                                                                                                                                gTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGAt                    <  1:2062883/148‑1 (MQ=255)
                                                                                                                                 tCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGAta                   <  2:2797280/148‑1 (MQ=255)
                                                                                                                                  cGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGAtag                  <  2:934776/148‑1 (MQ=255)
                                                                                                                                    cTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGAtagta                >  2:1639294/1‑148 (MQ=255)
                                                                                                                                     tGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGAtagtag               >  1:1552405/1‑148 (MQ=255)
                                                                                                                                       aaaGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGtt             <  2:3426237/148‑1 (MQ=255)
                                                                                                                                       aaaGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGtt             >  1:472611/1‑148 (MQ=255)
                                                                                                                                         aGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTACTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTAc           <  2:2881586/148‑1 (MQ=255)
                                                                                                                                            ccGAACTCACTATCATTAGCCAGTTCCAGTGCTTGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTg        >  1:2937472/1‑148 (MQ=255)
                                                                                                                                              gAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGcc      <  1:3159438/148‑1 (MQ=255)
                                                                                                                                               aaCTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGccc     >  2:2737267/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAG  >  CP009273/1607646‑1607939

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: