Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,187,748 T→C V160V (GTA→GTG U712_10960 ← Prephenate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,187,7480TC98.6% 224.6 / ‑5.2 71V160V (GTA→GTGU712_10960Prephenate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (44/26);  total (45/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTGATGAACTAGGCTTGCTGCGACAATATGCGGAAAATGGCTGATCACGCTTGTAACACCATCATGCTCTTCCGGCGACATTTCCACAAAATGGGCATTTGTCCCCTTCAGCAGGTTTTTCAATTGTTCTACCGCTTGTCTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAATTCCTTCGCAGCGGCCACTCCTGATTTATGTGAACCTGCCATCGGATGCCCCCCGACAAATTGGTATCGGC  >  CP006881/2187617‑2187881
                                                                                                                                   |                                                                                                                                     
gTTTGGTGAACAAGGCTAGCTGCTACAATATGCGGAAAATGACTGATTACGCTTGTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCtt                                                                                                                                >  8:155599/1‑139 (MQ=255)
        aaCAAGGCTAGCTGCTACAATATGCGGAAAATGACTGATTACGCTTGTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTcc                                                                                                                        <  1:67961/139‑1 (MQ=255)
          caaGGCTAGCTGCTACAATATGCGGAAAATGACTGATTACGCTTGTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGt                                                                                                                      >  8:50425/3‑139 (MQ=255)
             ggCTAGCTGCTACAATATGCGGAAAATGACTGATTACGCTTGTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCAGTCGTCCCCTGCAGCAGGTTTTTTAACTGTTCCCCAGCTTGTGTGTCCGtttt                                                                                                                   >  6:35394/1‑139 (MQ=255)
                tagctgctACAATATGCGGAAAATGACTGATTACGCTTGTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTctg                                                                                                                >  8:444560/3‑139 (MQ=255)
                  gctgctACAATATGCGGAAAATGACTGATTACGCTTGTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTctggc                                                                                                              >  3:385252/1‑139 (MQ=255)
                       taCAATATGCGGAAAATGACTGATTACGCTTGTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTcc                                                                                                                        >  3:346744/2‑124 (MQ=255)
                       taCAATATGCGGAAAATGACTGATTACGCTTGTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTcc                                                                                                                        <  4:346744/123‑1 (MQ=255)
                       taCAATATGCGGAAAATGACTGATTACGCTTGTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTctggcctggc                                                                                                         >  1:176303/2‑139 (MQ=255)
                       taCAATATGCGGAAAATGACTGATTACGCTTGTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTctggcctggc                                                                                                         >  6:132503/2‑139 (MQ=255)
                       taCAATATGCGGAAAATGACTGATTACGCTTGTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTctggcctggc                                                                                                         >  8:188766/2‑139 (MQ=255)
                          aaTATGCGGAAAATGACTGATTACGCTTGTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGtt                                                                                                      <  7:188766/139‑1 (MQ=255)
                           atatGCGGAAAATGACTGATTACGCTTGTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAACTTGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTa                                                                                                     <  8:34575/139‑1 (MQ=255)
                           atatGCGGAAAATGACTGATTACGCTTGTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTa                                                                                                     <  5:120770/139‑1 (MQ=255)
                            tatGCGGAAAATGACTGATTACGCTTGTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTa                                                                                                     >  5:421112/1‑138 (MQ=255)
                            tatGCGGAAAATGACTGATTACGCTTGTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTGTTCTGGCCTGGCGTTa                                                                                                     >  4:344610/1‑138 (MQ=255)
                            tatGCGGAAAATGACTGATTACGATTGTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGTTTATTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTa                                                                                                     >  5:140665/1‑138 (MQ=255)
                              tGCGGAAAATGACTGATTACGCTTGTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTtgtt                                                                                                                             <  3:385747/112‑2 (MQ=255)
                              tGCGGAAAATGACTGATTACGCTTGTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTtgtt                                                                                                                             >  4:385747/1‑111 (MQ=255)
                              tGCGGAAAATGACTGATTACGCTTGTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGtttt                                                                                                                   >  3:234468/1‑122 (MQ=255)
                              tGCGGAAAATGACTGATTACGCTTGTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGtttt                                                                                                                   <  4:234468/122‑1 (MQ=255)
                                cGGAAAATGACTGATTACGCTTGTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAata                                                                                                >  4:178229/1‑139 (MQ=255)
                                cGGAAAATGACTGATTACGCTTGTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTGAAAata                                                                                                >  7:221263/1‑139 (MQ=255)
                                                 cgcgTGTAACGCCATCAAGCTCCTCTCGCGACATATCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTATAACTGGTCCACAGCTTGTTAGTCCGTGTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCaa                                                                                <  6:92781/138‑1 (MQ=255)
                                                  gCTTGTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAAAGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCaaa                                                                               >  7:402034/1‑138 (MQ=255)
                                                      gTAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAgg                                                                           <  5:298915/138‑1 (MQ=255)
                                                       tAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGaa                                                                         >  3:251910/1‑139 (MQ=255)
                                                       tAACGCCATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGaa                                                                         >  2:433319/1‑139 (MQ=255)
                                                            ccATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGaa                                                                         >  7:87917/1‑134 (MQ=255)
                                                            ccATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGaa                                                                         <  8:87917/134‑1 (MQ=255)
                                                            ccATCATGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCt                                                                    >  1:229146/1‑139 (MQ=255)
                                                              atcatGCTCCTCTGGCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTc                                                                  >  7:338052/1‑139 (MQ=255)
                                                               tcatGCTCTTCCGGCGACATTTCCACAAAATGGGCATTTGTCCCCTTCAGCAGGTTTTTCAATTGTTCTACCGCTTGTCTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAATTCCTTCg                                                                 >  5:142899/1‑139 (MQ=255)
                                                                      cctctggCGACATTTCCTCAAAATGGGCATTCGTCCACTTCAGCATGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTTAAATATAAAATGCATTTTCAAACTGGAACTCCTTCGCAGCGGc                                                          <  6:239690/134‑1 (MQ=255)
                                                                       ctctggCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGcc                                                         >  3:241956/5‑139 (MQ=255)
                                                                       ctctggCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGc                                                          <  5:132503/134‑1 (MQ=255)
                                                                       ctctggCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGc                                                          <  3:162875/134‑1 (MQ=255)
                                                                       ctctggCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGc                                                          <  7:155599/134‑1 (MQ=255)
                                                                       ctctggCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGc                                                          <  3:178229/134‑1 (MQ=255)
                                                                       ctctggCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGc                                                          <  4:385252/134‑1 (MQ=255)
                                                                       ctctggCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTCAAACAGGAACTCCTTCGCAGCGGccc                                                        >  1:236594/5‑138 (MQ=255)
                                                                       ctctggCGACATTTCCACAAAATGGGCAATCGTCCCCTTCAGCAGGTTTTTTAACTGATCCACAGCTTGTATGTCCCTTATCTGGCCTGGCGTTATAATATAAAATGCATATTCAAACAGGAACTCCTTCTCATCGccc                                                         >  2:205747/5‑139 (MQ=255)
                                                                         ctggCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAAtata                                                                                              <  2:442016/98‑1 (MQ=255)
                                                                         ctggCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAAtata                                                                                              >  1:442016/3‑100 (MQ=255)
                                                                         ctggCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGCCAc                                                       >  3:115761/3‑139 (MQ=255)
                                                                         ctggCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCACCGGCCAc                                                       >  3:37991/3‑139 (MQ=255)
                                                                            gCGACATTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGCCACTc                                                     <  4:37991/138‑1 (MQ=255)
                                                                                 aTTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTa                                                                                                     <  5:192808/85‑1 (MQ=255)
                                                                                 aTTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTa                                                                                                     >  6:192808/1‑85 (MQ=255)
                                                                                 aTTTCCACAAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGCCACTCCTGAtt                                               >  4:141006/1‑139 (MQ=255)
                                                                                      cacaAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCa                                                               <  8:437532/118‑1 (MQ=255)
                                                                                      cacaAAATGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCa                                                               >  7:437532/1‑118 (MQ=255)
                                                                                           aaTGGGCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGCCACTCCTGATTTATGTGAAcc                                     >  8:430800/1‑139 (MQ=255)
                                                                                                gCATTCGTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGCCACTCCTGATTTATGTGAACCCGCCa                                >  7:239858/1‑139 (MQ=255)
                                                                                                   ttcgTCCCCTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGCCACTCCTGATTTATGTGAACCCGCCATCg                             >  8:440539/4‑139 (MQ=255)
                                                                                                           cTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGc                                                          >  5:45686/1‑102 (MQ=255)
                                                                                                           cTTCAGCAGGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGc                                                          <  6:45686/102‑1 (MQ=255)
                                                                                                              cagcagGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTTGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGCCACTCCTGATTTATGTGAACCCGCCATCGGATGccctcc                   <  8:235192/138‑1 (MQ=255)
                                                                                                                 cagGTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGCCACTCCTGATTTATGTGAACCCGCCATCGGATGCCCTCCGAc                >  2:107460/1‑138 (MQ=255)
                                                                                                                    gTTTTTTAACTGTTCCACAGCTTGTTTTTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGc                                                          <  5:51108/93‑1 (MQ=255)
                                                                                                                    gTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGc                                                          >  6:51108/1‑93 (MQ=255)
                                                                                                                    gTTTTTTAACTGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGCCACTCCTGATTTATGTGAACCCGCCATCGGATGCCCTCCGACAAAt            >  3:48194/1‑139 (MQ=255)
                                                                                                                           aactGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGCCACTCCTGATTTATGTGAACCCGCCATCGGATGCCCTCCGACAAATTGATag      <  7:444560/135‑2 (MQ=255)
                                                                                                                              tGTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGCCTCTCCTGATTTATGTGAACCCGCCATCGGAGGCCCTCCGACAAATTGATAGCGGc  <  6:140665/139‑1 (MQ=255)
                                                                                                                               gTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGCCACTCCTGATTTATGTGAACCCGCCATCGGATGCCCTCCGACAAATTgat        <  1:8335/132‑3 (MQ=255)
                                                                                                                               gTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGCCACTCCTGATTTATGTGAACCCGCCATCGGATGCCCTCCGACAAATTgat        >  2:8335/1‑130 (MQ=255)
                                                                                                                               gTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGCCACTCCTGATTTATGTGAACCCGCCATCGGATGCCCTCCGACAAATTGATAGCGGc  >  7:69353/1‑138 (MQ=255)
                                                                                                                               gTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGCCACTCCTGATTTATGTGAACCCGCCATCGGATGCCCTCCGACAAATTGATAGCGGc  >  7:126585/1‑138 (MQ=255)
                                                                                                                               gTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGCCACTCCTGATTTATGTGAACCCGCCATCGGATGCCCTCCGACAAATTGATAGCGGc  >  1:401862/1‑138 (MQ=255)
                                                                                                                               gTTCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGCCACTCCTGAATTATGTGAACCCGCCATCGGATGCCCTCCCACAAATTGATAGCgtc  >  6:14874/1‑136 (MQ=255)
                                                                                                                                ttCCACAGCTTGTTTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAACTCCTTCGCAGCGGCCACTCCTGATTTATGTGAACCCGCCATCGGATGCCCTCCGACAAATTGATAGCGGc  >  5:228584/1‑137 (MQ=255)
                                                                                                                                   |                                                                                                                                     
GTTTGATGAACTAGGCTTGCTGCGACAATATGCGGAAAATGGCTGATCACGCTTGTAACACCATCATGCTCTTCCGGCGACATTTCCACAAAATGGGCATTTGTCCCCTTCAGCAGGTTTTTCAATTGTTCTACCGCTTGTCTGTCCGTTTTCTGGCCTGGCGTTAAAATATAAAATGCATTTTCAAACAGGAATTCCTTCGCAGCGGCCACTCCTGATTTATGTGAACCTGCCATCGGATGCCCCCCGACAAATTGGTATCGGC  >  CP006881/2187617‑2187881

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: