Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,186,237 A→T G288G (GGT→GGA U712_10955 ← 3‑phosphoshikimate 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,186,2370AT98.9% 272.8 / NA 88G288G (GGT→GGAU712_109553‑phosphoshikimate 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (33/54);  total (34/54)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCCTTAATAACGGTGGTTCCTTCCGCCTGGGTCGCAAGAAGCGCAATGATAGGTATTTCATCAATTAAACGCGGAATGATATCTCCACCAATTTCAACCGCTTTTAGAGATGACGTTTCTATGACCAAATCACCGTAAGGCTCTGCACTGCTGTCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACATCAATGATACCTGTCCGAGTCGGATTTAAACCTACGTTTTTTAATACAATCTTACTGTTTG  >  CP006881/2186104‑2186367
                                                                                                                                     |                                                                                                                                  
cgTCCTTAATAACGGTGGTTCCTTCCGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTccat                                                                                                                                >  2:80363/1‑136 (MQ=255)
cgTCCTTAATAACGGTGGTTCCTTCCGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATa                                                                                                                               <  8:432315/139‑1 (MQ=255)
cgTCCTTAATAACGGTGGTTCCTTCCGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATa                                                                                                                               <  6:132644/139‑1 (MQ=255)
cgTCCTTAATAACGGTGGTTCCTTCCGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATa                                                                                                                               <  7:375752/139‑1 (MQ=255)
cgTCCTTAATAACGGTGGTGCCTTCCGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATa                                                                                                                               <  8:91860/139‑1 (MQ=255)
 gTCCTTAATAACGGTGGTTCCTTCCGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATaa                                                                                                                              <  3:144052/139‑1 (MQ=255)
  tCCTTAATAACGGTGGTTCCTTCCGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAg                                                                                                                             >  4:281100/1‑139 (MQ=255)
  tCCTTAATAACGGTGGTTCCTTCCGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAg                                                                                                                             >  8:116259/1‑139 (MQ=255)
  tCCTTAATAACGGTGGTTCCTTACGCCTGAGTCGCAAGAAGCGCCAAGATAGGCAACACAACAATTAAACCCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGAAGACGTTTCAATAATCAAATCACCATAAg                                                                                                                             >  4:277329/1‑139 (MQ=255)
    cTTAATAACGGTGGTTCCTTCCGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAgg                                                                                                                            <  3:409039/138‑1 (MQ=255)
       aataaCGGTGGTTCCTTCCGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGct                                                                                                                          <  8:57827/137‑1 (MQ=255)
           aCGGTGGTTCCTTCCGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTTAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCa                                                                                                                     <  2:99920/138‑1 (MQ=255)
           aCGGTGGTTCCTTCCGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCa                                                                                                                     <  3:277329/138‑1 (MQ=255)
           aCGGTGGTTCCTTCCGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCa                                                                                                                     <  5:419589/138‑1 (MQ=255)
           aCGGTGGTTCCTTCCGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCa                                                                                                                     <  3:375009/138‑1 (MQ=255)
           aCGGTGGTTCCTTCCGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCa                                                                                                                     <  7:49750/138‑1 (MQ=255)
            cGGTGGTTCCTTCCGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATa                                                                                                                               <  4:341653/127‑1 (MQ=255)
            cGGTGGTTCCTTCCGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATa                                                                                                                               >  3:341653/1‑127 (MQ=255)
              gtggtTCCTTCCGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCAccgc                                                                                                                 >  6:137222/1‑139 (MQ=255)
               tggtTCCTTCCGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCAccgct                                                                                                                <  7:116259/139‑1 (MQ=255)
               tggtTCCTTCCGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCAccgct                                                                                                                <  3:407960/139‑1 (MQ=255)
                 gtTCCTTCCGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCAccgctat                                                                                                              >  7:74636/1‑137 (MQ=255)
                 gtTCCTTCCGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCAccgctat                                                                                                              >  4:7244/1‑137 (MQ=255)
                 gtTCCTTCCGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCAccgctat                                                                                                              >  7:247706/1‑137 (MQ=255)
                 gtTCCTTCCGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCAccgctat                                                                                                              >  3:63361/1‑137 (MQ=255)
                       tccGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCa                                                                                                                     <  2:8066/126‑1 (MQ=255)
                       tccGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCa                                                                                                                     >  1:8066/1‑126 (MQ=255)
                       tccGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATcag                                                                                                           >  6:48266/1‑136 (MQ=255)
                        ccGCCTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATcagca                                                                                                         <  3:206543/137‑1 (MQ=255)
                           ccTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATa                                                                                                                               >  5:354100/1‑112 (MQ=255)
                           ccTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATa                                                                                                                               <  6:354100/112‑1 (MQ=255)
                            cTGAGTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGt                                                                                                   <  5:438726/139‑1 (MQ=255)
                               agTCGCAAGAAGCGCGATGATAGGGATCTCTTCAATTAAACGCGGTATGATATCTGCTCCGATTTCATCTGCCTTTAGAGATGACGTTTCTATATTCTTTTCTCCATAAGGGTCTGCTCCGCTATCAGCAGATGGTTTg                                                                                                >  5:359438/2‑139 (MQ=255)
                               agTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGttt                                                                                                 >  1:199016/2‑138 (MQ=255)
                               agTCGCAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTGAGAGAACACGTTTCTATAATCAAATCTCCATATGGCTCTGCACCGCTATCAGCAGATGGatt                                                                                                 >  5:365445/2‑138 (MQ=255)
                               agTCGCAAGAAGCGCGAGGATAGGGCTCTCATCAATTAAACGCGGAATCATATCTCCTCCGATTTCAACTGCCTTTAGAGATCACGTTTCTATAATCAAATCTCCATAATGCTCTGCACCGCTATCACAAGATGGttt                                                                                                 >  2:366445/2‑138 (MQ=255)
                                    cAAGAAGCGCGATGATAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTc                                                                                           <  7:57636/139‑1 (MQ=255)
                                                  aTAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGccccc                                                                              <  7:405517/138‑1 (MQ=255)
                                                  aTAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCGCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCa                                                                             <  2:289312/139‑1 (MQ=255)
                                                   tAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCAccgctat                                                                                                              <  3:366628/105‑3 (MQ=255)
                                                   tAGGGATCTCATCAATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCAccgctat                                                                                                              >  4:366628/1‑103 (MQ=255)
                                                          ctcatcaATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTg                                                                      >  3:132650/2‑138 (MQ=255)
                                                          ctcatcaATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTGTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTg                                                                      >  8:257366/2‑138 (MQ=255)
                                                            catcaATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATcagca                                                                                                         >  8:128292/1‑101 (MQ=255)
                                                            catcaATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATcagca                                                                                                         <  7:128292/101‑1 (MQ=255)
                                                            catcaATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAg                                                                   >  1:2527/1‑139 (MQ=255)
                                                               caATTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGAc                                                                <  1:288657/139‑1 (MQ=255)
                                                                aaTTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGAc                                                                <  4:41364/138‑1 (MQ=255)
                                                                aaTTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGAc                                                                <  2:2527/138‑1 (MQ=255)
                                                                aaTTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGAc                                                                <  5:137222/138‑1 (MQ=255)
                                                                aaTTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGAc                                                                <  3:7244/138‑1 (MQ=255)
                                                                aaTTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGAc                                                                <  4:20333/138‑1 (MQ=255)
                                                                aaTTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACa                                                               <  3:91410/139‑1 (MQ=255)
                                                                aaTTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACa                                                               <  2:199016/139‑1 (MQ=255)
                                                                aaTTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTAAAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGAc                                                                <  8:427272/138‑1 (MQ=255)
                                                                 aTTAAACGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACAt                                                              >  4:351206/1‑139 (MQ=255)
                                                                     aaCGCGGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACATCAAt                                                          <  4:63361/139‑1 (MQ=255)
                                                                       cgcgGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTg                                                                                   >  3:11333/1‑112 (MQ=255)
                                                                       cgcgGAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTg                                                                                   <  4:11333/112‑1 (MQ=255)
                                                                          ggAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGAc                                                                >  2:183402/1‑128 (MQ=255)
                                                                          ggAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGAc                                                                <  1:183402/128‑1 (MQ=255)
                                                                          ggAATGATATCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACATCAATAATAc                                                     <  4:132650/139‑1 (MQ=255)
                                                                                atatCTCCTCCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACATCAATAATACCTGTcc                                               <  1:80363/139‑1 (MQ=255)
                                                                                    ctcctcCGATTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACATCAATAATACCTGTCCGAGt                                           >  8:234980/1‑139 (MQ=255)
                                                                                      cctccgaTTTCATCTGCCTTTAGAGATGACGTTACTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGAAGCTTTGATTTCACGTTTTGCCCCCATGTTTTGAGGGCCCTCAATATTACCTGTCCTAGTCg                                         >  8:104782/7‑139 (MQ=255)
                                                                                      cctccgaTTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACATCAATAATACCTGTCCGAGTCg                                         >  1:18040/7‑139 (MQ=255)
                                                                                       ctgcgaTTTCTACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACAGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACATCAATAATACCTGTCCGAGt                                           <  1:366445/131‑1 (MQ=255)
                                                                                       ctccgaTTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACATCAATAATACCTGTCCGAGt                                           <  3:281100/131‑1 (MQ=255)
                                                                                            aTTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCa                                                                                                                     <  3:260694/57‑1 (MQ=255)
                                                                                            aTTTCAACTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCa                                                                                                                     >  4:260694/1‑57 (MQ=255)
                                                                                                 aaCTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGAc                                                                >  5:300602/1‑105 (MQ=255)
                                                                                                 aaCTGCCTTTAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGAc                                                                <  6:300602/105‑1 (MQ=255)
                                                                                                          tAGAGATGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACATCAATAATACCTGTCCGAGTCGGATTTAAACCTACGttttt                      <  2:385549/138‑1 (MQ=255)
                                                                                                            gagaTGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACATCAATAATACCTGTCCGAGTCGGATTTAAACCTACGTTTTTCaa                   <  8:247706/139‑1 (MQ=255)
                                                                                                            gagaTGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGGTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACATCAATAATACCTGTCCGAGTCGGATTTAAACCTACGTTTTTCaa                   <  8:74636/139‑1 (MQ=255)
                                                                                                              gaTGACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCATGTTTTGAAGGACATCAATAATACCTGTCCGAGTCGGATTTATACCTACGTTTTTCAATa                 <  3:169360/138‑1 (MQ=255)
                                                                                                                 gACGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACATCAATAATACCTGTCCGAGTCGGATTTAAACCTACGTTTTTCAATACaa              >  5:184697/1‑139 (MQ=255)
                                                                                                                  aCGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACATCaa                                                           <  8:380352/93‑1 (MQ=255)
                                                                                                                  aCGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACATCaa                                                           >  7:380352/1‑93 (MQ=255)
                                                                                                                  aCGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACATCATTAATACCTGTCCGAGTCGGATTTAAACCTACGTTTTTCAATACAAt             >  8:332967/1‑139 (MQ=255)
                                                                                                                  aCGTTTCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACATCAATAATACCTGTCCGAGTCGGATTTAAACCTACGTTTTTCAATACaa              <  6:365445/138‑1 (MQ=255)
                                                                                                                  aCGTTTCAATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTATGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACATCAATAATACCTGTCCGAGTCGGATTTAAACCTACGTTTTTCAATACaa              <  6:359438/138‑1 (MQ=255)
                                                                                                                       tCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACATCAATAATACCTGTCCGAGTCGGATTTAAACCTACGTTTTTCAATACAAttctg         <  7:257366/138‑5 (MQ=255)
                                                                                                                       tCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACATCAATAATACCTGTCCGAGTCGGATTTAAACCTACGTTTTTCAATACAAttctg         <  6:299554/138‑5 (MQ=255)
                                                                                                                       tCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACATCAATAATACCTGTCCGAGTCGGATTTAAACCTACGTTTTTCAATACAAttctg         <  5:48266/138‑5 (MQ=255)
                                                                                                                       tCTATAATCAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACATCAATAATACCTGTCCGAGTCGGATTTAAACCTACGTTTTTCAATACAAttctg         <  2:43648/138‑5 (MQ=255)
                                                                                                                       tCTATAATCAAATCTCCATAAGGCTCTGCACCGCTAGCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACATCAATAATACCTGTCCGAGTCGGATTTAAACCTACGTTTTTCAATACAAttctg         <  7:67996/138‑5 (MQ=255)
                                                                                                                              tcAAATCTCCATAAGGCTCTGCACCGCTATCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACATCAATAATACCTGTCCGAGTCGGATTTAAACCTACGTTTTTCAATACAATTCTGCTGTTTg  >  3:421283/2‑138 (MQ=255)
                                                                                                                                     |                                                                                                                                  
CGTCCTTAATAACGGTGGTTCCTTCCGCCTGGGTCGCAAGAAGCGCAATGATAGGTATTTCATCAATTAAACGCGGAATGATATCTCCACCAATTTCAACCGCTTTTAGAGATGACGTTTCTATGACCAAATCACCGTAAGGCTCTGCACTGCTGTCAGCAGATGGTTTGATTTCAAGTTTTGCCCCCATGTTTTGAAGGACATCAATGATACCTGTCCGAGTCGGATTTAAACCTACGTTTTTTAATACAATCTTACTGTTTG  >  CP006881/2186104‑2186367

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: