Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 579,606 A→T 5.9% T209S (ACG→TCG)  fdoI → formate dehydrogenase‑O subunit gamma

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947579,6060AT 5.9% 194.1 / 2.3 68T209S (ACG→TCG) fdoIformate dehydrogenase‑O subunit gamma
Reads supporting (aligned to +/- strand):  ref base A (30/32);  new base T (2/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

CCTTCGTGCTGATTCTCAGCATCATCGTGCACATCTACGCCGGCATCTGGATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGCAGATCGAAGCGTCGGCAGTCACGCCGCCATTTCTGCACCTGCCCGCCGCCAACCTGTTTGAACTGCGTGCCGCGCGCCTGGAGCAAC  >  NC_002947/579457‑579754
                                                                                                                                                     |                                                                                                                                                    
ccTTCGTGCTGATTCTCAGCATCATCGTGCACATCTACGCCGGCATCTGGATCAATGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCATGTGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGCCAAGa                                                                                                                                                      <  1:554345/150‑1 (MQ=255)
ccTTCGTGCTGATTCTCAGCATCATCGTGCACATCTACGCCGGCATCTGGATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGa                                                                                                                                                      <  2:138294/150‑1 (MQ=255)
     gTGCTGATTCTACGCATCATCGTGCACATCTACGCCGGCATCTGGATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCg                                                                                                                                                 <  1:264067/150‑1 (MQ=255)
      tGCTGATTCTCAGCATCATCGTTCACATCTACGCCGGCATCTGGATCAAGGGCTCAATCTGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCgg                                                                                                                                                <  1:13021/150‑1 (MQ=255)
      tGCTGATTCTCAGCATCATCGTGCACATCTACGCCGGCATCTGGATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCgg                                                                                                                                                <  1:946260/150‑1 (MQ=255)
       gCTGATTCTCAGCATCATCGTGCACATCTACGCCGGCATCTGGATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCg                                                                                                                                                 >  2:1192843/1‑148 (MQ=255)
       gCTGATTCTCAGCATCATCGTGCACATCTACGCCGGCATCTGGATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCg                                                                                                                                                 <  1:1192843/148‑1 (MQ=255)
        cTGATTCTCAGCATCATCGTGCACATCTACGCCGGCATCTGGATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGt                                                                                                                                              <  1:739856/150‑1 (MQ=255)
                 aGCATCATCGTGCACATCTACGCCGGCATCTGGATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGa                                                                                                                                     <  2:648367/150‑1 (MQ=255)
                      catcGTGCACATCTACGCCGGCATCTGGATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGaaa                                                                                                                                >  2:1081982/1‑150 (MQ=255)
                       atcGTGCACATCTACGCCGGCATCTGGATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGaaaa                                                                                                                               <  1:455888/150‑1 (MQ=255)
                        tcGTGCACATCTACGCCGGCATCTGGATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGAACTGGAAAGGCGAGATGGCCGCAGACAAGACACCGGTTGATCACGGACGAAaag                                                                                                                              >  2:977013/1‑150 (MQ=255)
                          gTGCACATCTACGCCGGCATCTGGATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAaagaa                                                                                                                            <  1:1216154/150‑1 (MQ=255)
                            gCACATCTACGCCGGCATCTGGATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAgg                                                                                                                          <  1:1081982/150‑1 (MQ=255)
                             cacaTCTACGCCGGCATCTGGATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGa                                                                                                                         >  1:193923/1‑150 (MQ=255)
                              acaTCTACGCCGGCATCTGGATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGAt                                                                                                                        <  2:135774/150‑1 (MQ=255)
                              acaTCTACGCCGGCATCTGGATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGAATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGAt                                                                                                                        >  2:442515/1‑150 (MQ=255)
                                  cTACGCCGGCATCTGGATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAgc                                                                                                                    >  2:483165/1‑150 (MQ=255)
                                      gCCGGCATCTGGATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGAATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCtt                                                                                                                <  1:442515/150‑1 (MQ=255)
                                            aTCTGGATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGt                                                                                                                                              >  2:699250/1‑114 (MQ=255)
                                            aTCTGGATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGt                                                                                                                                              <  1:699250/114‑1 (MQ=255)
                                              cTGGATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGcatca                                                                                                        <  2:193923/150‑1 (MQ=255)
                                               tGGATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGcatcat                                                                                                       >  2:382366/1‑150 (MQ=255)
                                                 gATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTCCCGCGAAGTGACCGGCGACCAGACGCCGGGTGATCACGGACGAAAAGACGGATGAGCGCTTGAGCATCATAc                                                                                                     <  2:1226115/150‑1 (MQ=255)
                                                 gATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAgcgc                                                                                                                  >  1:505136/1‑137 (MQ=255)
                                                 gATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAgcgc                                                                                                                  <  2:505136/137‑1 (MQ=255)
                                                 gATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATAc                                                                                                     >  2:953129/1‑150 (MQ=255)
                                                  aTCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACt                                                                                                    >  1:1098807/1‑150 (MQ=255)
                                                    cAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGaa                                                                                                                                 <  1:413369/119‑1 (MQ=255)
                                                    cAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGaa                                                                                                                                 >  2:413369/1‑119 (MQ=255)
                                                       gggCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAAc                                                                                               <  2:1007460/150‑1 (MQ=255)
                                                           tCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCgg                                                                                           >  2:467817/1‑150 (MQ=255)
                                                             gATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGc                                                                                         >  2:335367/1‑150 (MQ=255)
                                                               tCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGCAg                                                                                       >  1:49989/1‑150 (MQ=255)
                                                                cGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGCAGa                                                                                      <  1:953129/150‑1 (MQ=255)
                                                                     ccATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGACCCCGGGCAGATCGaa                                                                                 <  2:49989/150‑1 (MQ=255)
                                                                     ccATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGCAGATCGaa                                                                                 <  1:289376/150‑1 (MQ=255)
                                                                     ccATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGCAGATCGaa                                                                                 <  1:467817/150‑1 (MQ=255)
                                                                             cATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGCAGATCGAAGCGTCGGc                                                                         <  2:586152/150‑1 (MQ=255)
                                                                              aTGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGca                                                                                                           >  1:725416/1‑115 (MQ=255)
                                                                              aTGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGca                                                                                                           <  2:725416/115‑1 (MQ=255)
                                                                                    gggTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGCAGATCGAAGCGTCGGCAGTCAcg                                                                  >  1:1243273/1‑150 (MQ=255)
                                                                                      gTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGGGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGCAGATCGAAGCGTCGGCAGTCAcgcc                                                                >  1:1022467/1‑150 (MQ=255)
                                                                                       tCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGCAGATCGAAGCGTCGGCAGTCAcgccg                                                               >  1:919937/1‑150 (MQ=255)
                                                                                            cgcgCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGCAGATCGAAGCGTCGGCAGTCACGCCGCCAtt                                                          <  1:20302/150‑1 (MQ=255)
                                                                                             gcgcCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAgc                                                                                                                    <  1:765734/91‑1 (MQ=255)
                                                                                             gcgcCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAgc                                                                                                                    >  2:765734/1‑91 (MQ=255)
                                                                                             gcgcCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGCAGATCGAAGCGTCGGCAGTCACGCCGCCAttt                                                         >  1:1117018/1‑150 (MQ=255)
                                                                                                            aGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGCAGATCGAAGCGTCGGCAGTCACGCCGCCATTTCTGCACCTGCccgcc                                          >  2:1115447/1‑150 (MQ=255)
                                                                                                                  atGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCgg                                                                                           >  1:1053071/1‑95 (MQ=255)
                                                                                                                  atGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCgg                                                                                           <  2:1053071/95‑1 (MQ=255)
                                                                                                                  atGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGCAGATCGAAGCGTCGGCAGTCACGCCGCCATTTCTGCACCTGCCCGCCGCCAAc                                    <  1:382366/150‑1 (MQ=255)
                                                                                                                      ttcttttcgtccgtgatcacCCGGCGTCTTGTCGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGCAGATc                                                                                    >  2:408977‑M1/21‑98 (MQ=255)
                                                                                                                      ttcttttcgtccgtgatcacCCGGCGTCTTGTCGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGCAGATc                                                                                    <  1:408977‑M1/78‑1 (MQ=255)
                                                                                                                         gtgGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACATACACCCGGGCAGATCGAAGCGTCGGCAGTCACGCCGCCATTTCTGCACCTGCCCGCCGCCAACCTGTTTg                             <  1:730779/150‑1 (MQ=255)
                                                                                                                               ccGCGAAGTGACCGGCGACAAGGCGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGCAGATCGAAGCGTCGGCAGTCACGCCGCCATTTCTGCACCTGCCCGCCGCCAACCTGt                                <  1:1166498/141‑1 (MQ=255)
                                                                                                                               ccGCGAAGTGACCGGCGACAAGGCGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGCAGATCGAAGCGTCGGCAGTCACGCCGCCATTTCTGCACCTGCCCGCCGCCAACCTGt                                >  2:1166498/1‑141 (MQ=255)
                                                                                                                               ccGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGCAGATCGAAGCGTCGGCAGTCACGCCGCCATTTCTGCACCTGCCCGCCGCCAACCTGTTTGAACTGc                       >  2:299235/1‑150 (MQ=255)
                                                                                                                                cgcgAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGCAGATCGAAGCGTCGGCAGTCACGCCGCCATTTCTGCACCTGCCCGCCGCCAACCTGTTTGAACTGCg                      >  2:768711/1‑150 (MQ=255)
                                                                                                                                cgcgAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGAGCAGATCGAAGCGTCGGCAGTCACGCCGCCATTTCTGCACCTGCCCGCCGCAAACCTGTTTGAACTGCg                      >  2:135948/1‑150 (MQ=255)
                                                                                                                                   tgatcacCCGGCGTCTTGTCGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGaa                                                                                                <  2:202787‑M1/66‑1 (MQ=255)
                                                                                                                                   tgatcacCCGGCGTCTTGTCGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGaa                                                                                                >  1:202787‑M1/8‑73 (MQ=255)
                                                                                                                                           CGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGCAGATCGAAGCGTCGGCAGTCACGCCGCCATTTCTGCACCTGCCCGCCGCCAACCTGTTTGAACTGCGTGCCGCGCGcc           <  2:202799‑M1/150‑1 (MQ=255)
                                                                                                                                           CGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGCAGATCGAAGCGTCGGCAGTCACGCCGCCATTTCTGCACCTGCCCGCCGCCAACCTGTTTGAACTGCGTGCCGCGCGcc           <  2:1151470‑M1/150‑1 (MQ=255)
                                                                                                                                               GACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGCAGATCGAAGCGTCGGCAGTCACGCCGCCATTTCTGCACCTGCCCGCCGCCAACCTGTTTGAACTGCGTGCCGCGCGCCTGGa       <  2:758905‑M1/150‑1 (MQ=255)
                                                                                                                                                  aaGAGGCCGGGTGATCACGGACGAAAAGAACGATGAGCGCTTGAGCATCATACATGAACCCGGGCAGAGCGAAGAGTGGCCAGTCACGCCGCCCTTTCTGCACATGCCGGCCGCCAACCTGTATGAACTGCGTGTCGCCCGCCTgcagca    >  2:121425/1‑150 (MQ=255)
                                                                                                                                                  aaGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATAc                                                                                                     <  1:935752/53‑1 (MQ=255)
                                                                                                                                                  aaGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATAc                                                                                                     >  2:935752/1‑53 (MQ=255)
                                                                                                                                                    gACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGCAGATCGAAGCGTCGGCAGTCACGCCGACATTTCTCGCCCGCCCCGCCGCCAACCTGTGTGAACTGCGTGCCGCGCGCCTGGAGCAAc  >  2:997692/1‑150 (MQ=255)
                                                                                                                                                    gACGCCGGGTGATCACGGACGAAAACATGTAAGAGCGCTTGAGCATCATACTTGAACCCGGGCAGATCGAAGCGTCGGCAGTCACGCCGCCATTTCTGCACCTGCCCGCCGCCAACCTGTTTGAACTGCGTGCCGCGCGCCTGGAGCAAc  >  1:358717/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCTTCGTGCTGATTCTCAGCATCATCGTGCACATCTACGCCGGCATCTGGATCAAGGGCTCGATCGGCGCCATGCTGCATGGCTGGGTCAGCCGCGCCTGGGCCCGCAAGCATCATGAGTTGTGGTACCGCGAAGTGACCGGCGACAAGACGCCGGGTGATCACGGACGAAAAGAAGGATGAGCGCTTGAGCATCATACTTGAACCCGGGCAGATCGAAGCGTCGGCAGTCACGCCGCCATTTCTGCACCTGCCCGCCGCCAACCTGTTTGAACTGCGTGCCGCGCGCCTGGAGCAAC  >  NC_002947/579457‑579754

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