Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,361,441 G→A 6.8% V169M (GTG→ATG)  phoP → two‑component system DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,361,4410GA 6.8% 177.3 / 2.2 59V169M (GTG→ATG) phoPtwo‑component system DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (28/27);  new base A (2/2);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

GTTTCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGACACGGTGCGCGGCCTGGGC  >  NC_002947/1361298‑1361587
                                                                                                                                               |                                                                                                                                                  
gTTTCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGc                                                                                                                                                >  7:196533/1‑148 (MQ=255)
 tttCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCACCACCTGCTGGTGATGGcgc                                                                                                                                              <  1:196478/149‑3 (MQ=255)
   tCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCATCCGCAAGCAGGCCACCGTGGACGAAGAAGCCTTGGCCCTGTCCGCCTAAGAATACCGCATCCTCTATTTCCTCATGCGCCATCACCTGCAGGTGGTGGCCaa                                                                                                                                             <  6:225949/148‑1 (MQ=255)
     aCCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGa                                                                                                                                          >  6:202275/1‑149 (MQ=255)
     aCCCAGTCGACCATCCCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCACCACCTGCTGGTGATGGcgcatga                                                                                                                                          <  8:335607/149‑7 (MQ=255)
       ccAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATc                                                                                                                                        <  5:252182/149‑1 (MQ=255)
       ccAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATc                                                                                                                                        <  7:92452/149‑1 (MQ=255)
         aGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGc                                                                                                                                      >  2:110887/1‑149 (MQ=255)
         aGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGc                                                                                                                                      >  2:189671/1‑149 (MQ=255)
          gTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACTAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATTTTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGcc                                                                                                                                     <  1:87488/149‑1 (MQ=255)
             gACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGa                                                                                                                                  <  6:101513/149‑1 (MQ=255)
               ccATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATg                                                                                                                                >  3:257397/1‑149 (MQ=255)
               ccATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATg                                                                                                                                <  1:189671/149‑1 (MQ=255)
                 aTCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCTAGCAGGCCACCCGCGACGAACAAGCCTTGGCCCTGCCCGCCTACGAATACCGCAGCCTCGCATACCTCCTGCGCCCTCACCCGCAGGTGGTGGCCAAGGATCGCCGGATg                                                                                                                                >  2:125589/1‑147 (MQ=255)
                 aTCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGa                                                                                                                              <  6:186335/149‑1 (MQ=255)
                 aTCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGa                                                                                                                              <  1:110887/149‑1 (MQ=255)
                    gccgcTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCa                                                                                                                           <  7:348766/149‑1 (MQ=255)
                    gccgcTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCa                                                                                                                           <  7:343842/149‑1 (MQ=255)
                      cgcTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGt                                                                                                                         >  8:321223/1‑149 (MQ=255)
                       gcTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGtt                                                                                                                        <  8:303599/149‑1 (MQ=255)
                       gcTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGtt                                                                                                                        <  6:188133/149‑1 (MQ=255)
                       gcTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGtt                                                                                                                        <  5:202275/149‑1 (MQ=255)
                           gCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTa                                                                                                                     <  1:12223/148‑1 (MQ=255)
                             cccTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTAcc                                                                                                                   >  4:253944/1‑148 (MQ=255)
                              ccTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCa                                                                                                                 <  2:24158/149‑1 (MQ=255)
                                         gACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGaa                                                                                                        >  6:178101/1‑147 (MQ=255)
                                             tCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGCCCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGCAGCACTTGAACCCAGGAGACGAAGAgcgc                                                                                                  >  3:231957/1‑149 (MQ=255)
                                               aaCCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCATGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAaga                                                                                                      >  8:338945/1‑143 (MQ=255)
                                               aaCCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAaga                                                                                                      <  7:338945/143‑1 (MQ=255)
                                                   gCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGAccc                                                                                             >  2:251624/1‑148 (MQ=255)
                                                            ccACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATc                                                                                   <  2:186459/149‑1 (MQ=255)
                                                               cccTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGaa                                                                                <  8:196533/149‑1 (MQ=255)
                                                               cccTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGaa                                                                                >  8:355260/1‑149 (MQ=255)
                                                                    gacgaACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCt                                                                           >  3:171285/1‑149 (MQ=255)
                                                                             gCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGgccg                                                                  >  1:257902/1‑149 (MQ=255)
                                                                                  ggCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCt                                                               <  3:165640/147‑1 (MQ=255)
                                                                                     ccTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGcgcc                                                          >  7:193786/1‑149 (MQ=255)
                                                                                     ccTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGcgcc                                                          >  8:192463/1‑149 (MQ=255)
                                                                                       tGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGcgccg                                                         >  3:112722/1‑148 (MQ=255)
                                                                                         aCCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGcgccgc                                                        <  1:67838/147‑1 (MQ=255)
                                                                                         aCCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCaa                                                      <  8:302368/149‑1 (MQ=255)
                                                                                             ccTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCtt                                                  >  2:219441/1‑149 (MQ=255)
                                                                                               tACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTg                                                 >  3:276635/1‑148 (MQ=255)
                                                                                                  gAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTAGAAgg                                             >  5:317326/1‑149 (MQ=255)
                                                                                                   aaTACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGc                                            <  8:193786/149‑1 (MQ=255)
                                                                                                     tACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGc                                            >  7:243505/1‑147 (MQ=255)
                                                                                                       ccGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAgc                                        >  8:266687/1‑149 (MQ=255)
                                                                                                                       cTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGc                                            <  8:315618/129‑1 (MQ=255)
                                                                                                                       cTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGc                                            >  7:315618/1‑129 (MQ=255)
                                                                                                                         ccttggcCACCACCTGCTGGTGATGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATc                       >  3:332365/6‑148 (MQ=255)
                                                                                                                           ttggcCACCACCTGCTGGTGATGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGa                     >  3:249929/4‑148 (MQ=255)
                                                                                                                             cgcCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCt                                                                           <  6:321558/92‑1 (MQ=255)
                                                                                                                             cgcCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCt                                                                           >  5:321558/1‑92 (MQ=255)
                                                                                                                               ccATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGACACgg                <  4:112722/149‑1 (MQ=255)
                                                                                                                                 aTCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGACACGGt               >  5:109048/1‑148 (MQ=255)
                                                                                                                                    aCCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGACACGGTgcgc           <  4:171285/149‑1 (MQ=255)
                                                                                                                                        gcagGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGACACGGTGCGCGGcc       >  3:169523/1‑149 (MQ=255)
                                                                                                                                        gcagGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGACACGGTGCGCGGcc       >  1:175016/1‑149 (MQ=255)
                                                                                                                                          agGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGACACGGTGCGCGGCCTg     <  3:253944/149‑1 (MQ=255)
                                                                                                                                           ggtggtGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGACACGGTGCGCGGCCTgg    <  7:192463/149‑1 (MQ=255)
                                                                                                                                           ggtggtGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGACACGGTGCGCGGCCTgg    <  4:169523/149‑1 (MQ=255)
                                                                                                                                            gtggtgGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGACACGGTGCGCGGCCTggg   <  6:317326/149‑1 (MQ=255)
                                                                                                                                            gtggtgGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGACACGGTGCGCGGCCTggg   <  5:178101/149‑1 (MQ=255)
                                                                                                                                              ggtggCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGACACGGTGCGCGGCCTgcgc  >  4:241878/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GTTTCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGACACGGTGCGCGGCCTGGGC  >  NC_002947/1361298‑1361587

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