Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 802,770 T→C 14.3% intergenic (+28/‑58) obg → / → proB GTPase Obg/glutamate 5‑kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947802,7700TC14.3% 171.8 / 17.0 70intergenic (+28/‑58)obg/proBGTPase Obg/glutamate 5‑kinase
Reads supporting (aligned to +/- strand):  ref base T (29/31);  new base C (6/4);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTACCCTGCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGCCTCGACCGCGG  >  NC_002947/802625‑802916
                                                                                                                                                 |                                                                                                                                                  
ccGTACCCTGCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTg                                                                                                                                                 >  3:62071/1‑149 (MQ=255)
  gTACCCTGCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGa                                                                                                                                               <  5:245083/149‑1 (MQ=255)
      ccTGCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTTGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAgcgg                                                                                                                                           >  7:21362/1‑149 (MQ=255)
      ccTGCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAgcgg                                                                                                                                           >  3:300760/1‑149 (MQ=255)
      ccTGCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAgcgg                                                                                                                                           >  6:326303/1‑149 (MQ=255)
          cgtcgtACCGGCGGCAAGAGCGGGCATGACATCGGCGTCGACGATGGTTGGGACGATGATTTCGAGCACGACGAAGACGGTCCGGTAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCGCTCTGGAgcggcg                                                                                                                                         <  4:44317/147‑1 (MQ=255)
           gtcgtACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGttt                                                                                                                                      <  5:326303/149‑1 (MQ=255)
            tcgtACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCCAGAGAGCGGCGttta                                                                                                                                     <  7:83217/149‑2 (MQ=255)
              gtACCTGAGTCAAGAGCGTGCATGACTTCGGCGACGACGTTGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGttttt                                                                                                                                    <  3:214743/148‑1 (MQ=255)
              gtACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGATGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGttttt                                                                                                                                    <  6:197470/148‑1 (MQ=255)
              gtACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGttttt                                                                                                                                    <  2:6673/148‑1 (MQ=255)
              gtACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGttttt                                                                                                                                    <  5:282078/148‑1 (MQ=255)
              gtACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGttttt                                                                                                                                    <  2:1135/148‑1 (MQ=255)
              gtACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGttttt                                                                                                                                    <  8:43057/148‑1 (MQ=255)
               tACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTg                                                                                                                                   >  6:184323/1‑148 (MQ=255)
                 ccGGCGTCAAGAGCGTGCATGACATTGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCCAGAGAGCGGCGtt                                                                                                                                       >  6:197183/1‑142 (MQ=255)
                 ccGGCGTCAAGAGCGTGCATGACATTGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCCAGAGAGCGGCGtt                                                                                                                                       <  5:197183/142‑1 (MQ=255)
                 ccGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTa                                                                                                                                 >  8:87175/1‑148 (MQ=255)
                  cGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATc                                                                                                                               >  6:177086/1‑149 (MQ=255)
                         aaGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTa                                                                                                                        >  8:267171/1‑149 (MQ=255)
                          agagCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTa                                                                                                                        >  8:199473/1‑148 (MQ=255)
                              cGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTAc                                                                                                                   <  5:92404/149‑1 (MQ=255)
                                   aTGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGAtt                                                                                                              <  8:16818/149‑1 (MQ=255)
                                    tGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTc                                                                                                             <  3:170517/149‑1 (MQ=255)
                                    tGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTc                                                                                                             >  5:37813/1‑149 (MQ=255)
                                    tGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTc                                                                                                             <  7:199473/149‑1 (MQ=255)
                                             cgacgacgaTGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGt                                                                                                    <  4:300760/149‑1 (MQ=255)
                                                cgacgaTGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTgg                                                                                                 <  4:62071/149‑1 (MQ=255)
                                                cgacgaTGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTgg                                                                                                 <  5:177086/149‑1 (MQ=255)
                                                    gaTGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGgaag                                                                                              >  5:75831/1‑148 (MQ=255)
                                                    gaTGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACAGAGGTTGgaag                                                                                              >  6:306748/1‑148 (MQ=255)
                                                    gaTGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTAGACATAGGTTGgaag                                                                                              >  5:310432/1‑148 (MQ=255)
                                                     aTGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCCAGAGAGCGGCGttta                                                                                                                                     <  8:54795/108‑2 (MQ=39)
                                                     aTGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCCAGAGAGCGGCGttta                                                                                                                                     >  7:54795/1‑107 (MQ=39)
                                                        gTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGAt                                                                                         <  5:184323/149‑1 (MQ=255)
                                                         ttGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATg                                                                                        >  8:226242/1‑149 (MQ=255)
                                                             gACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGaa                                                                                    <  7:267171/149‑1 (MQ=255)
                                                                  tgatTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAg                                                                               <  5:212698/149‑1 (MQ=255)
                                                                  tgatTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAg                                                                               <  7:59528/149‑1 (MQ=255)
                                                                        cGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGAc                                                                          <  3:258262/148‑1 (MQ=255)
                                                                         gAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACgg                                                                        >  5:159139/1‑149 (MQ=255)
                                                                               gacgaAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGcc                                                                   >  8:86281/1‑148 (MQ=255)
                                                                                      aCGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTg                                                           >  2:312557/1‑149 (MQ=255)
                                                                                        ggTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTggg                                                         >  6:295404/1‑149 (MQ=255)
                                                                                              gAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGtt                                                    >  3:300576/1‑148 (MQ=255)
                                                                                              gAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTa                                                   <  8:21362/149‑1 (MQ=255)
                                                                                              gAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTa                                                   <  6:159139/149‑1 (MQ=255)
                                                                                                 aTCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGa                                                <  7:86281/149‑1 (MQ=255)
                                                                                                     tttACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTgg                                            <  7:87175/149‑1 (MQ=255)
                                                                                                         cGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAg                                         <  1:312557/148‑1 (MQ=255)
                                                                                                            gcgcGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGc                                      >  2:153157/1‑148 (MQ=255)
                                                                                                             cgcgACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGct                                     <  4:300576/148‑1 (MQ=255)
                                                                                                                   tGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGa                              >  5:284375/1‑149 (MQ=255)
                                                                                                                     aTCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCt                                                            <  1:112767/117‑1 (MQ=255)
                                                                                                                     aTCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCt                                                            >  2:112767/1‑117 (MQ=255)
                                                                                                                            gCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGAt                      <  7:226242/148‑1 (MQ=255)
                                                                                                                             cAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGgcgc                                                                    >  7:160016/1‑101 (MQ=255)
                                                                                                                             cAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGgcgc                                                                    <  8:160016/101‑1 (MQ=255)
                                                                                                                             cAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCTGTGCTCTGCTGACCGCCGATgg                    >  2:76457/1‑149 (MQ=255)
                                                                                                                              aGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATgg                    >  6:172173/1‑148 (MQ=255)
                                                                                                                               gTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGc                   >  6:253504/1‑148 (MQ=255)
                                                                                                                                tACACTAAACGCCGCTCTGTGGAGCGGCGGTTTTGTATGCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGCAGGAAGGTGACGGGCGCCAAGCGCGGGGGCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATTGc                   >  2:339485/1‑147 (MQ=255)
                                                                                                                                  cacTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAg                <  2:262002/148‑1 (MQ=255)
                                                                                                                                    aaaaaCGCCGCTCCAGAGAGCGGCGTTTTTGTATCCTCGGTATGTACAGATTCGACATAGGTTGGAAGATGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGc             >  3:5521/3‑149 (MQ=255)
                                                                                                                                    aaaaaCGCCGCTCCAGAGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAg                <  3:161759/144‑1 (MQ=255)
                                                                                                                                    aaaaaCGCCGCTCCAGAGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAg                >  4:161759/3‑146 (MQ=255)
                                                                                                                                     tAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGcc            >  3:63600/1‑149 (MQ=255)
                                                                                                                                      aaaCGCCGCTCCAGAGAGCGGCGTTTTTGTATCCACGGTTTCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGCCt           >  8:10152/1‑149 (MQ=255)
                                                                                                                                       aaCGCCGCTCTCTGGAGCGGCGTTTTTCTATCCACGGTATCTACAGAATCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGCCTc          >  2:224107/1‑149 (MQ=255)
                                                                                                                                       aaCGCCGCTCCAGAGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGCCTc          >  4:255387/1‑149 (MQ=255)
                                                                                                                                              ctctctGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGCCTCGACcgcg   >  1:192432/1‑149 (MQ=255)
                                                                                                                                                 tctGGAGCGGAGTTTTTGTATCCACGGTATCTACAGATTCGACATAGTTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGCCTCGACCGCgg  <  6:310432/147‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CCGTACCCTGCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGCCTCGACCGCGG  >  NC_002947/802625‑802916

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: