Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 503,832 A→T 7.9% intergenic (+35/‑116) PP_0414 → / → rpe polyamine ABC transporter permease/ribulose‑5‑phosphate 3‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947503,8320AT 7.9% 183.4 / 5.3 63intergenic (+35/‑116)PP_0414/rpepolyamine ABC transporter permease/ribulose‑5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (28/30);  new base T (2/3);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

GACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCATTCT  >  NC_002947/503686‑503979
                                                                                                                                                   |                                                                                                                                                    
gACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAg                                                                                                                                                    >  1:148618/1‑149 (MQ=255)
     ttGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑Ga                                                                                                                                              <  8:63724/149‑1 (MQ=255)
      tGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGACGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTGCCGGCCTCTTC‑CGGCTTTAGCCgcg                                                                                                                                               >  1:158966/1‑145 (MQ=255)
       gCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑Gaa                                                                                                                                             <  1:204688/148‑1 (MQ=255)
        ccGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTGCCGGCCTCTTC‑CGGCTTTAGCCGCGAaga                                                                                                                                           <  5:358492/149‑1 (MQ=39)
        ccGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑Gaa                                                                                                                                             >  6:160569/1‑147 (MQ=255)
          gcggcggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGACGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTGCCGGCCTCTTC‑CGGCTTTAGCCGCGAAGAgg                                                                                                                                         <  2:158966/149‑1 (MQ=38)
            ggcggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGcc                                                                                                                                       >  3:225387/1‑149 (MQ=255)
                gcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGca                                                                                                                                   <  3:67875/149‑1 (MQ=255)
                     cTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGt                                                                                                                              >  2:322391/1‑149 (MQ=255)
                        gCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTACCTGTGCCGGCCTCTTC‑CGGCTTTAGCCGCGAAGAGGCCGGTACAggcaa                                                                                                                            <  4:185837/148‑4 (MQ=37)
                             ttGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTgg                                                                                                                      <  1:322391/149‑1 (MQ=255)
                                cGTTCTCGGTGTTGTTGCTGCTGACCCTGGTGGGTTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTc                                                                                                                    <  4:225387/148‑1 (MQ=255)
                                       ggTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGGGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGacac                                                                                                               <  8:245995/146‑1 (MQ=255)
                                         tgttgtTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAtt                                                                                                          >  8:295564/1‑149 (MQ=255)
                                         tgttgtTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAtt                                                                                                          <  3:349278/149‑1 (MQ=255)
                                             gttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCAGGGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGa                                                                                                        >  1:303442/1‑148 (MQ=255)
                                              ttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATa                                                                                                      <  2:285819/148‑1 (MQ=255)
                                              ttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATAc                                                                                                     <  6:326779/149‑1 (MQ=255)
                                                   gctgACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCCACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGACGCCGGCACACGTTGCGTTGCTCGAGACCAATTGATACCAAt                                                                                                 >  6:303411/1‑148 (MQ=255)
                                                    ctgACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAAt                                                                                                 >  4:36966/1‑147 (MQ=255)
                                                    ctgACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCa                                                                                               <  6:122936/149‑1 (MQ=255)
                                                         ccTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAag                                                                                                                                            >  3:287910/1‑99 (MQ=255)
                                                         ccTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAag                                                                                                                                            <  4:287910/99‑1 (MQ=255)
                                                         ccTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGa                                                                                                        >  1:380457/1‑135 (MQ=255)
                                                          cTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGt                                                                                                                     >  7:299836/1‑121 (MQ=255)
                                                          cTGGAGTGGTGGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGt                                                                                                                     <  8:299836/121‑1 (MQ=255)
                                                             gAGGGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCAt                                                                                       <  8:324202/148‑1 (MQ=255)
                                                                   ttGCGAGGGCGTAGCGAGAAGCTGAGGACTCATCCGCCCGC‑TTGATGATGGCTTGCCTGTAAAGGCCTCATCGCGGCTAAAGCCG‑GAAGAGCCCGTGACAGGTAGCGTTCGTCGAGACCAAATCATACCAATCTGCATTCATAGCtata                                                                                >  6:22930/1‑146 (MQ=255)
                                                                        tGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGa                                                                           <  4:287366/149‑1 (MQ=255)
                                                                          ggCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCGTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGt                                                                         <  2:321401/149‑1 (MQ=255)
                                                                             gTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTa                                                                       >  6:282612/1‑148 (MQ=255)
                                                                                 cGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATg                                                                  >  4:289171/1‑149 (MQ=255)
                                                                                   agaAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATgcg                                                                >  4:226934/1‑149 (MQ=255)
                                                                                    gaAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCCTCGCTTTCCCTGAGTTACAATgcgc                                                               <  2:362382/149‑1 (MQ=255)
                                                                                        cTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCa                                                           >  4:162630/1‑149 (MQ=255)
                                                                                          gAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCAcc                                                         >  4:287059/1‑149 (MQ=255)
                                                                                                tCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAAt                                                   >  4:148222/1‑149 (MQ=255)
                                                                                                tCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAAt                                                   >  3:223001/1‑149 (MQ=255)
                                                                                                  aaCAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAAt                                                   <  2:303442/147‑1 (MQ=255)
                                                                                                    cAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGc                                               <  5:282612/149‑1 (MQ=255)
                                                                                                                atgatgGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTgc                                     <  3:226934/148‑1 (MQ=255)
                                                                                                                       cTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCa                             <  3:289171/149‑1 (MQ=255)
                                                                                                                        ttGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACGAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCGAACAGGTGCGCCATGcgg                            >  3:212607/1‑147 (MQ=255)
                                                                                                                           ccTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATgcgc                                                               <  5:34392/111‑1 (MQ=255)
                                                                                                                           ccTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATgcgc                                                               >  6:34392/1‑111 (MQ=255)
                                                                                                                            cTGTACCGGCCTCTTCGCGGCTAAAGCCT‑GATGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCTGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCt                        <  5:22930/149‑1 (MQ=255)
                                                                                                                            cTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAAc                                           >  4:193836/1‑130 (MQ=255)
                                                                                                                            cTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAAc                                           <  3:193836/130‑1 (MQ=255)
                                                                                                                               tACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACg                     <  3:36966/149‑1 (MQ=255)
                                                                                                                                 ccGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGc                    >  2:264536/1‑148 (MQ=255)
                                                                                                                                   ggCCTCTTC‑CGGCTTTAGCCGCGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTAt                 >  8:246894/1‑149 (MQ=39)
                                                                                                                                    gCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGAAACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTAtt                <  7:295564/149‑1 (MQ=255)
                                                                                                                                       tctTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGcc             >  5:361083/1‑149 (MQ=255)
                                                                                                                                          tCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGccccc          >  1:263485/1‑149 (MQ=255)
                                                                                                                                            gcgGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCt         <  5:303411/148‑1 (MQ=255)
                                                                                                                                             cgGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGc                                                                                             >  7:335716/1‑63 (MQ=255)
                                                                                                                                             cgGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGc                                                                                             <  8:335716/63‑1 (MQ=255)
                                                                                                                                              ggCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCa      >  2:35282/1‑149 (MQ=255)
                                                                                                                                              ggCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCa      <  3:287059/149‑1 (MQ=255)
                                                                                                                                                cTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTa                       >  7:44445/1‑130 (MQ=255)
                                                                                                                                                cTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTa                       <  8:44445/130‑1 (MQ=255)
                                                                                                                                                 tAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCt                        <  8:229032/128‑1 (MQ=255)
                                                                                                                                                  aaaGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCATtct  >  4:239770/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGC‑TTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCATTCT  >  NC_002947/503686‑503979

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