Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 503,831 A→T 13.4% intergenic (+34/‑117) 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,8310AT13.4% 175.6 / 16.6 73intergenic (+34/‑117)PP_0414/rpepolyamine ABC transporter permease/ribulose‑5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (25/38);  new base T (4/6);  total (29/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

CCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTC  >  NC_002947/503684‑503973
                                                                                                                                                   |                                                                                                                                                
ccGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTa                                                                                                                                                  >  8:108907/1‑148 (MQ=255)
 cGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTa                                                                                                                                                  >  1:199820/1‑147 (MQ=255)
  gACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAg                                                                                                                                               <  8:162694/149‑1 (MQ=255)
  gACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAg                                                                                                                                               >  5:76097/1‑149 (MQ=255)
      aTTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCg                                                                                                                                            >  2:29729/1‑148 (MQ=255)
       ttGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGa                                                                                                                                         <  3:353615/149‑1 (MQ=255)
         gCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGGTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGaa                                                                                                                                        <  8:80726/148‑1 (MQ=255)
         gCCGCGGCGGCGACCTTGCTGATGGCGGTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGATGAGGACTCAACAGCCGGCTTTATGATGGCTTGCGTGTACCGGCCTCTTCGCGGCTAAAGCCGGa                                                                                                                                         <  8:353123/147‑1 (MQ=255)
          ccGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTACCTGTGCCGGCCTCTTCCGGCTTTAGCCGCGAaga                                                                                                                                      <  7:256344/149‑1 (MQ=38)
           cgcgGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTATTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAagag                                                                                                                                     >  2:30787/1‑149 (MQ=255)
            gcggcggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTGGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTGCCGGCCTCTTCCGGCTTTAGCCGCGAAGAgg                                                                                                                                    <  1:280113/149‑1 (MQ=38)
                 ggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGc                                                                                                                              <  8:221065/149‑1 (MQ=255)
                      ccTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTGCCGGCCTCTTCCGGCTTTAGCCGCGAAGAGGCCGGTACAgg                                                                                                                         <  3:291222/149‑1 (MQ=37)
                            tGATTGCTTTCTCGGTGTTATTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGt                                                                                                                   <  4:318925/149‑1 (MQ=255)
                              aTTGCGTTCTCGGTGTTTTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTg                                                                                                                 <  1:30787/149‑1 (MQ=255)
                              aTTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTg                                                                                                                 <  6:76097/149‑1 (MQ=255)
                                 gCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGt                                                                                                               <  1:29729/148‑1 (MQ=255)
                                 gCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGt                                                                                                               <  2:29659/148‑1 (MQ=255)
                                      ctcGGGGTTGTTGCTGCGGACCCTGGAGGGGTGGCGGGGGCGTAGCGATAAGCGGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAACGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGaca                                                                                                          <  3:364761/148‑1 (MQ=255)
                                       tcGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGacac                                                                                                         >  2:345661/1‑148 (MQ=255)
                                         ggTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCaa                                                                                                      >  4:153144/1‑149 (MQ=255)
                                          gtgtTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCaaa                                                                                                     <  2:50764/149‑2 (MQ=255)
                                           tgttgtTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAtt                                                                                                    >  4:122093/1‑149 (MQ=255)
                                                 tgctgcTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATAc                                                                                               >  1:37195/1‑148 (MQ=255)
                                                   ctgctgACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCACGGCACAGGTTGCGTTGGTCGACACCAATTGATACCa                                                                                             <  1:69479/149‑1 (MQ=255)
                                                     gctgACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGATGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGACTCTTCGCGGCTAAAGCCGGAAGAGGCCCGCACAGGTTGCTTTGGTCGACACCAATTGATACCAAt                                                                                           >  5:289598/1‑148 (MQ=255)
                                                      ctgACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAAt                                                                                           <  3:153144/147‑1 (MQ=255)
                                                        gACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCACATCGCCGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATAAGc                                                                                       >  5:262770/1‑149 (MQ=255)
                                                            cTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTc                                                                                   <  1:345661/149‑1 (MQ=255)
                                                             tGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAAt                                                                                           >  3:189963/1‑140 (MQ=255)
                                                             tGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAAt                                                                                           <  4:189963/140‑1 (MQ=255)
                                                               gAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCAt                                                                                 <  3:122093/148‑1 (MQ=255)
                                                               gAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCAt                                                                                 <  2:199820/148‑1 (MQ=255)
                                                               gAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATc                                                                                <  6:92239/149‑1 (MQ=255)
                                                                         gTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATGAGCCGTCATCGCTTTCCCt                                                                       >  5:355656/1‑148 (MQ=255)
                                                                         gTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCt                                                                       >  3:214870/1‑148 (MQ=255)
                                                                          tGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGa                                                                     <  5:123964/149‑1 (MQ=255)
                                                                            ggCGTAGCGAGTCGCTGAGGAAGCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCGTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCAGAGGTTGCGTTGGTCGACACCAAATGATACCAATCAGCCTTCATCGCTTTCCCTGAGt                                                                   <  6:196574/149‑1 (MQ=255)
                                                                                        aGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATgcgcgc                                                       <  2:37195/149‑1 (MQ=255)
                                                                                            gAGGACTCAACAGCCGGCGTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCTCCAcc                                                   <  8:119664/149‑1 (MQ=255)
                                                                                                 cTCAACAGCCGGCTTGAAGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGaa                                              <  3:319066/149‑1 (MQ=255)
                                                                                                  tCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCa                                                                                  >  6:71810/1‑112 (MQ=255)
                                                                                                  tCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGGCGACACCAATTGATACCAATCAGCCTTCa                                                                                  <  5:71810/112‑1 (MQ=255)
                                                                                                    aaCAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCCCGCCACCGTGAAttt                                           >  2:1902/1‑149 (MQ=255)
                                                                                                    aaCAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCATCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCCTTGGTCGACGCCAATTGAAACCAATCAGCCTTCATCCATTTCCCTGAGGTACAATGCGCGCCACCGTGAAttt                                           >  5:7638/1‑149 (MQ=255)
                                                                                                     aCAGCCGGCGTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTg                                          <  6:262770/149‑1 (MQ=255)
                                                                                                      cAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGTCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGc                                         <  6:309918/149‑1 (MQ=255)
                                                                                                          cGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTACCTGAGTTACAATGCGCGCCACCGTGAATTTGCCaa                                      >  7:145195/1‑148 (MQ=255)
                                                                                                            gCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAg                                   >  2:122166/1‑149 (MQ=255)
                                                                                                            gCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAg                                   <  7:108907/149‑1 (MQ=255)
                                                                                                                 atgatgGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTgc                               <  1:1902/148‑1 (MQ=255)
                                                                                                                   gttgGCTTGCCTGTACCGGCCTCTTCGCGTCTAAATCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTgcgc                             <  8:145195/146‑1 (MQ=255)
                                                                                                                   gatgGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTgcgc                             >  7:41109/1‑148 (MQ=255)
                                                                                                                      ggCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCa                                                     >  3:218174/1‑121 (MQ=255)
                                                                                                                      ggCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCa                                                     <  4:218174/121‑1 (MQ=255)
                                                                                                                      ggCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATg                         >  5:121185/1‑149 (MQ=255)
                                                                                                                       gCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACGATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATg                         >  5:388413/1‑148 (MQ=255)
                                                                                                                       gCTTGCCTGTACCGGCCTCTTCCCGGCTAAAGCCGGAAGAGGCCGGCTCAGGTTGCGTTGGTCGACGCCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATg                         <  6:7638/148‑1 (MQ=255)
                                                                                                                        cTTGCCTGTACCGGCCTCTTCGAGGATAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTTGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCTCCACCGTGAATTTGCCAACAGGTGCGCCATGc                        <  6:289598/148‑1 (MQ=255)
                                                                                                                         ttGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAg                      >  1:60689/1‑149 (MQ=255)
                                                                                                                         ttGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAg                      >  7:318441/1‑149 (MQ=255)
                                                                                                                                tACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACg               <  6:121185/149‑1 (MQ=255)
                                                                                                                                   cGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGcc                    <  5:146814/141‑1 (MQ=255)
                                                                                                                                   cGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGcc                    >  6:146814/1‑141 (MQ=255)
                                                                                                                                   cGGCCTCTTCCGGCTTTAGCCGCGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTa            >  5:356738/1‑149 (MQ=39)
                                                                                                                                     gCCTCTTCCGGCTTTAGCCGCGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTAc                                                                <  4:288327/95‑1 (MQ=37)
                                                                                                                                     gCCTCTTCCGGCTTTAGCCGCGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTAc                                                                >  3:288327/1‑95 (MQ=37)
                                                                                                                                      ccTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTg         <  4:214870/149‑1 (MQ=255)
                                                                                                                                         tcttccgGCTTTAGCCGCGAAGAGGCCGGTACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGcc       >  8:35825/6‑149 (MQ=39)
                                                                                                                                          ttCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGa                                                                     <  4:283064/85‑1 (MQ=255)
                                                                                                                                          ttCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGa                                                                     >  3:283064/1‑85 (MQ=255)
                                                                                                                                          ttCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGcc       <  8:41109/147‑1 (MQ=255)
                                                                                                                                             ccgGCTTTAGCCGCGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTAtt          >  4:100845/2‑142 (MQ=255)
                                                                                                                                             ccgGCTTTAGCCGCGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTAtt          <  3:100845/141‑1 (MQ=255)
                                                                                                                                              cgGCTTTAGCCGCGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTc  <  6:356738/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTC  >  NC_002947/503684‑503973

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: