Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,823,843 A→C 9.1% F50C (TTT→TGT)  astD ← succinylglutamic semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,823,8430AC 9.1% 1083.4 / 21.1 379F50C (TTT→TGT) astDsuccinylglutamic semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (225/108);  new base C (3/43);  total (228/151)
Fisher's exact test for biased strand distribution p-value = 4.97e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 391 total aligned reads displayed.

CCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATG  >  CP009273/1823694‑1823992
                                                                                                                                                     |                                                                                                                                                     
  gTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAAc                                                                                                                                                     >  2:284008/1‑150 (MQ=255)
  gTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTGCAGCAAATGAAAGCCGCGCCCAGCGCGGAAAc                                                                                                                                                     >  2:4771478/1‑150 (MQ=255)
   tCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAAcg                                                                                                                                                    >  1:2784834/1‑150 (MQ=255)
     aCTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGAATGACGTTCAGCAATTGAAAGCCGCGCCCAGCGCGGAAAcgcc                                                                                                                                                  <  1:2059957/150‑1 (MQ=255)
        tCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCa                                                                                                                                               >  1:2394465/1‑150 (MQ=255)
          ggTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACg                                                                                                                                             >  1:4870646/1‑150 (MQ=255)
          ggTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACg                                                                                                                                             >  2:2698893/1‑150 (MQ=255)
                cGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGTCTTTATTGATTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCAGCGCCCAGCGCGGAAACGCCGCACGGGCTGc                                                                                                                                       >  2:2791055/1‑150 (MQ=255)
                            cgGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACACGCCTg                                                                                                                           >  2:3274303/1‑150 (MQ=255)
                             ggCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGACACGCCGCACGGGCTGCCCGACAAGCCTGc                                                                                                                          >  2:254313/1‑150 (MQ=255)
                             ggCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGc                                                                                                                          >  1:2839369/1‑150 (MQ=255)
                             ggCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGc                                                                                                                          >  1:1365170/1‑150 (MQ=255)
                                ttACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCg                                                                                                                       >  1:307922/1‑150 (MQ=255)
                                  aCCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGAc                                                                                                                     >  1:5247398/1‑150 (MQ=255)
                                   cccGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGAcc                                                                                                                    >  1:1590180/1‑150 (MQ=255)
                                     cGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTg                                                                                                                  >  1:2991630/1‑150 (MQ=255)
                                       tttCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAg                                                                                                                >  1:5128539/1‑150 (MQ=255)
                                         tctcTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcg                                                                                                              >  1:3674468/1‑150 (MQ=255)
                                         tctcTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcg                                                                                                              >  1:145324/1‑150 (MQ=255)
                                         tctcTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcg                                                                                                              >  2:1740618/1‑150 (MQ=255)
                                           tctGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggc                                                                                                            >  2:2108269/1‑150 (MQ=255)
                                           tctGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggc                                                                                                            >  1:639398/1‑150 (MQ=255)
                                              ggCAATAATCGCGGTTAATTCGTCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACTGCATGACGTTCAGCAAATGAAAGCAGCGCCCAGCTCGTAAACGACGCACTGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcatc                                                                                                         >  2:1977087/1‑150 (MQ=255)
                                               gCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcatcg                                                                                                        >  2:2181929/1‑150 (MQ=255)
                                               caaataatCGCGTTTTATCCGGCTTTTTGGCTTCCCAGCGGTGCGCCACTGGGTTCGACACCGGCATGACTTTCACAAAAGGAACCCCGCCCCCCGCGCGCACACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcatcg                                                                                                        <  1:3739074/148‑1 (MQ=255)
                                                cAATATTGGCGTTTAATTCGGCTTTATTGCTTTCCCGCAGTGCGGCAACGCTTTCGACACCGCCATGACGTTCAGAAAATGAAAGCCGCGCCCTGCGCGTAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcatcgg                                                                                                       <  1:4733543/150‑1 (MQ=255)
                                                cAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCATCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcatcgg                                                                                                       >  1:2226567/1‑150 (MQ=255)
                                                    aaacgcgGTTAATACGGCTTTATTGCTTACCAGCAGTGAGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCCCCAAGCGCTGAAACGCCGCACGGGCTGCCCGACAAGCCAGCTCGACCTGAGAgtcatcggaatc                                                                                                   >  2:1462516/4‑150 (MQ=255)
                                                        gggtttAATTCGGCTTTTTTGCTTTCCAGCAGTGCGGCAACGCGTTCGACACCGGCATGACGTTCAGCAAATGAAAGCCGCGCCCCGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttg                                                                                               <  2:4123703/146‑1 (MQ=255)
                                                            ttAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttgcctt                                                                                           >  2:3378290/1‑150 (MQ=255)
                                                             tAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttgccttg                                                                                          >  1:3232936/1‑150 (MQ=255)
                                                                    gCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAAGGGCATGACGTTCAGCAAAGGAAAGCCGCGCCCCGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCCTTGCCTTGCCATAac                                                                                   <  2:2476283/150‑1 (MQ=255)
                                                                     cTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCCGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAaca                                                                                  <  2:2614155/150‑1 (MQ=255)
                                                                      tttATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGACACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAacac                                                                                 >  2:1871733/1‑150 (MQ=255)
                                                                       ttATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACAcc                                                                                >  1:804486/1‑150 (MQ=255)
                                                                        tATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCt                                                                               >  2:1180100/1‑150 (MQ=255)
                                                                         tttGCTTTCCAGCAGTGCGGCAAAGGGTTCGACAACGGAATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTc                                                                              <  1:2565205/149‑1 (MQ=255)
                                                                          ttGCTTTCCCGCAGTGCGGAAATGCGTTCGAAAACGGCATGACGTTCAGAAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCg                                                                             <  2:632304/150‑1 (MQ=255)
                                                                          ttGCTTTCCCGCAGTGCGGAAACGCGTTCGACACCGGCATGACGTTCAGCAAATGAAAGCCGCGCCCCGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCg                                                                             <  2:3114278/150‑1 (MQ=255)
                                                                          ttGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCg                                                                             >  2:4606105/1‑150 (MQ=255)
                                                                          ttGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCg                                                                             >  2:1782214/1‑150 (MQ=255)
                                                                           tGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGc                                                                            >  1:3319353/1‑150 (MQ=255)
                                                                             cTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGccc                                                                          >  1:425090/1‑150 (MQ=255)
                                                                               ttCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGa                                                                        >  2:3454412/1‑150 (MQ=255)
                                                                                 cccgctgTGCGGCAAAGCTTTCGACAACGGCTTGACGTTCAGCAAATGAACGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATa                                                                      <  2:606661/144‑1 (MQ=255)
                                                                                 ccAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATa                                                                      >  1:2689938/1‑150 (MQ=255)
                                                                                  cagcagTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATAc                                                                     >  1:1884172/1‑150 (MQ=255)
                                                                                  cagcagTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATAc                                                                     >  1:1647491/1‑150 (MQ=255)
                                                                                   agccgggCGGCAAAGCGTTCGACAACGGAATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATAcc                                                                    <  2:2252307/144‑1 (MQ=255)
                                                                                    gcagTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAGTGCCTTGCCATAACACCTCGCCCGATACCg                                                                   <  2:2784815/150‑1 (MQ=255)
                                                                                      agTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGa                                                                 >  2:3741036/1‑150 (MQ=255)
                                                                                       gggCGGCAATGGGTTCGACACGGGAATGAGTTTCAGCAATTGAAAGCCGCGCCCAGGGGGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGAt                                                                <  1:4771478/148‑1 (MQ=255)
                                                                                       gTGCGGCAACGGGTTCGACACGGGATTTAGGTTCAGAAAATCAAAGCCGCGCCCCGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGAt                                                                <  1:729206/150‑1 (MQ=255)
                                                                                         ggggaaaaGGGTTCGAAATCGGATTTACGTTCAGAAAATGAATGCCGCGCCCCGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTa                                                              <  1:1167128/145‑1 (MQ=255)
                                                                                         gCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATTAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTa                                                              <  1:1262865/150‑1 (MQ=255)
                                                                                           ggCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACg                                                            >  1:1539257/1‑150 (MQ=255)
                                                                                           ggCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACg                                                            <  2:2731928/150‑1 (MQ=255)
                                                                                             cAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCt                                                          >  1:3974583/1‑150 (MQ=255)
                                                                                              aaaGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCtt                                                         >  2:3233232/1‑150 (MQ=255)
                                                                                               aaGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTc                                                        >  2:4598614/1‑150 (MQ=255)
                                                                                               aaGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTc                                                        >  1:2596490/1‑150 (MQ=255)
                                                                                                 gCGTTCGACACCGGCATGAGGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTcac                                                      <  2:2318557/150‑1 (MQ=255)
                                                                                                  ggtgcGACACGGGCTTGACGTGCATAAAATGAAAGCCGCGCCCCGCGGGGACACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTcaca                                                     <  1:2911031/146‑1 (MQ=255)
                                                                                                  cGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTcaca                                                     >  1:573720/1‑150 (MQ=255)
                                                                                                  cGTTCGACAACGGCATGACGTTCAGCAAAGGAACGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTcaca                                                     <  2:4403238/150‑1 (MQ=255)
                                                                                                   gTTCTACAACGGATTGACGTTCACCAAATGAACGCCCCCCCCCCCGCGTAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCacac                                                    <  1:2493114/150‑1 (MQ=255)
                                                                                                    ttCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACAcg                                                   <  2:4171395/150‑1 (MQ=255)
                                                                                                       gACACCGGCATGACGTTCAGCAAATGAAAGCCGCGCCCCGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGa                                                <  2:123714/150‑1 (MQ=255)
                                                                                                       gACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGa                                                >  2:5085224/1‑150 (MQ=255)
                                                                                                             ggCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGcc                                          >  1:4574126/1‑150 (MQ=255)
                                                                                                                aTGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTg                                       >  2:1085393/1‑150 (MQ=255)
                                                                                                                aTGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTg                                       >  2:4980061/1‑150 (MQ=255)
                                                                                                                   aCGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGcc                                    >  1:113861/1‑150 (MQ=255)
                                                                                                                        cAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTa                               >  2:2021328/1‑150 (MQ=255)
                                                                                                                        cAGAAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTa                               <  1:4648077/150‑1 (MQ=255)
                                                                                                                          gCAAATGACCGCCGCGCCCCGCGCGGACACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATc                             <  2:4506285/150‑1 (MQ=255)
                                                                                                                          gCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACTAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATc                             <  2:1816152/150‑1 (MQ=255)
                                                                                                                           cAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATcc                            >  1:2288752/1‑150 (MQ=255)
                                                                                                                            aaaTGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCa                           >  1:3040131/1‑150 (MQ=255)
                                                                                                                             aaTGAAAGCCGCGCCCAGCGATGAAACTCCGCACGGGCTGCCCGACAAGCATGCTCGACCTGAGCGGCATCGGAATCATAGCCATGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGAGCCCAGGCCCGTTATCCAg                          >  2:2493866/1‑150 (MQ=255)
                                                                                                                              aTGAAAGCCGCGCCCAGCGCGGACACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGt                         <  2:4453238/150‑1 (MQ=255)
                                                                                                                                  aatccGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCAcc                     >  1:3440517/4‑150 (MQ=255)
                                                                                                                                  aaGCCGCGCCCAGCGCGGACACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCAcc                     <  2:3587666/150‑1 (MQ=255)
                                                                                                                                      cgccCCCCGCGCGCACACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTa                 <  1:4425709/150‑1 (MQ=255)
                                                                                                                                       gcgcCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTaa                <  1:2108269/150‑1 (MQ=255)
                                                                                                                                       gcgcCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTaa                <  2:842568/150‑1 (MQ=255)
                                                                                                                                       gcgcCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTaa                <  2:2404760/150‑1 (MQ=255)
                                                                                                                                        cgcCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAAt               <  1:5084668/150‑1 (MQ=255)
                                                                                                                                          cccAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATcc             >  1:357987/1‑150 (MQ=255)
                                                                                                                                           ccAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCa            <  2:397585/150‑1 (MQ=255)
                                                                                                                                            cAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGCTGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCAt           <  2:2967311/150‑1 (MQ=255)
                                                                                                                                             aGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATa          >  2:503925/1‑150 (MQ=255)
                                                                                                                                              gcgcGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATaa         <  2:1545009/150‑1 (MQ=255)
                                                                                                                                                gcgGACACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTGTCCAGTCACCGTTAATCCATAAAg       <  1:613068/150‑1 (MQ=255)
                                                                                                                                                 cgGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGt      >  2:4215308/1‑150 (MQ=255)
                                                                                                                                                 cgGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGt      >  2:759693/1‑150 (MQ=255)
                                                                                                                                                  ggACACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTc     <  1:704883/150‑1 (MQ=255)
                                                                                                                                                    aaaCGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCAt   <  2:4291574/150‑1 (MQ=255)
                                                                                                                                                    aaaCGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCAt   >  2:4877525/1‑150 (MQ=255)
                                                                                                                                                     aaCGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATg  >  2:347370/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATG  >  CP009273/1823694‑1823992

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑