Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,607,793 C→A 8.6% 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 8.6% 624.4 / 23.4 174E390* (GAG→TAG) sadsuccinate semialdehyde dehydrogenase, NAD(P)+‑dependent
Reads supporting (aligned to +/- strand):  ref base C (78/81);  new base A (6/9);  total (84/90)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 179 total aligned reads displayed.

AGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAGC  >  CP009273/1607646‑1607940
                                                                                                                                                   |                                                                                                                                                   
aGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTc                                                                                                                                                     >  2:641632/1‑148 (MQ=255)
 gCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTCAAAATGGGCGCTGAAAGGCCGAAGTct                                                                                                                                                    >  1:2861825/1‑147 (MQ=255)
     ctctTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTAt                                                                                                                                                >  1:1556513/1‑148 (MQ=255)
      tctTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATc                                                                                                                                               <  2:2926212/148‑1 (MQ=255)
        tttttCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCAt                                                                                                                                             >  1:1772267/1‑148 (MQ=255)
          tttCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTa                                                                                                                                           <  1:204221/148‑1 (MQ=255)
           ttCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAg                                                                                                                                          <  1:820824/148‑1 (MQ=255)
           ttCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAg                                                                                                                                          <  2:1201486/148‑1 (MQ=255)
             cACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGcc                                                                                                                                        <  2:2846313/148‑1 (MQ=255)
             cACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGcc                                                                                                                                        >  1:2739629/1‑148 (MQ=255)
               cGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAg                                                                                                                                      <  1:1304998/148‑1 (MQ=255)
                   accaAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTcc                                                                                                                                  <  2:2569276/148‑1 (MQ=255)
                    ccaAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCa                                                                                                                                 >  1:2175160/1‑148 (MQ=255)
                    ccaAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCa                                                                                                                                 >  1:2173981/1‑148 (MQ=255)
                      aaaGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGt                                                                                                                               >  1:1879633/1‑148 (MQ=255)
                      aaaGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGt                                                                                                                               >  2:730526/1‑148 (MQ=255)
                        aGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGc                                                                                                                             <  1:479840/148‑1 (MQ=255)
                         ggCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCa                                                                                                                            >  1:796253/1‑148 (MQ=255)
                            cACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGt                                                                                                                         <  1:1371384/148‑1 (MQ=255)
                            cACTCGCGCGGCGCTGGCACAATAACCATTGATAAACACCCCCCCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGt                                                                                                                         <  2:3154283/148‑1 (MQ=255)
                             aCTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGtt                                                                                                                        <  2:2801458/148‑1 (MQ=255)
                             aCTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGtt                                                                                                                        <  2:1151226/148‑1 (MQ=255)
                             aCTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTAACTATCATTAGCCAGTTCCAGTGCATGtt                                                                                                                        >  2:1836226/1‑148 (MQ=255)
                              cTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCTTGTTc                                                                                                                       <  1:3376/148‑1 (MQ=255)
                              cTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTc                                                                                                                       >  2:1922711/1‑148 (MQ=255)
                                cgcgcgTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTg                                                                                                                     >  1:181615/1‑148 (MQ=255)
                                    cgTCGCTGGCACAATAACCCTTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATTTTCTGCATc                                                                                                                 >  1:210281/1‑148 (MQ=255)
                                    cgTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATc                                                                                                                 <  1:340987/148‑1 (MQ=255)
                                             cacaATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAAt                                                                                                        >  2:2303584/1‑148 (MQ=255)
                                               caATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCTTGTTCTGCATCTTTCGCAATgg                                                                                                      <  2:1491464/148‑1 (MQ=255)
                                               caATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATgg                                                                                                      <  1:1662418/148‑1 (MQ=255)
                                                 aTAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTg                                                                                                    <  1:2075753/148‑1 (MQ=255)
                                                  tAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTAACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGa                                                                                                   >  2:2375941/1‑148 (MQ=255)
                                                   aaCCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGAt                                                                                                  >  1:2919649/1‑148 (MQ=255)
                                                   aaCCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGAt                                                                                                  >  1:393297/1‑148 (MQ=255)
                                                    aCCCTTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGAtt                                                                                                 <  1:2029611/148‑1 (MQ=255)
                                                     ccATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTAACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTg                                                                                                >  2:575490/1‑148 (MQ=255)
                                                            tAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAAc                                                                                         >  2:1080143/1‑148 (MQ=255)
                                                             aaaCACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACg                                                                                        <  2:387906/148‑1 (MQ=255)
                                                             aaaCACCCCACCGCAGTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACg                                                                                        <  1:1339490/148‑1 (MQ=255)
                                                               acacCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACggg                                                                                      <  1:3215328/148‑1 (MQ=255)
                                                               acacCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACggg                                                                                      <  1:3214007/148‑1 (MQ=255)
                                                               acacCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACggg                                                                                      >  2:744932/1‑148 (MQ=255)
                                                                cacCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACgggg                                                                                     >  1:6537/1‑148 (MQ=255)
                                                                     cACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCaaa                                                                                >  2:89862/1‑148 (MQ=255)
                                                                       ccGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCTTCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACa                                                                              <  2:20976/148‑1 (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)
                                                                             tCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATttctt                                                                        <  1:2886163/148‑1 (MQ=255)
                                                                              ccAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTAACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATttcttc                                                                       <  1:1472634/148‑1 (MQ=255)
                                                                               cAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTcc                                                                      >  2:2408675/1‑148 (MQ=255)
                                                                               cAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTcc                                                                      <  2:683696/148‑1 (MQ=255)
                                                                                 gACGTGCAGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACAATCATTAGCCAGATCCAGTGCATGTTCTGAATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTccct                                                                    >  2:929148/1‑147 (MQ=255)
                                                                                   cGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCTTGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGaa                                                                  <  1:3034622/148‑1 (MQ=255)
                                                                                   cGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGaa                                                                  <  1:1410944/148‑1 (MQ=255)
                                                                                    gTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGaaa                                                                 <  2:2433123/148‑1 (MQ=255)
                                                                                     tGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAAc                                                                >  2:477609/1‑148 (MQ=255)
                                                                                          ccATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGt                                                           >  2:1539031/1‑148 (MQ=255)
                                                                                          ccATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTAACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGt                                                           <  2:421006/148‑1 (MQ=255)
                                                                                           cATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTc                                                          >  2:3162507/1‑148 (MQ=255)
                                                                                            aTCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCTTGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCa                                                         <  1:1532384/148‑1 (MQ=255)
                                                                                            aTCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCa                                                         >  1:395734/1‑148 (MQ=255)
                                                                                             tctgtctgGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCAt                                                        >  1:456037/1‑148 (MQ=255)
                                                                                               tgtctgGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCAttt                                                      <  2:1590184/148‑1 (MQ=255)
                                                                                                  ctgGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTg                                                   >  2:1379382/1‑148 (MQ=255)
                                                                                                   tgGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTgg                                                  >  2:1236016/1‑148 (MQ=255)
                                                                                                   tgGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTgg                                                  >  2:2886164/1‑148 (MQ=255)
                                                                                                      ccTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGt                                               >  2:2751966/1‑148 (MQ=255)
                                                                                                        tgtgTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTATGGGGTaa                                             <  1:944565/148‑1 (MQ=255)
                                                                                                         gtgtTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAAc                                            >  2:648036/1‑148 (MQ=255)
                                                                                                          tgtTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACa                                           <  1:833290/148‑1 (MQ=255)
                                                                                                           gtTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCTTGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACAt                                          <  1:644816/148‑1 (MQ=255)
                                                                                                            tttCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACAtt                                         >  1:100690/1‑148 (MQ=255)
                                                                                                            tttCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACAtt                                         >  1:100160/1‑148 (MQ=255)
                                                                                                                gTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGcc                                     <  1:2975848/148‑1 (MQ=255)
                                                                                                                 tCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCa                                    >  2:1932114/1‑148 (MQ=255)
                                                                                                                 tCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCa                                    >  2:2625061/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)
                                                                                                                       gTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCg                              <  1:1625649/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)
                                                                                                                                 tCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCTTGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGAta                    >  2:12851/1‑148 (MQ=255)
                                                                                                                                 tCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCTTGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGAta                    >  2:13262/1‑148 (MQ=255)
                                                                                                                                   gCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGAtagt                  <  1:2761798/148‑1 (MQ=255)
                                                                                                                                       aaaGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGtt              >  2:588242/1‑148 (MQ=255)
                                                                                                                                       aaaGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGtt              >  2:588253/1‑148 (MQ=255)
                                                                                                                                        aaGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTa             >  1:2625109/1‑148 (MQ=255)
                                                                                                                                          ggCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTAcc           >  1:315385/1‑148 (MQ=255)
                                                                                                                                             cGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGc        >  1:575186/1‑148 (MQ=255)
                                                                                                                                             cGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGc        <  2:3099359/148‑1 (MQ=255)
                                                                                                                                               aaCTCACTATCATTAGCCAGTTCCAGTGCTTGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGccc      <  1:2510777/148‑1 (MQ=255)
                                                                                                                                               aaCTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGccc      >  1:2631280/1‑148 (MQ=255)
                                                                                                                                                aCTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGcccc     <  2:2975058/148‑1 (MQ=255)
                                                                                                                                                  tCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATTTCTTCCCGAAACGCGGTCATTTCTGGGGTAACATTCGCCAGAACCGTTGGCGGATAGTAGTTACCTGCCCCAg   >  1:2921034/1‑148 (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: ‑