Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,948,683 T→A 7.1% intergenic (‑64/+18) PP_mr62 ← / ← trxA Pseudomon‑Rho regulatory region/thioredoxin I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,948,6830TA 7.1% 248.6 / 7.5 84intergenic (‑64/+18)PP_mr62/trxAPseudomon‑Rho regulatory region/thioredoxin I
Reads supporting (aligned to +/- strand):  ref base T (41/37);  new base A (4/2);  total (45/39)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

GATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGCCGGGG  >  NC_002947/5948541‑5948828
                                                                                                                                               |                                                                                                                                                 
gATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCgg                                                                                                                                              <  7:25857/148‑1 (MQ=255)
      gACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCttt                                                                                                                                       >  3:414655/1‑149 (MQ=255)
       aCAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCtttt                                                                                                                                      <  7:17210/149‑1 (MQ=255)
         aGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTtca                                                                                                                                    <  5:95772/149‑1 (MQ=255)
              cgcgTACTGATTCGAATTATGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGGCCCCGCAAATGCGGGGCttttcttttt                                                                                                                                 <  4:297216/148‑7 (MQ=255)
                    cTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGt                                                                                                                          <  4:414655/148‑1 (MQ=255)
                     tGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTgg                                                                                                                        >  3:82501/1‑149 (MQ=255)
                       aTTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGTCACTGGATAAGCAGTGCGAGGGCGAATGTAACTCTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCTCAGGTg                                                                                                                         <  4:285477/146‑1 (MQ=255)
                        ttCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTgg                                                                                                                        <  1:217330/146‑1 (MQ=255)
                        ttCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTgg                                                                                                                        >  2:217330/1‑146 (MQ=255)
                           gAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTc                                                                                                                   <  4:112605/148‑1 (MQ=255)
                             aTTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAg                                                                                                                <  5:133931/149‑1 (MQ=255)
                                aGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGaa                                                                                                             >  7:305878/1‑149 (MQ=255)
                                 gcgATGCTGCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAAt                                                                                                            <  8:193198/147‑1 (MQ=255)
                                      gCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGc                                                                                                       >  5:434422/1‑149 (MQ=255)
                                      gCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGc                                                                                                       >  1:170875/1‑149 (MQ=255)
                                          cGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAg                                                                                                    <  3:421968/148‑1 (MQ=255)
                                             cGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGc                                                                                                >  1:2224/1‑149 (MQ=255)
                                                cGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGAt                                                                                             >  2:381169/1‑149 (MQ=255)
                                                    aTGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCggg                                                                                                                                             >  1:372866/1‑97 (MQ=255)
                                                    aTGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCggg                                                                                                                                             <  2:372866/97‑1 (MQ=255)
                                                      gTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGa                                                                                        >  6:193474/1‑148 (MQ=255)
                                                       tAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGa                                                                                        >  6:309902/1‑147 (MQ=255)
                                                         gagGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAGAGAAAAGACCCGCATTTGCGGGGACTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGt                                                                                    >  7:379291/1‑149 (MQ=255)
                                                         gagGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTATCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGt                                                                                    >  1:97386/1‑149 (MQ=255)
                                                           gggCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTgc                                                                                  >  8:86080/1‑149 (MQ=255)
                                                                   gAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTtc                                                                                                                                     <  1:133711/90‑1 (MQ=255)
                                                                   gAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTtc                                                                                                                                     >  2:133711/1‑90 (MQ=255)
                                                                   gAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCt                                                                          >  1:383566/1‑149 (MQ=255)
                                                                     aaaGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTtg                                                                        >  3:301029/1‑149 (MQ=255)
                                                                       aGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTtggt                                                                      <  1:381169/149‑1 (MQ=255)
                                                                        gCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTtggtg                                                                     >  7:210918/1‑149 (MQ=255)
                                                                         cAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTtggtgg                                                                    >  4:43507/1‑149 (MQ=255)
                                                                                  gCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGt                                                           >  2:11082/1‑149 (MQ=255)
                                                                                   ccGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGtt                                                          >  2:293606/1‑149 (MQ=255)
                                                                                      aaTGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGcc                                                       <  3:43507/149‑1 (MQ=255)
                                                                                          tAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTc                                                   >  4:138792/1‑149 (MQ=255)
                                                                                            acacTTGAGTTTTTCTGCGTCTAGTGATCTGAAATAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTTGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCt                                                  <  8:379291/148‑1 (MQ=255)
                                                                                              acTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTc                                                   <  6:358521/145‑1 (MQ=255)
                                                                                              acTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTc                                                   >  5:358521/1‑145 (MQ=255)
                                                                                              acTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGa                                               <  8:210918/149‑1 (MQ=255)
                                                                                               cTTGCTTTTTTCAGCGTCTAGTGCTATGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGaa                                              <  8:173520/149‑1 (MQ=255)
                                                                                                ttGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGaa                                              >  4:169339/1‑148 (MQ=255)
                                                                                                  gCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACa                                            <  8:305878/148‑1 (MQ=255)
                                                                                                  gCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACa                                            <  1:258748/148‑1 (MQ=255)
                                                                                                      ttttCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCAt                                        >  6:87504/1‑148 (MQ=255)
                                                                                                           tGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTtg                                                                        >  1:273404/1‑111 (MQ=255)
                                                                                                           tGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTtg                                                                        <  2:273404/111‑1 (MQ=255)
                                                                                                            gCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAATCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTc                                 >  1:127474/1‑149 (MQ=255)
                                                                                                                  aGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCgg                               <  7:121694/145‑1 (MQ=255)
                                                                                                                  aGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCgg                               <  7:127930/145‑1 (MQ=255)
                                                                                                                  aGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCgg                               >  8:121694/1‑145 (MQ=255)
                                                                                                                  aGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCgg                               >  8:127930/1‑145 (MQ=255)
                                                                                                                   gTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGAt                            <  6:434422/147‑1 (MQ=255)
                                                                                                                    tGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATAc                          >  3:402641/1‑148 (MQ=255)
                                                                                                                    tGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATAc                          >  3:396264/1‑148 (MQ=255)
                                                                                                                      ctctGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAATCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCAc                       <  2:127474/149‑1 (MQ=255)
                                                                                                                      ctctGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCAc                       <  2:383566/149‑1 (MQ=255)
                                                                                                                      ctctGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCAc                       <  3:138792/149‑1 (MQ=255)
                                                                                                                         tGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCgg                               <  1:129817/138‑1 (MQ=255)
                                                                                                                         tGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCgg                               >  2:129817/1‑138 (MQ=255)
                                                                                                                          gAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGca                    >  6:438543/1‑148 (MQ=255)
                                                                                                                           aaaaaGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACCGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGca                    <  5:309902/147‑1 (MQ=255)
                                                                                                                           aaaaaGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGcaca                  >  2:25127/1‑149 (MQ=255)
                                                                                                                           aaaaaGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGcaca                  >  5:398326/1‑149 (MQ=255)
                                                                                                                           aaaaaGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGcaca                  >  1:417952/1‑149 (MQ=255)
                                                                                                                             aaagaaaagCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTc                                                                                                                   <  2:219680/50‑1 (MQ=255)
                                                                                                                             aaagaaaagCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTc                                                                                                                   >  1:219680/1‑50 (MQ=255)
                                                                                                                             aaagaaaagCCCCGAATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCAAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGAATCAGAGTCGGGATACAACGCACAcc                >  6:112048/1‑149 (MQ=255)
                                                                                                                              aagaaaagCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCg               >  7:260092/1‑149 (MQ=255)
                                                                                                                               gaaaaggCCCCGCAAATGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCATCCTTGGTGGCTTCGACGTTGCCGTTCTTGaa                                              <  6:332894/115‑1 (MQ=255)
                                                                                                                               gaaaaggCCCCGCAAATGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCACCCTTGGTGGCTTCGACGTTGCCGTTCTTGaa                                              >  5:332894/3‑117 (MQ=255)
                                                                                                                               gaaaaggCCCCGCAAATGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGt              >  4:427009/3‑149 (MQ=255)
                                                                                                                               gaaaaggCCCCGCAAATGCGGGGCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTg             >  2:33248/3‑149 (MQ=255)
                                                                                                                                gaaaagCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTg             >  7:189563/1‑149 (MQ=255)
                                                                                                                                aaaaggCCCCGCAAATGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTg             >  5:195377/2‑149 (MQ=255)
                                                                                                                                  aaagCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCt           >  5:159547/1‑149 (MQ=255)
                                                                                                                                         ccGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGCCg     <  4:396264/149‑1 (MQ=255)
                                                                                                                                         ccGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGCCg     <  4:402641/149‑1 (MQ=255)
                                                                                                                                          cGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGCCgg    <  3:169339/149‑1 (MQ=255)
                                                                                                                                          cGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGCCgg    <  1:11082/149‑1 (MQ=255)
                                                                                                                                           gCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAAGCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGCCggg   <  4:423918/149‑1 (MQ=255)
                                                                                                                                            cATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGCCgggg  <  3:77958/149‑1 (MQ=255)
                                                                                                                                            cATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGCCgggg  <  2:97386/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGCCGGGG  >  NC_002947/5948541‑5948828

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