Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,189,679 G→C G165G (GGC→GGG U712_10970 ← Tryptophan synthase alpha chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,189,6790GC98.7% 256.1 / ‑3.9 79G165G (GGC→GGGU712_10970Tryptophan synthase alpha chain
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base C (49/29);  total (49/30)
Fisher's exact test for biased strand distribution p-value = 3.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

GCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCTAGAGAAGATACACAGTAGACGAAGCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTTACCTCATGGTTTTTGCAATCCGCTTGAAGGCTGGCGCTTTCTTCTAATGGCAGATCCGGAAC  >  CP006881/2189544‑2189810
                                                                                                                                       |                                                                                                                                   
acgGTTTGATATACCGAACCCTACAGCAACCGGAAGAGTGCTGAGATTGTTCTCAGTACGATTGAACGGGTACACGTATGAATTGAACTCATTGCGGAGACCGGTCACAGGAAGAGAAGATACACAGTAGACGAAcccc                                                                                                                                  <  7:289158/138‑1 (MQ=255)
  ggTTTGATATACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGCCACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTc                                                                                                                                >  7:251090/1‑139 (MQ=255)
     ttGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCTAGAGAAGATACACAGTAGACGAAGCCCTCGGc                                                                                                                             <  6:80128/139‑1 (MQ=255)
      tGATATCCCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGc                                                                                                                             <  7:32619/138‑1 (MQ=255)
      tGATATACCGAACCCTACAGCAACCGGACCAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGc                                                                                                                             <  6:282578/138‑1 (MQ=255)
      tGATATACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGc                                                                                                                             <  1:198319/138‑1 (MQ=255)
      tGATATACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGc                                                                                                                             <  8:71465/138‑1 (MQ=255)
      tGATATACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGc                                                                                                                             <  1:231628/138‑1 (MQ=255)
      tGATATACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCt                                                                                                                            >  5:259535/1‑139 (MQ=255)
      tGATAAACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACCGTACGAATGACCGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGc                                                                                                                             <  4:381088/138‑1 (MQ=255)
       gATATACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTCCACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCtt                                                                                                                           >  6:227369/1‑139 (MQ=255)
       gATATACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCtt                                                                                                                           >  7:126009/1‑139 (MQ=255)
       gATATACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCtt                                                                                                                           >  8:285306/1‑139 (MQ=255)
              cGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACCCCGGTCACACCCAGAGAAGATACACAGTAGACGACCCCCTCGGCTTGTTCaa                                                                                                                     >  1:160186/1‑138 (MQ=255)
              cGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCaa                                                                                                                     >  4:391786/1‑138 (MQ=255)
              cGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCaa                                                                                                                     >  1:131987/1‑138 (MQ=255)
                            aCCGGAACAGTGCTGAGATTCTTCACAGTACGAATGCACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCCCACCCAGAGAAGATACACAGTAGACGAACCCCACGGCTAGTTCAATAATGGTTTTCAAAc                                                                                                      >  5:152654/1‑139 (MQ=255)
                               ggAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAACCGGc                                                                                                   >  3:150293/1‑139 (MQ=255)
                                  aCAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCttt                                                                                                >  6:364196/1‑139 (MQ=255)
                                      tGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCt                                                                                            >  4:303512/1‑139 (MQ=255)
                                      tGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCt                                                                                            >  4:90471/1‑139 (MQ=255)
                                      tGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCt                                                                                            >  2:33223/1‑139 (MQ=255)
                                      tGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCt                                                                                            >  7:7724/1‑139 (MQ=255)
                                              atcTTCACACTACGCATGAACGGGGACCCGGATGAATTGAACTCATTGCGGACACCGGGCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGGTCAATAATGGTTTTCAAACGGCTTTCGCGTGTCGgcg                                                                                    <  4:235366/137‑1 (MQ=255)
                                                cTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAACCGGCTTTCGCTTGTCGgcgc                                                                                   >  5:76160/1‑138 (MQ=255)
                                                cTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGgcgc                                                                                   >  5:33938/1‑138 (MQ=255)
                                                       gtacgaaTGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAg                                                                           <  8:126009/133‑1 (MQ=255)
                                                        tacgaaTGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACCCAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCACACGGCTTTCGCTTGTCGGCGCAACTAAAg                                                                           >  1:414863/6‑138 (MQ=255)
                                                        tacgaaTGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGCCGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAg                                                                           >  7:299851/6‑138 (MQ=255)
                                                                aaCGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACg                                                                  >  5:193147/1‑139 (MQ=255)
                                                                aaCGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACg                                                                  >  2:401316/1‑139 (MQ=255)
                                                                 aggggTTCACGGATGAATTGAACTCATTGCGGAAACCGGTCACACCCAGATAAGATACACAGTAGTCGAACCCCTCGGCTTGTTAAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGt                                                                 <  4:400628/137‑1 (MQ=255)
                                                                 aCGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGt                                                                 <  3:90471/139‑1 (MQ=255)
                                                                  cGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGtc                                                                >  6:233864/1‑138 (MQ=255)
                                                                        cacGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTc                                                           <  8:92306/138‑1 (MQ=255)
                                                                         acGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCAt                                                         <  7:285306/139‑1 (MQ=255)
                                                                         acGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCAt                                                         <  6:259535/139‑1 (MQ=255)
                                                                             aTGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATggct                                                     >  2:114206/1‑137 (MQ=255)
                                                                               gAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCATTATTGGGTTTCAAACGGCTTTGGCTTGTCGGCGCTCCTAAAGCATTATACGTCTCCTCATGGcttg                                                   >  2:416212/1‑138 (MQ=255)
                                                                               gAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGt                                                                 >  2:423505/1‑125 (MQ=255)
                                                                               gAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGcttt                                                   >  3:280436/1‑139 (MQ=255)
                                                                               gAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGcttt                                                   >  2:49695/1‑139 (MQ=255)
                                                                               gAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGAAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGt                                                                 <  1:423505/125‑1 (MQ=255)
                                                                                     aaCTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCAt                                                         >  6:321970/1‑127 (MQ=255)
                                                                                     aaCTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCAt                                                         <  5:321970/127‑1 (MQ=255)
                                                                                     aaCTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACatt                                             >  6:126814/1‑137 (MQ=255)
                                                                                     aaCTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTACACGAACCCCTCTGCTTGTTCAATAATGGTTATCAAACGGCTTTCGCTTGTCGGCTCAACTAAAGAAATATACGTCACCTCATGTCTTTTACat                                              >  7:228640/1‑137 (MQ=255)
                                                                                      aCTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTc                                            <  5:364196/139‑1 (MQ=255)
                                                                                       cTCATTGCGGACACCGGTCACCCCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTACAAACGGCTATCGCATGTCGGCGCAACTAAAGAAATATACGTCACCTCATGCCTTTTACATTcc                                           >  3:394411/1‑139 (MQ=255)
                                                                                       cTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATTGTTTTCAAACGGCTTTCGCTTGTCGGTGCAACTAAAGAAATATACGTCACCTCATGGCATTTACATTcc                                           >  5:227986/1‑139 (MQ=255)
                                                                                            tGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTcatctt                                       >  4:342461/1‑133 (MQ=255)
                                                                                            tGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTCCTCtt                                       >  4:374980/1‑138 (MQ=255)
                                                                                                   acCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTCCTCTTGAAGGCTg                               <  3:374980/139‑1 (MQ=255)
                                                                                                   acCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTCCTCTTGAAGGCTg                               <  5:120515/139‑1 (MQ=255)
                                                                                                        tCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTCCTCTTGAAGGCTgttgct                          <  4:280436/139‑6 (MQ=255)
                                                                                                         cacaCCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTCCTCTTGAAGGCTGTTGCtt                         <  2:17270/139‑1 (MQ=255)
                                                                                                          acacCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTCCTCTTGAAGGCTGTTGCttt                        <  1:421212/139‑1 (MQ=255)
                                                                                                            acCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTCCTCTTGAAGGCTGTTGCTttct                      >  4:419154/1‑139 (MQ=255)
                                                                                                                 gagaAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTCCTCTTGCAGGCTGGTGCTTTCTTCTaa                 >  1:330362/1‑139 (MQ=255)
                                                                                                                  agaagaTACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATCAACCTCACCTCATGGCTTTTACATTCATCTTGAAGGGTGTTGCTTTCTTCTAAt                >  3:51700/1‑139 (MQ=255)
                                                                                                                     agaTACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTCCTCTTGAAGGCTGTTGCTTTCTTCTAATGGc             <  8:7724/139‑1 (MQ=255)
                                                                                                                     agaTACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTCCTCTTGAAGGCTGTTGCTTTCTTCTAATGGc             <  2:131987/139‑1 (MQ=255)
                                                                                                                      gaTACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTa                                                                              <  7:147002/73‑1 (MQ=255)
                                                                                                                      gaTACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTa                                                                              >  8:147002/1‑73 (MQ=255)
                                                                                                                      gaTACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTCCTCTTGAAGGCTGTTGCTTTCTTCTAATgg              >  2:268158/1‑137 (MQ=255)
                                                                                                                      gaTACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCATTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCAGGGCTTTTACATTACTCTTGAAGGCTGTTGCTTTCTTCTAATgg              >  5:42311/1‑137 (MQ=255)
                                                                                                                        tACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGt                                                                 >  8:386651/1‑84 (MQ=255)
                                                                                                                        tACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGt                                                                 <  7:386651/84‑1 (MQ=255)
                                                                                                                        tACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTCCTCTTGAAGGCTGTTGCTTTCTTCTAATgg              >  2:328567/1‑135 (MQ=255)
                                                                                                                          cacaGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTCCTCTTGAAGGCTGTTGCTTTCTTCTAATGGCAGAt         <  8:228640/138‑1 (MQ=255)
                                                                                                                          cacaGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTCCTCTTGAAGGCTGTTGCTTTCTTCTAATGGCAGATc        >  5:297120/1‑139 (MQ=255)
                                                                                                                              gTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGt                                                                 >  8:213326/1‑78 (MQ=255)
                                                                                                                              gTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGt                                                                 <  7:213326/78‑1 (MQ=255)
                                                                                                                              gTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTACTCTTGAAGGCTGTTGCTTTCTTCTAATTGCAGATCCgg     >  3:110052/1‑138 (MQ=255)
                                                                                                                              gTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTATCGCTTGTCGGCGCAACTAAAGAAATATACGACACCTCAAGGCTTTTAGATTCCTCTTGAAGGCAGAATCTTACTTCAAATGGCACATCCgg     >  7:215323/1‑138 (MQ=255)
                                                                                                                                  aCGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCAt                                                         >  8:252975/1‑82 (MQ=255)
                                                                                                                                  aCGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCAt                                                         <  7:252975/82‑1 (MQ=255)
                                                                                                                                  aCGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTCCTCTTg                                      <  2:208624/101‑1 (MQ=255)
                                                                                                                                  aCGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTCCTCTTg                                      >  1:208624/1‑101 (MQ=255)
                                                                                                                                  aCGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTCCTCTTGAAGGCTGTTGCTTTCTTCTAATGGCAGATCCGGAAc  >  3:304560/1‑137 (MQ=255)
                                                                                                                                       |                                                                                                                                   
GCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCTAGAGAAGATACACAGTAGACGAAGCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTTACCTCATGGTTTTTGCAATCCGCTTGAAGGCTGGCGCTTTCTTCTAATGGCAGATCCGGAAC  >  CP006881/2189544‑2189810

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: