Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,043,666 T→A 5.6% intergenic (+18/+16) gmk → / ← PP_5297 deoxyguanylate kinase/guanylate kinase/amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,043,6660TA 5.6% 217.9 / 2.7 71intergenic (+18/+16)gmk/PP_5297deoxyguanylate kinase/guanylate kinase/amino acid transporter
Reads supporting (aligned to +/- strand):  ref base T (30/37);  new base A (2/2);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

TACGTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGAA  >  NC_002947/6043520‑6043807
                                                                                                                                                  |                                                                                                                                             
tACGTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTg                                                                                                                                             <  6:74235/149‑1 (MQ=255)
  cGTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTg                                                                                                                                             <  8:202937/147‑1 (MQ=255)
  cGTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCa                                                                                                                                           <  2:214619/149‑1 (MQ=255)
   gTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAg                                                                                                                                          >  3:174247/1‑149 (MQ=255)
    tcatcatcaATGATGATTTCGACACTGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGc                                                                                                                                         <  2:285491/149‑1 (MQ=255)
      atcatcaATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGccc                                                                                                                                       <  2:164297/149‑1 (MQ=255)
            aaTGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATg                                                                                                                                  <  4:172451/148‑1 (MQ=255)
              tgatgaTTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGtt                                                                                                                               <  1:333592/149‑1 (MQ=255)
              tgatgaTTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGtt                                                                                                                               <  3:170710/149‑1 (MQ=255)
               gatgatTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGtt                                                                                                                               >  1:302727/1‑148 (MQ=255)
                  gatTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGt                                                                                                                           >  6:36378/1‑149 (MQ=255)
                   atTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCAATGCAGCCCCAATGGTTTAGTc                                                                                                                          >  4:334395/1‑149 (MQ=255)
                        gACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCa                                                                                                                     <  8:250150/149‑1 (MQ=255)
                        gACACGGCGCTAGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTg                                                                                                                       >  5:187986/1‑147 (MQ=255)
                           acGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCAATGCAGCCCCatcgat                                                                                                                                <  7:13372/135‑6 (MQ=255)
                           acGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCAATGCAGCCCCatcgat                                                                                                                                >  8:13372/1‑130 (MQ=255)
                             ggCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCgg                                                                                                                 >  5:34447/1‑148 (MQ=255)
                                gcTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAg                                                                                                             <  7:67925/149‑1 (MQ=255)
                                  tGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAgg                                                                                                            >  5:251779/1‑148 (MQ=255)
                                    gCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTa                                                                                                          <  3:307446/148‑1 (MQ=255)
                                            gAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGa                                                                                                 >  1:142931/1‑149 (MQ=255)
                                                gCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTa                                                                                                                             >  8:79938/1‑117 (MQ=255)
                                                gCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCAATGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCa                                                                                             <  3:334395/149‑1 (MQ=255)
                                                  aGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAgg                                                                                           <  6:317542/149‑1 (MQ=255)
                                                       tCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGc                                                                                      >  2:234842/1‑149 (MQ=255)
                                                          gCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGcccc                                                                                   >  8:242014/1‑149 (MQ=255)
                                                            tCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTa                                                                                                                             >  3:205423/1‑105 (MQ=255)
                                                            tCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTa                                                                                                                             <  4:205423/105‑1 (MQ=255)
                                                              gAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGa                                                                                 <  6:187986/147‑1 (MQ=255)
                                                              gAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAg                                                                                <  6:303634/148‑1 (MQ=255)
                                                              gAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAg                                                                                >  8:295865/1‑148 (MQ=255)
                                                              gAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGGCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAg                                                                                <  5:36378/148‑1 (MQ=255)
                                                                aTCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAgg                                                                               >  8:245028/1‑147 (MQ=255)
                                                                aTCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGaa                                                                             <  4:174247/149‑1 (MQ=255)
                                                                   ggCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGGAGt                                                                           <  3:234657/148‑1 (MQ=255)
                                                                   ggCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTa                                                                          >  3:151267/1‑149 (MQ=255)
                                                                    gctgctgctGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTAc                                                                         <  6:34447/149‑1 (MQ=255)
                                                                        ctgctgAAGAAGCAGCAGCAGCGTCCTGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGc                                                                       <  7:242170/147‑1 (MQ=255)
                                                                        ctgctgAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCg                                                                     <  2:302727/149‑1 (MQ=255)
                                                                           ctgAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGgcg                                                                  >  5:323025/1‑149 (MQ=255)
                                                                            tgAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGgcg                                                                  >  4:207412/1‑148 (MQ=255)
                                                                            tgAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGgcg                                                                  >  3:196299/1‑148 (MQ=255)
                                                                              aagaagCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCa                                                                                             <  6:92160/119‑1 (MQ=255)
                                                                              aagaagCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCa                                                                                             >  5:92160/1‑119 (MQ=255)
                                                                              aagaagCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGcgca                                                                <  6:251779/148‑1 (MQ=255)
                                                                                   gcagcagcagcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATg                                                          >  5:69671/1‑149 (MQ=255)
                                                                                   gcagcagcagcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATg                                                          >  4:137475/1‑149 (MQ=255)
                                                                                        agcagcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGa                                                     <  7:245028/149‑1 (MQ=255)
                                                                                           agcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAg                                                  >  8:270059/1‑149 (MQ=255)
                                                                                                 aTGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGcatcat                                            <  6:69671/149‑1 (MQ=255)
                                                                                                 aTGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGcatcat                                            <  3:207412/149‑1 (MQ=255)
                                                                                                     ggCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAgg                                        <  3:137475/149‑1 (MQ=255)
                                                                                                     ggCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAgg                                        <  1:234842/149‑1 (MQ=255)
                                                                                                         ttGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGa                                    <  7:242014/149‑1 (MQ=255)
                                                                                                             tCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCgg                                                                                                                 >  2:300512/1‑68 (MQ=255)
                                                                                                             tCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCgg                                                                                                                 <  1:300512/68‑1 (MQ=255)
                                                                                                             tCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACa                                                   >  6:346978/1‑130 (MQ=255)
                                                                                                             tCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACa                                                   <  5:346978/130‑1 (MQ=255)
                                                                                                              cAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATgc                                                         >  6:223163/1‑123 (MQ=255)
                                                                                                              cAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATgc                                                         <  5:223163/123‑1 (MQ=255)
                                                                                                                      tGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCga                       <  4:151267/149‑1 (MQ=255)
                                                                                                                           ggCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGc                            >  3:331483/1‑139 (MQ=255)
                                                                                                                           ggCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGc                            <  4:331483/139‑1 (MQ=255)
                                                                                                                           ggCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAgatga                  <  4:196299/149‑1 (MQ=255)
                                                                                                                            gCTGAAGATTGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCgaga                     >  6:312509/1‑145 (MQ=255)
                                                                                                                            gCTGAAGATTGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCgaga                     <  5:312509/145‑1 (MQ=255)
                                                                                                                                           ggCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGa                                                                              <  6:123136/73‑1 (MQ=255)
                                                                                                                                           ggCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGa                                                                              >  5:123136/1‑73 (MQ=255)
                                                                                                                                           ggCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGaa  >  4:336119/1‑149 (MQ=255)
                                                                                                                                            gCTGCATTGCCGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGCCGCAGGGTAGACCCGAGGAAGTACGCCGGCGCACTCATGCTCGACAGCATCATCAGGCCGATGCCCAGAACCGAGATGATCGCGGCAATccagac  >  3:33206/1‑143 (MQ=255)
                                                                                                                                            gCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACCCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGaa  >  1:110062/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
TACGTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGAA  >  NC_002947/6043520‑6043807

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