Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,422,841 C→A 100% D293Y (GAT→TAT)  purF ← amidophosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,422,8410CA98.6% 734.6 / 1.3 213D293Y (GAT→TAT) purFamidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base A (98/112);  total (100/113)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 213 total aligned reads displayed.

GCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAGCAC  >  CP009273/2422694‑2422987
                                                                                                                                                   |                                                                                                                                                  
gctgctGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATAtt                                                                                                                                                  <  2:1601397/150‑1 (MQ=255)
gctgctGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATAtt                                                                                                                                                  <  1:197670/150‑1 (MQ=255)
 ctgctgGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTc                                                                                                                                                 >  2:128621/1‑150 (MQ=255)
  tgctgGCCCGGCATGAAAAAGGTGCGGCCAACATAGCGGATATTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAAATTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACTTCGATATCCAGATATTcc                                                                                                                                                >  2:1932166/1‑150 (MQ=255)
    ctgGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCa                                                                                                                                              <  2:4622915/150‑1 (MQ=255)
      ggCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCAtt                                                                                                                                            >  1:331459/1‑150 (MQ=255)
      ggCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCAtt                                                                                                                                            <  1:937077/150‑1 (MQ=255)
       gCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTc                                                                                                                                           <  2:4027584/150‑1 (MQ=255)
       gCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTc                                                                                                                                           <  2:3820425/150‑1 (MQ=255)
       gCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCGGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTc                                                                                                                                           <  2:1585081/150‑1 (MQ=255)
          cGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTcgcg                                                                                                                                        <  2:53660/150‑1 (MQ=255)
             cATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGc                                                                                                                                     <  1:2572414/150‑1 (MQ=255)
              aTGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCa                                                                                                                                    >  1:138864/1‑150 (MQ=255)
               tGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCaa                                                                                                                                   >  1:4323248/1‑150 (MQ=255)
                gATAAATGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAAt                                                                                                                                  >  2:3664405/1‑150 (MQ=255)
                gATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAAt                                                                                                                                  <  1:933104/150‑1 (MQ=255)
                gATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAAt                                                                                                                                  >  1:1769923/1‑150 (MQ=255)
                 aTAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAAtt                                                                                                                                 >  1:887893/1‑150 (MQ=255)
                  tAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAAttt                                                                                                                                >  1:2905096/1‑150 (MQ=255)
                  tAAAGGTGCGGCCAACATAGCGGATTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAAttt                                                                                                                                <  1:115110/150‑1 (MQ=255)
                    aaGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTtc                                                                                                                              <  2:3219885/150‑1 (MQ=255)
                     aGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCACTTTtct                                                                                                                             >  1:780704/1‑150 (MQ=255)
                     aGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTtct                                                                                                                             <  1:1417752/150‑1 (MQ=255)
                     aGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTtct                                                                                                                             <  2:3809600/150‑1 (MQ=255)
                     aGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTtct                                                                                                                             >  1:3493740/1‑150 (MQ=255)
                         gCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGcc                                                                                                                         <  2:1009489/150‑1 (MQ=255)
                           ggCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCa                                                                                                                       <  1:16449/150‑1 (MQ=255)
                            gCCAACATAGCGGTTTTTAATGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAg                                                                                                                      <  1:1486731/150‑1 (MQ=255)
                               aaCATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGttt                                                                                                                   <  1:3835295/150‑1 (MQ=255)
                                aCATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTc                                                                                                                  <  1:4563496/150‑1 (MQ=255)
                                aCATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTc                                                                                                                  <  1:296547/150‑1 (MQ=255)
                                       ggTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCa                                                                                                           <  2:2693711/150‑1 (MQ=255)
                                       ggTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCa                                                                                                           <  2:1714899/150‑1 (MQ=255)
                                              aaCGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTcac                                                                                                    >  2:2319833/1‑150 (MQ=255)
                                                 gAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACAcg                                                                                                 >  1:706914/1‑150 (MQ=255)
                                                 gAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACAcg                                                                                                 >  1:2420148/1‑150 (MQ=255)
                                                  aaGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACAcgc                                                                                                >  1:1020930/1‑150 (MQ=255)
                                                    gCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACTCGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACAcgcgc                                                                                              <  2:4023170/150‑1 (MQ=255)
                                                      ccTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACgcgcgc                                                                                            >  1:2041551/1‑150 (MQ=255)
                                                       cTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCt                                                                                           <  2:2489554/150‑1 (MQ=255)
                                                       cTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCt                                                                                           >  1:2706249/1‑150 (MQ=255)
                                                         ggCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGt                                                                                         >  1:2928703/1‑150 (MQ=255)
                                                            ggTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTaaa                                                                                      >  1:904314/1‑150 (MQ=255)
                                                                 ggTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGaa                                                                                 <  1:1590805/150‑1 (MQ=255)
                                                                   tttGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAAt                                                                               <  2:3269326/150‑1 (MQ=255)
                                                                      gCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATCTTCTCACGCGCGCTGTAAACGGAAAtttt                                                                            <  1:3954762/150‑1 (MQ=255)
                                                                        ccAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGt                                                                          >  2:2244028/1‑150 (MQ=255)
                                                                        ccAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGt                                                                          >  2:2761612/1‑150 (MQ=255)
                                                                        ccAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGt                                                                          >  1:2126990/1‑150 (MQ=255)
                                                                           gAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGa                                                                       <  1:1084509/150‑1 (MQ=255)
                                                                           gAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGa                                                                       <  2:109879/150‑1 (MQ=255)
                                                                             aTACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCCTTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCAAATTCACACGCGCGCTGTCAACGGAAATTTTGTCGATa                                                                     >  1:363631/1‑150 (MQ=255)
                                                                             aTACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATa                                                                     >  2:2056019/1‑150 (MQ=255)
                                                                               aCGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATaaa                                                                   >  2:3983736/1‑150 (MQ=255)
                                                                                  aGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGa                                                                >  2:3832250/1‑150 (MQ=255)
                                                                                   gCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAg                                                               <  2:267375/150‑1 (MQ=255)
                                                                                    cAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGt                                                              <  2:3690428/150‑1 (MQ=255)
                                                                                    cAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGt                                                              <  1:303516/150‑1 (MQ=255)
                                                                                      aTTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTcc                                                            >  1:1294244/1‑150 (MQ=255)
                                                                                      aTTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTcc                                                            <  1:2858214/150‑1 (MQ=255)
                                                                                       tttCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCcg                                                           <  2:1854285/150‑1 (MQ=255)
                                                                                       tttCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCcg                                                           <  2:2545339/150‑1 (MQ=255)
                                                                                       tttCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCcg                                                           <  2:4181017/150‑1 (MQ=255)
                                                                                          ccAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGATATTTTGTCGATAAACGAGTCcgggc                                                        <  2:1512794/150‑1 (MQ=255)
                                                                                               gcgATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCa                                                   <  2:2194565/150‑1 (MQ=255)
                                                                                                 gATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCaaa                                                 >  2:3384931/1‑150 (MQ=255)
                                                                                                 gATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCaaa                                                 <  1:30715/150‑1 (MQ=255)
                                                                                                  atatCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAg                                                <  2:728641/150‑1 (MQ=255)
                                                                                                   tatCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGt                                               >  1:4167807/1‑150 (MQ=255)
                                                                                                   tatCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGt                                               >  2:839956/1‑150 (MQ=255)
                                                                                                     tCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTGTGTCGATAAACGAGTCCGGGCGGGCAAAGtat                                             >  2:3398796/1‑150 (MQ=255)
                                                                                                      cacaCGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGtata                                            <  1:4452981/150‑1 (MQ=255)
                                                                                                         acGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTatacat                                         <  2:1669856/150‑1 (MQ=255)
                                                                                                             ggTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTc                                     <  1:3978683/150‑1 (MQ=255)
                                                                                                               tttCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCaa                                   <  1:3439225/150‑1 (MQ=255)
                                                                                                                 tCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAAc                                 <  2:1767474/150‑1 (MQ=255)
                                                                                                                 tCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCAAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAAc                                 <  1:2736529/150‑1 (MQ=255)
                                                                                                                   tGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAg                               >  1:3164019/1‑150 (MQ=255)
                                                                                                                      gatcgggatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCa                            >  2:2472721/1‑150 (MQ=255)
                                                                                                                         cgggatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACgg                         >  2:851562/1‑150 (MQ=255)
                                                                                                                           ggatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGAt                       <  2:3419662/150‑1 (MQ=255)
                                                                                                                           ggatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGAt                       >  1:535469/1‑150 (MQ=255)
                                                                                                                            gatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGAtt                      <  2:2596385/150‑1 (MQ=255)
                                                                                                                             atcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTTTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTg                     >  2:4146152/1‑150 (MQ=255)
                                                                                                                              tcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGc                    <  2:1130173/150‑1 (MQ=255)
                                                                                                                               caccacGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCt                   <  2:3873002/150‑1 (MQ=255)
                                                                                                                                 ccacGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGa                 >  1:4591796/1‑150 (MQ=255)
                                                                                                                                  cacGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGAc                >  2:201577/1‑150 (MQ=255)
                                                                                                                                     gTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCgg             <  2:2357097/150‑1 (MQ=255)
                                                                                                                                       cGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGAt           >  1:1843750/1‑150 (MQ=255)
                                                                                                                                        gATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGTAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGAtt          >  2:1804614/1‑150 (MQ=255)
                                                                                                                                        gATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGAtt          <  1:3338093/150‑1 (MQ=255)
                                                                                                                                        gATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGAtt          >  1:1928700/1‑150 (MQ=255)
                                                                                                                                           atCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTc       >  2:3977376/1‑150 (MQ=255)
                                                                                                                                            tCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCa      <  2:4350351/150‑1 (MQ=255)
                                                                                                                                             ccAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAg     >  2:3303036/1‑150 (MQ=255)
                                                                                                                                             ccAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAg     <  2:3611911/150‑1 (MQ=255)
                                                                                                                                             ccAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAg     >  2:4395455/1‑150 (MQ=255)
                                                                                                                                               aGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAGca   >  2:2785333/1‑150 (MQ=255)
                                                                                                                                                gATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTTTCGATAAACGAGTCCGGGAGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAGcac  >  1:628745/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAGCAC  >  CP009273/2422694‑2422987

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