Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 802,773 G→A 13.9% intergenic (+31/‑55) 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,7730GA13.9% 174.6 / 17.4 72intergenic (+31/‑55)obg/proBGTPase Obg/glutamate 5‑kinase
Reads supporting (aligned to +/- strand):  ref base G (30/32);  new base A (6/4);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 7.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.79e-01

CCGTACCCTGCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGCCTCGACCGCGGTGCAA  >  NC_002947/802625‑802921
                                                                                                                                                    |                                                                                                                                                    
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                                                                                                                  <  7:199473/149‑1 (MQ=255)
                                    tGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTc                                                                                                                  >  5:37813/1‑149 (MQ=255)
                                    tGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTc                                                                                                                  <  3:170517/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                                                                                                                                          >  7:54795/1‑107 (MQ=39)
                                                     aTGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCCAGAGAGCGGCGttta                                                                                                                                          <  8:54795/108‑2 (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                                                        <  6:159139/149‑1 (MQ=255)
                                                                                              gAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTa                                                        <  8:21362/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)
                                                                                                                                                   tGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGCCTCGACCGCGGTGc    >  4:233875/1‑148 (MQ=255)
                                                                                                                                                    ggAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGCCTCGACCGCGGTGCaa  >  8:253091/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCGTACCCTGCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGCCTCGACCGCGGTGCAA  >  NC_002947/802625‑802921

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: ‑