Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,077,237 A→C 9.1% intergenic (+32/+76) PP_4470 → / ← mgtE Arc domain‑containing transcriptional regulator/magnesium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,077,2370AC 9.1% 216.3 / 10.3 74intergenic (+32/+76)PP_4470/mgtEArc domain‑containing transcriptional regulator/magnesium transporter
Reads supporting (aligned to +/- strand):  ref base A (28/39);  new base C (3/4);  total (31/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGT  >  NC_002947/5077094‑5077378
                                                                                                                                                  |                                                                                                                                             
ccTGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTgg                                                                                                                                          >  8:130690/1‑149 (MQ=255)
ccTGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTg                                                                                                                                           <  4:64886/148‑1 (MQ=255)
   gCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCtt                                                                                                                                       >  1:62382/1‑149 (MQ=255)
    cATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCttt                                                                                                                                      <  6:168496/149‑1 (MQ=255)
     aTGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCtttt                                                                                                                                     >  3:26543/1‑149 (MQ=255)
     aTGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCtttt                                                                                                                                     <  3:207058/149‑1 (MQ=255)
        aGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCGCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGTGTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCtttttt                                                                                                                                   <  6:346998/148‑1 (MQ=255)
         gcgcGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑TTACGCTGGCTTTTTcgt                                                                                                                                 <  4:116861/149‑1 (MQ=255)
         gcgcGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGACAGCAGAACGCCTTGTTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGTCTGCAGCGAACACAAAAAAGCCAGC‑‑‑TTACGCTGGCTTTTTcgt                                                                                                                                 <  8:355590/149‑1 (MQ=255)
          cgcgAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTtgtg                                                                                                                                <  3:154277/149‑1 (MQ=255)
             gAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGc                                                                                                                             <  3:109034/149‑1 (MQ=255)
              aGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCg                                                                                                                            >  7:272701/1‑149 (MQ=255)
                       aGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTgg                                                                                                                   >  6:319881/1‑149 (MQ=255)
                       aGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTgg                                                                                                                   <  2:62382/149‑1 (MQ=255)
                        tcttttCCGCGAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGGCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTgg                                                                                                                   <  6:362900/145‑1 (MQ=255)
                            ttCCGCCAACTGTCACACCGCCAGCAGAACGCCTTGGACGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGCCTGCAGCGAACACGAAAAAGCCACC‑‑‑GTAAGCTGGATTTTTTGTGGGCGTCTGACTgg                                                                                                                   >  5:362900/1‑144 (MQ=255)
                              ccGCCAGCTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAg                                                                                                            >  4:222945/1‑149 (MQ=255)
                               cGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGGCGCAGACGCGTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGGTTTTTTGTGGGCGTCTGACTGGCAGCAAgg                                                                                                           <  1:62985/149‑1 (MQ=255)
                                 ccAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGc                                                                                                          >  6:265228/1‑148 (MQ=255)
                                          ccACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCATGGGGAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTcc                                                                                                   >  1:316030/1‑149 (MQ=255)
                                          ccACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCt                                                                                                 <  1:1786/148‑1 (MQ=255)
                                           cACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGc                                                                                               >  8:196215/1‑149 (MQ=255)
                                           cACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGc                                                                                               >  5:84794/1‑149 (MQ=255)
                                           cACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGc                                                                                               <  1:12972/149‑1 (MQ=255)
                                           cACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACAAAAAAGCCAGC‑‑‑TTACGCTGGCt                                                                                                                                        >  3:146476/1‑108 (MQ=255)
                                           cACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACAAAAAAGCCAGC‑‑‑TTACGCTGGCt                                                                                                                                        <  4:146476/108‑1 (MQ=255)
                                           cACCCCCAGCTCAACGCCTTGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCg                                                                                              <  8:221551/149‑1 (MQ=255)
                                            aCCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCg                                                                                              <  5:265228/149‑1 (MQ=255)
                                                 cagcagAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGAtt                                                                                         <  6:84794/149‑1 (MQ=255)
                                                    cagAACGCCTTGGTCTCCCTAATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACTAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCTGTATTGtt                                                                                      <  6:133488/149‑1 (MQ=255)
                                                       aaCGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTcc                                                                                    >  1:74777/1‑148 (MQ=255)
                                                         cGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCca                                                                                  >  7:179812/1‑148 (MQ=255)
                                                           ccTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAg                                                                               <  7:130690/149‑1 (MQ=255)
                                                            cTTGGTCGCCCTCATCGCTCACGATGTGGATATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAgg                                                                              >  1:33100/1‑149 (MQ=255)
                                                            cTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAgg                                                                              >  1:276662/1‑149 (MQ=255)
                                                               ggTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGggag                                                                           >  3:106554/1‑149 (MQ=255)
                                                                    cccTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTg                                                                      <  4:106554/149‑1 (MQ=255)
                                                                     ccTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTgg                                                                     <  8:272701/149‑1 (MQ=255)
                                                                         aTCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCg                                                                                              <  6:371429/120‑1 (MQ=255)
                                                                         aTCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCg                                                                                              >  5:371429/1‑120 (MQ=255)
                                                                                  gATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGaa                                                        <  4:166657/149‑1 (MQ=255)
                                                                                  gATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGaa                                                        <  4:91122/149‑1 (MQ=255)
                                                                                   aTGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAAc                                                       <  2:33100/149‑1 (MQ=255)
                                                                                    tgtgGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACa                                                      <  7:196215/149‑1 (MQ=255)
                                                                                        gAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGt                                                  <  8:179812/149‑1 (MQ=255)
                                                                                           aTGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTacaaca                                                            >  5:62739/1‑136 (MQ=255)
                                                                                           aTGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTacaaca                                                            <  6:62739/136‑1 (MQ=255)
                                                                                           aTGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGc                                               <  2:74777/149‑1 (MQ=255)
                                                                                           aTGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGc                                               <  3:222945/149‑1 (MQ=255)
                                                                                               ccGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGtt                                                                                      >  7:173483/1‑106 (MQ=255)
                                                                                               ccGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGtt                                                                                      <  8:173483/106‑1 (MQ=255)
                                                                                                  cAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACt                                                                                                                     >  4:206413/1‑72 (MQ=255)
                                                                                                  cAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACt                                                                                                                     <  3:206413/72‑1 (MQ=255)
                                                                                                  cAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGccc                                         >  4:360628/1‑148 (MQ=255)
                                                                                                   aGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGccc                                         >  6:45575/1‑147 (MQ=255)
                                                                                                   aGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCaa                                       <  3:366658/149‑1 (MQ=255)
                                                                                                     aCGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAg                                      >  4:204949/1‑148 (MQ=255)
                                                                                                       gCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTacaa                                                              <  5:151059/122‑1 (MQ=255)
                                                                                                       gCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTacaa                                                              >  6:151059/1‑122 (MQ=255)
                                                                                                       gCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGaa                                   <  3:360628/149‑1 (MQ=255)
                                                                                                         ctcctGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGa                                 <  1:231088/149‑1 (MQ=255)
                                                                                                         ctcctGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGa                                 <  6:194988/149‑1 (MQ=255)
                                                                                                         ctcctGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGa                                 <  5:319881/149‑1 (MQ=255)
                                                                                                            ctGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATg                               <  2:281794/148‑1 (MQ=255)
                                                                                                                  gCTGCAGCGATCACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAac                        <  3:288573/149‑1 (MQ=255)
                                                                                                                          gAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTg                >  1:299643/1‑149 (MQ=255)
                                                                                                                            ccacaAAAAAGCCAGC‑‑‑TTACGCTGGCTTTTTCGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTc              >  8:345259/2‑149 (MQ=255)
                                                                                                                                gAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAgg                                                          <  1:44335/101‑1 (MQ=255)
                                                                                                                                gAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAgg                                                          >  2:44335/1‑101 (MQ=255)
                                                                                                                                 aaaaaGCCAGC‑‑‑TTACGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTca         >  7:182220/1‑149 (MQ=255)
                                                                                                                                 aaaaaGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTca         <  2:276662/149‑1 (MQ=255)
                                                                                                                                  aaaaGCCAGC‑‑‑TTACGCTGGCTTTTTCGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTcac        <  7:345259/149‑1 (MQ=255)
                                                                                                                                      gCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCATGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACAACCGCTGTCGCTCAcag      >  8:203618/1‑147 (MQ=255)
                                                                                                                                        cAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGt  >  3:61132/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CCTGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGC‑‑‑GTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGT  >  NC_002947/5077094‑5077378

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: