Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,829,328 C→G 5.0% G136A (GGC→GCC)  astD ← succinylglutamic semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,829,3280CG 5.0% 346.8 / 6.3 99G136A (GGC→GCC) astDsuccinylglutamic semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (52/42);  new base G (3/2);  total (55/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TGCCATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCG  >  NC_000913/1829181‑1829475
                                                                                                                                                   |                                                                                                                                                   
tgcCATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGcc                                                                                                                                                    >  2:802841/1‑149 (MQ=255)
  ccATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGt                                                                                                                                                  <  4:676588/149‑1 (MQ=255)
   cATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTg                                                                                                                                                 >  2:82447/1‑149 (MQ=255)
    aTAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGc                                                                                                                                                >  3:5685/1‑149 (MQ=255)
      aaaCGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCgg                                                                                                                                              >  3:105179/1‑149 (MQ=255)
      aaaCGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCgg                                                                                                                                              >  2:87175/1‑149 (MQ=255)
       aaCGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTTCggg                                                                                                                                             >  4:220340/1‑149 (MQ=255)
       aaCGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCggg                                                                                                                                             >  1:853975/1‑149 (MQ=255)
       aaCGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCggg                                                                                                                                             <  3:847106/149‑1 (MQ=255)
       aaCGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCggg                                                                                                                                             >  4:421663/1‑149 (MQ=255)
        aCGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCCGCACGGCGTGCTGGc                                                                                                                                            <  4:840677/149‑1 (MQ=255)
           cATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGAt                                                                                                                                         >  3:63814/1‑149 (MQ=255)
           cATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGAt                                                                                                                                         <  1:583838/149‑1 (MQ=255)
             ttACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGt                                                                                                                                       >  1:292755/1‑149 (MQ=255)
              tACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTc                                                                                                                                      <  4:575572/149‑1 (MQ=255)
               aCCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCg                                                                                                                                     >  2:642866/1‑149 (MQ=255)
                 cGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCACACACCGCCAGCACGCCGTGCGGGCGATGtttgg                                                                                                                                   <  3:584566/149‑5 (MQ=255)
                 cGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCa                                                                                                                                   <  4:832376/149‑1 (MQ=255)
                 cGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCa                                                                                                                                   <  2:327151/149‑1 (MQ=255)
                  gCTTCGCCACTCCACGTTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACATCGCCCGCACGGCGTGCTGGCGgtgtttggg                                                                                                                                  <  3:863054/149‑6 (MQ=255)
                  gCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAg                                                                                                                                  <  2:765342/149‑1 (MQ=255)
                      cGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTc                                                                                                                              >  2:538245/1‑149 (MQ=255)
                            tCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTcgccgcg                                                                                                                        <  3:421663/149‑1 (MQ=255)
                            tCCACGGTGTCAGTTCGCTGGGTTTAAAGATAAAGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTcgccgcg                                                                                                                        >  4:857297/1‑149 (MQ=255)
                              cACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTcgccgcgcc                                                                                                                      >  3:705058/1‑149 (MQ=255)
                                 ggTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTc                                                                                                                   <  4:246576/149‑1 (MQ=255)
                                  gtgtCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTcc                                                                                                                  >  2:920118/1‑149 (MQ=255)
                                  gtgtCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTcc                                                                                                                  >  3:681477/1‑149 (MQ=255)
                                     tCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGc                                                                                                               <  1:753326/149‑1 (MQ=255)
                                         ttCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCAttt                                                                                                           >  2:147354/1‑149 (MQ=255)
                                         ttCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCAttt                                                                                                           >  4:331430/1‑149 (MQ=255)
                                            gCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCAc                                                                                                        <  1:571967/149‑1 (MQ=255)
                                              tGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTa                                                                                                      >  2:356162/1‑149 (MQ=255)
                                              tGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTa                                                                                                      >  1:597416/1‑149 (MQ=255)
                                                           aaTGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCgg                                                                                         >  4:387752/1‑149 (MQ=255)
                                                            aTGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGt                                                                                        >  3:273983/1‑149 (MQ=255)
                                                             tGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGagc                                                                                                                                        <  3:588793/100‑3 (MQ=255)
                                                             tGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGagc                                                                                                                                        >  4:588793/1‑98 (MQ=255)
                                                             tGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTa                                                                                       <  2:597416/149‑1 (MQ=255)
                                                             tGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTa                                                                                       >  1:612832/1‑149 (MQ=255)
                                                             tGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTa                                                                                       >  1:389913/1‑149 (MQ=255)
                                                               gtgtTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACg                                                                                     >  4:141064/1‑149 (MQ=255)
                                                                  ttACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCa                                                                                                                                   <  4:474123/100‑1 (MQ=255)
                                                                  ttACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCa                                                                                                                                   >  3:474123/1‑100 (MQ=255)
                                                                  ttACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAAc                                                                                  >  2:560371/1‑149 (MQ=255)
                                                                     ccTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAGATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTg                                                                               >  1:766642/1‑149 (MQ=255)
                                                                      cTGCCAGCAATGCCGGAACGCTATGTACGTTACGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGa                                                                              >  2:752330/1‑149 (MQ=255)
                                                                        gCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATa                                                                            <  4:105179/149‑1 (MQ=255)
                                                                        gCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATa                                                                            >  4:865861/1‑149 (MQ=255)
                                                                          cAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACg                                                                          >  1:231832/1‑149 (MQ=255)
                                                                            gCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGcc                                                                        <  3:857297/149‑1 (MQ=255)
                                                                              aaTGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCtt                                                                      >  4:412749/1‑149 (MQ=255)
                                                                                         gATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATcg                                                           >  2:503009/1‑149 (MQ=255)
                                                                                          atatGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATcgc                                                          <  3:517396/149‑1 (MQ=255)
                                                                                          atatGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATcgc                                                          >  2:687271/1‑149 (MQ=255)
                                                                                             tGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGAt                                                       <  4:273983/149‑1 (MQ=255)
                                                                                              gTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGAtt                                                      <  1:793905/149‑1 (MQ=255)
                                                                                               tCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAAt                                                                   >  4:608381/1‑135 (MQ=255)
                                                                                               tCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGAttt                                                     >  3:151789/1‑149 (MQ=255)
                                                                                               tCCGTTCGGCAAATGACCAGGGAAAGTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGGTACGCCTTAAt                                                                   <  3:608381/135‑1 (MQ=255)
                                                                                                 cGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTa                                                   <  4:516204/149‑1 (MQ=255)
                                                                                                  gTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTAt                                                  <  1:356162/149‑1 (MQ=255)
                                                                                                  gTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTAt                                                  <  1:82447/149‑1 (MQ=255)
                                                                                                   ttCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTAtt                                                 <  1:920118/149‑1 (MQ=255)
                                                                                                   ttCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTAtt                                                 <  1:538245/149‑1 (MQ=255)
                                                                                                   ttCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTAtt                                                 <  2:555246/149‑1 (MQ=255)
                                                                                                   ttCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTAtt                                                 <  2:74702/149‑1 (MQ=255)
                                                                                                   ttCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGACCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTAtt                                                 <  4:63814/149‑1 (MQ=255)
                                                                                                     cGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGa                                               <  4:384679/149‑1 (MQ=255)
                                                                                                      ggCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGat                                              >  3:869885/1‑149 (MQ=255)
                                                                                                       gCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGatc                                             <  2:878127/149‑1 (MQ=255)
                                                                                                       gCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGatc                                             >  1:878127/1‑149 (MQ=255)
                                                                                                         aaaTGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGa                                                        >  3:51981/1‑136 (MQ=255)
                                                                                                         aaaTGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGa                                                        <  4:51981/136‑1 (MQ=255)
                                                                                                           aTGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCg                                         <  3:865861/149‑1 (MQ=255)
                                                                                                              aCCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCg                                      >  1:190763/1‑149 (MQ=255)
                                                                                                              aCCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCg                                      >  4:555315/1‑149 (MQ=255)
                                                                                                                 aGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCa                                   <  4:849765/149‑1 (MQ=255)
                                                                                                                 aGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCa                                   >  2:819550/1‑149 (MQ=255)
                                                                                                                 aGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCa                                   <  4:127572/149‑1 (MQ=255)
                                                                                                                  gggAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCAc                                  >  3:476680/1‑149 (MQ=255)
                                                                                                                   ggAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACt                                 >  2:419773/1‑149 (MQ=255)
                                                                                                                    gAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTACTTGATATCGCGATTTTATTGCTCATCGCCGTCACtt                                >  1:883475/1‑149 (MQ=255)
                                                                                                                      aaTTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCg                              <  4:681477/149‑1 (MQ=255)
                                                                                                                       aTTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCgg                             <  1:642866/149‑1 (MQ=255)
                                                                                                                        ttATATGGCCCAAACACCGCCAGCACGCCGTGCCGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGATCGTGATACGCCTTAATTGATATCGCGATTTTAtt                                                 >  1:268028/1‑128 (MQ=255)
                                                                                                                        ttATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTAtt                                                 <  2:268028/128‑1 (MQ=255)
                                                                                                                        ttATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGt                            >  4:727869/1‑149 (MQ=255)
                                                                                                                         tataCGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCTCTTCGGtt                           <  1:355243/149‑1 (MQ=255)
                                                                                                                          ataCGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTg                          <  1:87175/149‑1 (MQ=255)
                                                                                                                           taCGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGc                         >  1:503155/1‑149 (MQ=255)
                                                                                                                             cGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCg                       <  3:730804/149‑1 (MQ=255)
                                                                                                                               gCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCt                     >  2:795844/1‑149 (MQ=255)
                                                                                                                                 ccAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTc                   >  4:204276/1‑149 (MQ=255)
                                                                                                                                 ccAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTc                   >  2:426028/1‑149 (MQ=255)
                                                                                                                                     acacCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAg               <  2:766642/149‑1 (MQ=255)
                                                                                                                                       atcGCCCGCACGGCGTGCTGGCGGTGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAgcg             >  2:936354/3‑149 (MQ=255)
                                                                                                                                         cGCCCGCACGGCGTGCTGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGcgcg           >  4:563638/1‑149 (MQ=255)
                                                                                                                                           cccgCACGGCGTGCTGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGc         >  4:882703/4‑149 (MQ=255)
                                                                                                                                                  gCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCg  >  4:560145/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGCCATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCG  >  NC_000913/1829181‑1829475

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