Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 10,024 C→A G157G (GGG→GGT)  satP ← acetate/succinate:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091310,0240CA97.6% 122.2 / NA 41G157G (GGG→GGT) satPacetate/succinate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (24/16);  total (25/16)
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

GCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACGGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAGCGT  >  NC_000913/9783‑10263
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
gCCAGCGTACCGTACTGCATTCAGTGGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCa                                                                                                                                                                                                                                         <  1:435838/250‑1 (MQ=255)
 ccAGCGTACCGTACGGCATTCAGTGGCGGGAAGGGCCATACGTTGAAACGGGACCAGAAGTGGTTGCTGCATCCATAACGAAGAGTGAAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGTGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCCCCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGc                                                                                                                                                                                                                                       <  2:230994/251‑1 (MQ=255)
       taccgtacTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCa                                                                                                                                                                                                                                  >  1:15869/1‑250 (MQ=255)
               tGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGgatg                                                                                                                                                                                                                         >  1:479824/1‑251 (MQ=255)
                      aGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTgcgg                                                                                                                                                                                                                  <  1:127282/251‑1 (MQ=255)
                           atGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTAAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCATCCATCGCCAGATAGATTGCACTGGCACCAAAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGtt                                                                                                                                                                                                              >  2:451944/2‑250 (MQ=255)
                             ggAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTAACCAATCGGCAGAACGGTGCGACAAAACTGCTAGTTCAGTACTTCACCCATCGACAAATAGATTGCACTGGCCCCACAGATCAGACCAATCAAGCCGGCAAAGGGGAAGATTGCGGCGTAAcc                                                                                                                                                                                                           >  2:59016/1‑251 (MQ=255)
                                     cATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGCCCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGCCCATTCCAGCCTGCACCGTGGATGATTGCGGCGTTACCGGCAATGt                                                                                                                                                                                                   >  2:190326/1‑251 (MQ=255)
                                       tACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGtt                                                                                                                                                                                                  >  2:445676/1‑250 (MQ=255)
                                                               gTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGc                                                                                                                                                                         >  1:168033/1‑251 (MQ=255)
                                                                            gACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTa                                                                                                                                                             >  1:213061/1‑250 (MQ=255)
                                                                               cGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGACCGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGCTTGCACTGGCCCCACAGATCAGACCAATCCAGCCGGCACAGTGGATAATTGAGGCGTTACCGGCAATGATACCGATCGCCAACAGGGACAACAGCAcagtgagacgaa                                                                                                                                                            >  2:390383/1‑237 (MQ=255)
                                                                                 aaGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGaa                                                                                                                                                        >  1:296798/1‑250 (MQ=255)
                                                                                         aaGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTg                                                                                                                                               <  1:59818/251‑1 (MQ=255)
                                                                                                gACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAAcg                                                                                                                                        <  1:74885/251‑1 (MQ=255)
                                                                                                      aGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTgcg                                                                                                                                  <  1:445676/251‑1 (MQ=255)
                                                                                                      aGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTgcg                                                                                                                                  >  1:150242/1‑251 (MQ=255)
                                                                                                           aTAAAAAAATCCCCCCGGGCGGGGGGATCGCAAAACAATTAGTGTGATTCCCCAATCGGCAGAACGGTGAGCCCAAACGGCTCGTTTAGTTCTTCACCCATCGCCATATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGGGGATGATTGCTGCGTTACCGTCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAAGTGCAGAATGCGTGCGCCtgt                                                                                                                             <  2:294361/251‑3 (MQ=255)
                                                                                                            aaaaaaaaCCCCCCCGAGCGGGGGGATCTCAAAACAAGTAGTGGGATTCACCAAACGGCAGAACGGTGCGACCAAACTGCGCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTc                                                                                                                            <  2:213061/250‑1 (MQ=255)
                                                                                                               aaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAAAACGCGGGCCCCTTTAAGc                                                                                                                         >  2:431460/1‑251 (MQ=255)
                                                                                                                aaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAAc                                                                                                                                                     >  1:376422/1‑222 (MQ=255)
                                                                                                                aaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGGGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAAc                                                                                                                                                     <  2:376422/222‑1 (MQ=255)
                                                                                                                     ccccccGAGCGGGGGGATCACAAAACAATTAGTGGGATTCCCCAAGCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCg                                                                                                                   <  2:15869/251‑1 (MQ=255)
                                                                                                                       ccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGAGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGAGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCga                                                                                                                  >  1:358446/1‑250 (MQ=255)
                                                                                                                              cGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACAt                                                                                                           >  2:34216/1‑250 (MQ=255)
                                                                                                                               ggggggATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGCAGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGAAAACAGCACAGTCAGGCTAAAGAAAACAAATTGAAGAACGCGAGCGCCTTTCAGCGTGCAGAAGAACATa                                                                                                          >  2:344237/1‑250 (MQ=255)
                                                                                                                               gggggaATCTCAAAACAATTAGTGGGATTCAACAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGACAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACTCGTGAGCCGTTACGCGTGCagaagaccaca                                                                                                          >  2:414197/1‑245 (MQ=255)
                                                                                                                                gggggATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATaa                                                                                                         >  1:163278/1‑250 (MQ=255)
                                                                                                                                                             aCCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGAGCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGt                                                                           <  2:163278/251‑1 (MQ=255)
                                                                                                                                                                aaTCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGt                                                                           >  1:196703/1‑248 (MQ=255)
                                                                                                                                                                    ggCATAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACCGATCAGACCAATCCAGCCGGCAAAGTGGATGCTTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGAGCACACCGCACAGTCAGGCTAAAGAAAACGCATTGCAGAACGCGTGCGCCTTTCCGCGTGCCGAAGAACATAAACAGCCTAAATCCACCCCCCAGCCCCCGGTAGCCAcc                                                                    >  2:56208/1‑251 (MQ=255)
                                                                                                                                                                    ggCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACAcc                                                                    <  1:152323/251‑1 (MQ=255)
                                                                                                                                                                    ggCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACAcc                                                                    <  1:34216/251‑1 (MQ=255)
                                                                                                                                                                    ggCAGAACAGGGCGACCAAACTACTCGATCAGCACTTCACCCATCGCAAGATAGATTGCACTGGCACCACAGCTCAGACAAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCCGCAGGGCAAAAAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAAAGCGAGCGCCTTTCCGAGTGCGGTAGAACATAAACAGCGTAAATACGCCCCACATATCCAGGGAGACAcc                                                                    >  2:33572/1‑251 (MQ=255)
                                                                                                                                                                            ggTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACt                                                            <  1:124018/251‑1 (MQ=255)
                                                                                                                                                                                  aCCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGACTTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCAt                                                      >  2:108776/1‑251 (MQ=255)
                                                                                                                                                                                    cAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCAttt                                                    >  2:397184/1‑251 (MQ=255)
                                                                                                                                                                                     aaaCTGCTCGTTCAGTACTTCCCCCATCGCCAGATAGATTGCCCTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGGGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTg                                                   <  2:255332/251‑1 (MQ=255)
                                                                                                                                                                                                                      atagatTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATcagc                   >  1:279688/1‑250 (MQ=255)
                                                                                                                                                                                                                              cacTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGc           >  1:176041/3‑250 (MQ=255)
                                                                                                                                                                                                                                      ccacAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCTCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACTCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGt  <  1:321023/248‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
GCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACGGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAGCGT  >  NC_000913/9783‑10263

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑