Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,724,587 C→T 100% G80G (GGG→GGA U712_19200 ← SacPA operon antiterminator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,724,5870CT99.9% 225.5 / ‑7.5 71G80G (GGG→GGAU712_19200SacPA operon antiterminator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/41);  minor base G (2/0);  total (30/41)
Fisher's exact test for biased strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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.

TTTTTTATAGAGCGCTTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCAGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCCCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCAATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTAAACTTTTCATTCTCGTCGTGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGA  >  CP006881/3724452‑3724717
                                                                                                                                       |                                                                                                                                  
ttttttATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCGCTCCt                                                                                                                                 >  5:153379/1‑139 (MQ=255)
  ttttATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTc                                                                                                                               >  1:166186/1‑139 (MQ=255)
  ttttATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTc                                                                                                                               >  6:78537/1‑139 (MQ=255)
  ttttATAGAGCGACTTTATCTCATGCAGCAATTTATTCTTCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCtt                                                                                                                                <  3:220840/138‑1 (MQ=255)
     tataGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCg                                                                                                                             <  2:78298/138‑1 (MQ=255)
      ataGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAATCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGc                                                                                                                            <  8:225398/138‑1 (MQ=255)
      ataGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGc                                                                                                                            <  2:52481/138‑1 (MQ=255)
       taGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCa                                                                                                                           <  2:209412/138‑1 (MQ=255)
          agctcctttATCTCCTGCTGCAATTTATTCTGCACATGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAg                                                                                                                        <  7:133997/132‑1 (MQ=255)
          agCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGc                                                                                                                       <  6:131128/139‑1 (MQ=255)
            cgcCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGCGCTCTCCTTCCGCATAACt                                                                                                                      >  3:218472/1‑138 (MQ=255)
            cgcCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCt                                                                                                                                 >  6:199978/1‑127 (MQ=255)
            cgcCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCt                                                                                                                                 <  5:199978/127‑1 (MQ=255)
            cgcCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCa                                                                                                                           >  6:234960/1‑133 (MQ=255)
            cgcCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCa                                                                                                                           <  5:234960/133‑1 (MQ=255)
            cgcCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCtg                                                                                                                     <  3:288217/139‑1 (MQ=255)
            cgcCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCt                                                                                                                      >  6:216827/1‑138 (MQ=255)
                  tATCTCATGCCGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGGCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTgatgat                                                                                                                <  6:16579/138‑1 (MQ=255)
                  tATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTgatgat                                                                                                                <  8:163533/138‑1 (MQ=255)
                        aTGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGAGCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCa                                                                                                         <  4:141074/139‑1 (MQ=255)
                         tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCt                                                                                                           >  7:254244/1‑136 (MQ=255)
                         tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCt                                                                                                           <  8:254244/136‑1 (MQ=255)
                         tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAg                                                                                                        >  2:35534/1‑139 (MQ=255)
                         tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAg                                                                                                        >  2:233157/1‑139 (MQ=255)
                         tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAg                                                                                                        <  6:153379/139‑1 (MQ=255)
                         tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCAGCAAAGCTGATGATATACTCAg                                                                                                        >  2:34391/1‑139 (MQ=255)
                         tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTACGCATAGCTGATGATATCcacac                                                                                                        >  8:50228/1‑138 (MQ=255)
                         tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCGTCCGCATAGCTGATGAAATCCTCAg                                                                                                        >  6:173664/1‑139 (MQ=255)
                            agcaATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTa                                                                                                     <  7:170057/139‑1 (MQ=255)
                            agcaATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTa                                                                                                     <  1:155765/139‑1 (MQ=255)
                                          cacaAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGCGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGctcc                                                                                       >  1:257538/1‑139 (MQ=255)
                                          cacaAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTgatgat                                                                                                                >  7:52474/1‑114 (MQ=255)
                                          cacaAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTgatgat                                                                                                                <  8:52474/114‑1 (MQ=255)
                                          cacaAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGc                                                                                          <  1:34391/136‑1 (MQ=255)
                                          cacaAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGc                                                                                          <  1:233157/136‑1 (MQ=255)
                                             aaGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCGGGCCTTCCGCATAGCTGATGAGATCCTCAGCGATTTATATATGCtccgcc                                                                                    >  5:181356/1‑139 (MQ=255)
                                                cagcCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGc                                                                                 <  1:162173/139‑1 (MQ=255)
                                                 agcCCATTTGGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCGCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCa                                                                                <  4:36383/139‑1 (MQ=255)
                                                  gcCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAg                                                                               >  8:29778/1‑139 (MQ=255)
                                                      aTTTTGAATCCTTTCGATCGCAAAGGACAATTGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGtt                                                                           <  7:32842/139‑1 (MQ=255)
                                                       ttttGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCGACTTCCCCATAGCTGCAGAAATCCTCAGCAAATTATATAAGCTCCTCCGGCAGTGgtt                                                                          >  8:32842/1‑139 (MQ=255)
                                                         ttGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGc                                                                        >  3:91942/1‑139 (MQ=255)
                                                             aTCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGa                                                                    >  7:277600/1‑139 (MQ=255)
                                                             aTCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGa                                                                    >  7:277597/1‑139 (MQ=255)
                                                                  ttCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTGTCCGTCCGCATAGCTGATGATATCCTCAGCGATGTCTATATTCTCCTCCGGCAGGTTTTTAAGGATTTGt                                                               >  3:118492/1‑139 (MQ=255)
                                                                     gATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGtttg                                                            <  2:166186/139‑2 (MQ=255)
                                                                           aaaGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACttt                                                      >  6:145936/1‑139 (MQ=255)
                                                                                       gTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGc                                                                                                                       >  8:199839/1‑62 (MQ=255)
                                                                                       gTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGc                                                                                                                       <  7:199839/62‑1 (MQ=255)
                                                                                         cggAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCtcgtcgcg                                        <  5:78537/137‑3 (MQ=255)
                                                                                           gAAAGCGCGATATGAATGTTATCTCTCAGCGGCGCCGAGAGCTCTCCTTCCTCATAGCTTATGATATCCTCAGCTATTTCTATATGCTCCTCCTGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCa                                      <  8:230342/139‑1 (MQ=255)
                                                                                           gAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCa                                      <  7:29778/139‑1 (MQ=255)
                                                                                           gAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCa                                      <  6:181356/139‑1 (MQ=255)
                                                                                                cgcgATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACa                                 <  6:227535/139‑1 (MQ=255)
                                                                                                   gATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAg                              <  1:35534/139‑1 (MQ=255)
                                                                                                   gATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAg                              <  1:221490/139‑1 (MQ=255)
                                                                                                             tgATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCa                                                                                                                           >  8:287504/1‑36 (MQ=255)
                                                                                                             tgATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCa                                                                                                                           <  7:287504/36‑1 (MQ=255)
                                                                                                             tgATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTg                                                                         >  4:131170/1‑86 (MQ=255)
                                                                                                             tgAGCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTg                                                                         <  3:131170/86‑1 (MQ=255)
                                                                                                               aTCGCTCCGCGGCGCCGCGCGCGGTCCTTCCGCATAGCTGCTGATATCCTCAGCTATGTCTTTATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTGGAACTTTTCATTCTCGTCGCGCACGCCAAAGGTCTTTTCCACt                  >  5:282652/1‑139 (MQ=255)
                                                                                                               aTCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGc                                                                                                                       >  7:272864/1‑38 (MQ=255)
                                                                                                               aTCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGc                                                                                                                       <  8:272864/38‑1 (MQ=255)
                                                                                                               aTCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACt                  <  7:94893/139‑1 (MQ=255)
                                                                                                               aTCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACt                  <  6:56214/139‑1 (MQ=255)
                                                                                                               aTCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACt                  >  1:177600/1‑139 (MQ=255)
                                                                                                                      aGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCa           <  2:177600/139‑1 (MQ=255)
                                                                                                                      aGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCa           <  2:257538/139‑1 (MQ=255)
                                                                                                                      aGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTATCGTCGCGCACGACAAAGATCTTTTACACTTTAttaa           >  4:133530/1‑137 (MQ=255)
                                                                                                                               gcgAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGa  <  4:118492/139‑1 (MQ=255)
                                                                                                                               gcgAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGa  <  3:139685/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                  
TTTTTTATAGAGCGCTTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCAGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCCCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCAATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTAAACTTTTCATTCTCGTCGTGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGA  >  CP006881/3724452‑3724717

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