Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,829,334 G→T 5.2% P134Q (CCG→CAG)  astD ← succinylglutamic semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,829,3340GT 5.2% 336.8 / 7.3 96P134Q (CCG→CAG) astDsuccinylglutamic semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (39/52);  new base T (3/2);  total (42/54)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

AAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCT  >  NC_000913/1829187‑1829482
                                                                                                                                                   |                                                                                                                                                    
aaaCGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCgg                                                                                                                                                     >  2:745461/1‑149 (MQ=255)
aaaCGCATTACCGCATCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCgg                                                                                                                                                     <  3:160500/149‑1 (MQ=255)
 aaCGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCggg                                                                                                                                                    >  1:727987/1‑149 (MQ=255)
    gCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGa                                                                                                                                                 <  2:319763/149‑1 (MQ=255)
     cATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGAt                                                                                                                                                <  3:977661/149‑1 (MQ=255)
          ccGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGc                                                                                                                                           >  4:924411/1‑149 (MQ=255)
           cGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCCGCACGGCGTGCTGGCGgtgtttgg                                                                                                                                          <  1:253740/149‑5 (MQ=255)
           cGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAACACGCCGTGCGGGCGATGTcaca                                                                                                                                          <  3:845361/149‑1 (MQ=255)
            gCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAg                                                                                                                                         >  3:59449/1‑149 (MQ=255)
            gCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCCATGTCGCAg                                                                                                                                         <  1:916092/149‑1 (MQ=255)
             cTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCGAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAgg                                                                                                                                        <  1:674819/149‑1 (MQ=255)
              ttCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGc                                                                                                                                       >  3:42296/1‑149 (MQ=255)
                  ccACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTcgc                                                                                                                                   >  4:418102/1‑149 (MQ=255)
                   cACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTcgcc                                                                                                                                  >  3:254659/1‑149 (MQ=255)
                    aCTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTcgccg                                                                                                                                 >  1:978315/1‑149 (MQ=255)
                    aCTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTcgcc                                                                                                                                  >  2:106288/1‑148 (MQ=255)
                    aCTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTcgcc                                                                                                                                  <  1:106288/148‑1 (MQ=255)
                     cTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTcgccgc                                                                                                                                <  2:824446/149‑1 (MQ=255)
                     cTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTcgccgc                                                                                                                                <  1:405035/149‑1 (MQ=255)
                     cTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTcgccgc                                                                                                                                <  2:1017683/149‑1 (MQ=255)
                        cACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTcgccgcgcc                                                                                                                             >  1:88480/1‑149 (MQ=255)
                        cACGGTGGCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTcgccgcgcc                                                                                                                             <  3:418102/149‑1 (MQ=255)
                            gtgtCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTcc                                                                                                                         <  3:259432/149‑1 (MQ=255)
                                    tCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTc                                                                                                                 <  2:100628/149‑1 (MQ=255)
                                      gCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCAc                                                                                                               >  3:424521/1‑149 (MQ=255)
                                       ctaggTTTAAAGATAATGGTGTGACCTGCCAGCAATGCCGGAACGATAGGTCCGTTCGGCAAATGGCCAGGGAAATTATACGGCCCAAACCCCGCCAGCACGCCGTGCGAGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACt                                                                                                              <  4:102990/146‑1 (MQ=255)
                                       cTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACt                                                                                                              <  3:890548/149‑1 (MQ=255)
                                         gggTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGGGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTAc                                                                                                            >  4:901842/1‑149 (MQ=255)
                                            tttAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACgct                                                                                                         >  2:167336/1‑149 (MQ=255)
                                            tttAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACgct                                                                                                         >  2:224695/1‑149 (MQ=255)
                                             ttAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCACGCCGTCCGGCATTTCACTACgctg                                                                                                        >  3:688233/1‑149 (MQ=255)
                                                aaGATAATGGTGTTACCTGCCAGCAATGCCGGCACGATATGACCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGCCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTc                                                                                                     <  4:853840/149‑1 (MQ=255)
                                                 aGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCg                                                                                                    >  3:619618/1‑149 (MQ=255)
                                                  gATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGc                                                                                                   >  4:239828/1‑149 (MQ=255)
                                                  gATAATGGTGTTACATGCCAGCATTGCCGGAACGAAATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCCCCGTCCGGTATTTCACTACGCTTCTCGc                                                                                                   >  4:780484/1‑149 (MQ=255)
                                                      aTGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGt                                                                                               >  1:220164/1‑149 (MQ=255)
                                                       tGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTa                                                                                              <  4:640448/149‑1 (MQ=255)
                                                         gtgtTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACg                                                                                            >  4:440535/1‑149 (MQ=255)
                                                            ttACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAAc                                                                                         <  3:202209/149‑1 (MQ=255)
                                                            ttACCTGCCAACAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAAc                                                                                         <  1:260246/149‑1 (MQ=255)
                                                               ccTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTTCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTg                                                                                      <  1:1014637/149‑1 (MQ=255)
                                                                cTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACATCGCCCGCACGGCGTGCTGGCGgtgttt                                                                                                                                            >  4:832992/1‑92 (MQ=37)
                                                                cTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACATCGCCCGCACGGCGTGCTGGCGgtgttt                                                                                                                                            <  3:832992/94‑3 (MQ=37)
                                                                cTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACgctg                                                                                                        <  4:388725/130‑1 (MQ=255)
                                                                cTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACgctg                                                                                                        >  3:388725/1‑130 (MQ=255)
                                                                cTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGa                                                                                     >  1:661832/1‑149 (MQ=255)
                                                                cTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGa                                                                                     <  2:604007/149‑1 (MQ=255)
                                                                   ccAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATAc                                                                                  >  4:763248/1‑149 (MQ=255)
                                                                   ccAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATAc                                                                                  >  1:396052/1‑149 (MQ=255)
                                                                      gCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGcc                                                                               <  4:232134/149‑1 (MQ=255)
                                                                          tGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTaa                                                                           <  4:424521/149‑1 (MQ=255)
                                                                           gCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAAt                                                                          <  3:763248/149‑1 (MQ=255)
                                                                           gCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAAt                                                                          <  4:447048/149‑1 (MQ=255)
                                                                            ccGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGGCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAAtt                                                                         <  3:851307/149‑1 (MQ=255)
                                                                              ggAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGa                                                                       <  1:171092/149‑1 (MQ=255)
                                                                              ggAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGa                                                                       >  4:241120/1‑149 (MQ=255)
                                                                              ggAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGa                                                                       >  3:1027013/1‑149 (MQ=255)
                                                                               gAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACgc                                                                                                          >  1:175042/1‑113 (MQ=255)
                                                                               gAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACgc                                                                                                          <  2:175042/113‑1 (MQ=255)
                                                                               gAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGat                                                                      <  2:911277/149‑1 (MQ=255)
                                                                               gAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGat                                                                      <  3:625918/149‑1 (MQ=255)
                                                                                aaCGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGata                                                                     >  2:745882/1‑149 (MQ=255)
                                                                                aaCGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCACCGGTACGAACGTGATACGCCTTAATTGata                                                                     <  2:765967/149‑1 (MQ=255)
                                                                                  cGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATc                                                                   <  3:634611/149‑1 (MQ=255)
                                                                                  cGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATc                                                                   <  1:943759/149‑1 (MQ=255)
                                                                                    atatGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATcgc                                                                 <  3:241120/149‑1 (MQ=255)
                                                                                      atGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAg                                                                                                                                         <  2:251316/75‑1 (MQ=255)
                                                                                      atGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAg                                                                                                                                         >  1:251316/1‑75 (MQ=255)
                                                                                      acgcccGTTCGGCAAATGCCCAGGGAAATTATACGGCCCAAACCCCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCCCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATGTCGCGa                                                               <  2:761361/145‑1 (MQ=255)
                                                                                        gTCCGTTCGGCAAATTACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACgc                                                                                                          <  1:82177/104‑1 (MQ=255)
                                                                                        gTCCGTTCGGCAAATTACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACgc                                                                                                          >  2:82177/1‑104 (MQ=255)
                                                                                        gTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTTATTGATATCGCGAtt                                                             <  4:1004395/149‑1 (MQ=255)
                                                                                           cGTTCGGCAAAAGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTa                                                          <  3:909980/149‑1 (MQ=255)
                                                                                             ttCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTAtt                                                        <  1:965088/149‑1 (MQ=255)
                                                                                                ggCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGat                                                     <  1:769648/149‑1 (MQ=255)
                                                                                                ggCAAATGACCAGGGAAATTAGACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACTTGATACGCCTTAATTGATATCGCGATTTTATTGat                                                     <  4:441872/149‑1 (MQ=255)
                                                                                                 gCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGatc                                                    <  1:900526/149‑1 (MQ=255)
                                                                                                         ccAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGt                                            <  2:978315/149‑1 (MQ=255)
                                                                                                          cAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTc                                           >  4:1039806/1‑149 (MQ=255)
                                                                                                             ggAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACt                                        >  3:1030955/1‑149 (MQ=255)
                                                                                                                aaTTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTc                                      >  2:361028/1‑148 (MQ=255)
                                                                                                                aaTTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTc                                      <  1:361028/148‑1 (MQ=255)
                                                                                                                  ttATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGt                                   <  1:1559/149‑1 (MQ=255)
                                                                                                                     taCGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGc                                <  4:619618/149‑1 (MQ=255)
                                                                                                                     taCGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGAGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGc                                <  2:915835/149‑1 (MQ=255)
                                                                                                                          cccAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCtt                           >  2:509125/1‑149 (MQ=255)
                                                                                                                           ccAAACACCGCCAGCACGCCGTGCGGGCTATTTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACATGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTc                          >  4:644559/1‑149 (MQ=255)
                                                                                                                           ccAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTc                          <  4:1027013/149‑1 (MQ=255)
                                                                                                                            cAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTcc                         >  1:435512/1‑149 (MQ=255)
                                                                                                                              aaCACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCa                       >  1:624045/1‑149 (MQ=255)
                                                                                                                              aaCACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCa                       >  1:942482/1‑149 (MQ=255)
                                                                                                                               acatCGCCCGCACGGCGTGCTGGCGGTGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAg                      >  4:183334/1‑149 (MQ=255)
                                                                                                                                     ccAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGc                <  2:208242/149‑1 (MQ=255)
                                                                                                                                      ccgCACGGCGTGCTGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCt               >  4:753094/3‑149 (MQ=255)
                                                                                                                                            gCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGat                                                                      <  3:399811/88‑1 (MQ=255)
                                                                                                                                            gCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGat                                                                      >  4:399811/1‑88 (MQ=255)
                                                                                                                                             ccGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGt        >  2:807977/1‑149 (MQ=255)
                                                                                                                                                  cGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTtctc   <  4:59449/149‑1 (MQ=255)
                                                                                                                                                  cGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTtctc   >  2:992944/1‑149 (MQ=255)
                                                                                                                                                   gggCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTctct  >  3:834895/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCT  >  NC_000913/1829187‑1829482

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: