Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,248,705 A→G 6.4% V441A (GTG→GCG)  PP_3724 ← acyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,248,7050AG 6.4% 176.1 / 2.1 62V441A (GTG→GCG) PP_3724acyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base A (33/25);  new base G (2/2);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

TTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTG‑GCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCC  >  NC_002947/4248558‑4248851
                                                                                                                                                   |                                                                                                                                                   
ttGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCca                                                                                                                                                     >  5:138830/1‑148 (MQ=255)
    aTTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACAccc                                                                                                                                                >  5:323701/1‑149 (MQ=255)
        ccATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCa                                                                                                                                            <  7:193214/149‑1 (MQ=255)
         cATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAg                                                                                                                                           >  6:38292/1‑149 (MQ=255)
              gTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCcaca                                                                                                                                                   >  7:150429/1‑136 (MQ=255)
              gTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCcaca                                                                                                                                                   <  8:150429/136‑1 (MQ=255)
                ccAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGggcggc                                                                                                                                     >  7:118557/1‑148 (MQ=255)
                   aCAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATg                                                                                                                                 >  8:275907/1‑149 (MQ=255)
                         ggCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAggg                                                                                                                                         >  1:324704/1‑135 (MQ=255)
                         ggCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAggg                                                                                                                                         <  2:324704/135‑1 (MQ=255)
                               cGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATg                                                                                                                      >  8:226312/1‑148 (MQ=255)
                                gCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCa                                                                                                                    <  5:375463/149‑1 (MQ=255)
                                 cccGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAg                                                                                                                   <  6:291988/149‑1 (MQ=255)
                                 cccGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAg                                                                                                                   >  4:360514/1‑149 (MQ=255)
                                    gTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAg                                                                                                                                      >  1:116042/1‑127 (MQ=255)
                                    gTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAg                                                                                                                                      <  2:116042/127‑1 (MQ=255)
                                        gggCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCa                                                                                                            >  1:95077/1‑149 (MQ=255)
                                            cTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAgg                                                                                                                                          >  6:336341/1‑115 (MQ=255)
                                             tGCGGGTACCGCTGTGCCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCAcgcg                                                                                                        >  4:299628/1‑148 (MQ=255)
                                              gCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCa                                                                                                      >  1:407213/1‑149 (MQ=255)
                                                gggCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACt                                                                                                    <  8:198842/149‑1 (MQ=255)
                                                  gCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTgg                                                                                                  >  7:207690/1‑149 (MQ=255)
                                                   cacaGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTggc                                                                                                 <  6:407664/149‑1 (MQ=255)
                                                   cacaGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTggc                                                                                                 <  2:407213/149‑1 (MQ=255)
                                                    acaGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTggca                                                                                                <  6:323701/149‑1 (MQ=255)
                                                     caGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTggcag                                                                                               <  7:275907/149‑1 (MQ=255)
                                                       gCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCTCTCACGCGCACTggcaggc                                                                                             <  7:226312/149‑1 (MQ=255)
                                                              ccGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGtt                                                                                      >  8:22425/1‑149 (MQ=255)
                                                                gCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCGCACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTcc                                                                                    >  5:9528/1‑149 (MQ=255)
                                                                gCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCGCACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTcc                                                                                    >  5:3615/1‑149 (MQ=255)
                                                                gCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTcc                                                                                    <  2:115985/149‑1 (MQ=255)
                                                                  acaGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCgg                                                                                  <  5:379564/149‑1 (MQ=255)
                                                                        tCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGc                                                                            >  5:242795/1‑149 (MQ=255)
                                                                           cccTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGa                                                                         <  8:207690/149‑1 (MQ=255)
                                                                              tGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGt                                                                      >  7:204266/1‑149 (MQ=255)
                                                                                gTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGcccc                                                                                                                             <  7:390042/92‑1 (MQ=255)
                                                                                gTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGcccc                                                                                                                             >  8:390042/1‑92 (MQ=255)
                                                                                gTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGt                                                                      <  4:305128/147‑1 (MQ=255)
                                                                                 tCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCGGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGAc                                                                   <  8:118557/149‑1 (MQ=255)
                                                                                   cATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACg                                                                  >  2:339365/1‑148 (MQ=255)
                                                                                    aTCAACCGCTGGCGGTTCTTTTGGTAAGCGGCGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACg                                                                  >  2:235795/1‑147 (MQ=255)
                                                                                                     cTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTG‑Gc                                              >  2:166989/1‑149 (MQ=255)
                                                                                                       tttGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGaaa                                                                              >  8:266269/1‑116 (MQ=255)
                                                                                                       tttGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGaaa                                                                              <  7:266269/116‑1 (MQ=255)
                                                                                                        ttGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTG‑GCGt                                            >  4:349720/1‑148 (MQ=255)
                                                                                                        ttGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGACCTTGCCGTG‑GCGt                                            <  8:204266/148‑1 (MQ=255)
                                                                                                             aaGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTG‑GCGTTCAACg                                      >  3:191385/1‑149 (MQ=255)
                                                                                                              aGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTG‑GCGTTCAACgg                                     >  4:324800/1‑149 (MQ=255)
                                                                                                               gCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTG‑GCGTTCAACGGc                                    <  2:95077/149‑1 (MQ=255)
                                                                                                                cgacgaTGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCCCGCGCACTGGCAGGCAACAGTTCCGGGAACGCCGAGGTGACGAAGCCGATCTTGCCGTG‑GCGTTCAACGGCa                                   >  7:360198/1‑149 (MQ=255)
                                                                                                                  acgaTGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCtt                                                       <  3:133405/128‑1 (MQ=255)
                                                                                                                  acgaTGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCtt                                                       >  4:133405/1‑128 (MQ=255)
                                                                                                                     aTGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTG‑GCGTTCAACGGCATCGCa                              <  5:38292/149‑1 (MQ=255)
                                                                                                                         gcCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTG‑GCGTTCAACGGCATCGCAGTa                           >  8:72073/1‑148 (MQ=255)
                                                                                                                         gcCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTG‑GCGTTCAACGGCATCGCAGTa                           >  5:109693/1‑148 (MQ=255)
                                                                                                                         gcCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTG‑GCGTTCAACGGCATCGCAGTa                           >  1:209001/1‑148 (MQ=255)
                                                                                                                           ccccAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTG‑GCGTTCAACGGCATCGCAGTACAg                        >  7:127357/1‑149 (MQ=255)
                                                                                                                            cccAGGCCGCCAGGCCGTCTTCCGCACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTG‑GCGTTCAACGGCATCGCAGTACAGc                       <  6:9528/149‑1 (MQ=255)
                                                                                                                            cccAGGCCGCCAGGCCGTCTTCCGCACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTG‑GCGTTCAACGGCATCGCAGTACAGc                       <  6:3615/149‑1 (MQ=255)
                                                                                                                             ccaggccgccaggccgTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCaa                                                                                           >  7:262004/1‑81 (MQ=255)
                                                                                                                             ccaggccgccaggccgTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCaa                                                                                           <  8:262004/81‑1 (MQ=255)
                                                                                                                                ggccgccaggccgTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTG‑GCGTTCAACGGCATCGCAGTACAGCCCca                   <  2:209001/149‑1 (MQ=255)
                                                                                                                                   cgccaggccgTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTG‑GCGTTCAACGGCATCGCAGTACAGCCCcagca                <  1:235795/149‑1 (MQ=255)
                                                                                                                                            gTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTG‑GCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCg        >  2:10610/1‑148 (MQ=255)
                                                                                                                                                 ccACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTG‑GCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTc   >  2:90699/1‑148 (MQ=255)
                                                                                                                                                  cacaCCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGCGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTcc  <  1:90699/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTG‑GCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCC  >  NC_002947/4248558‑4248851

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