Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 591,082 T→C 8.8% V110A (GTG→GCG)  PP_0501 → NAD‑dependent epimerase/dehydratase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947591,0820TC 8.8% 218.8 / 8.4 79V110A (GTG→GCG) PP_0501NAD‑dependent epimerase/dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base T (35/37);  new base C (3/4);  total (38/41)
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

CAGGGCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGAC  >  NC_002947/590935‑591227
                                                                                                                                                   |                                                                                                                                                 
cAGGGCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTg                                                                                                                                                  <  3:166382/149‑1 (MQ=255)
 aGGGCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTg                                                                                                                                                  >  4:243147/1‑148 (MQ=255)
   ggCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGGTCAGGCATAGGTCGAAGACCCGGGGCGGACCCCCCAGAGTTACGTGATCGGTTCCCGCAACGGTTGCGAAGCCCTGCACGCGCAt                                                                                                                                               >  6:190155/1‑149 (MQ=255)
        gTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCACGCGCATGGct                                                                                                                                           <  5:190155/148‑2 (MQ=255)
         tCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTg                                                                                                                                         <  3:37524/149‑1 (MQ=255)
         tCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCAGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCACGCGCATGGcttc                                                                                                                                         <  6:29000/149‑4 (MQ=255)
           aCTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTgcg                                                                                                                                       >  8:5042/1‑149 (MQ=255)
             tCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGcc                                                                                                                                     >  7:130747/1‑149 (MQ=255)
             tCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGcc                                                                                                                                     >  4:388779/1‑149 (MQ=255)
               aGGCTGCGGCTGGCTGCCGTGCTGTGGTGGATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCgg                                                                                                                                   <  6:47306/149‑1 (MQ=255)
               aGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCCTCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCgg                                                                                                                                   <  6:53318/149‑1 (MQ=255)
                ggCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCggg                                                                                                                                  <  2:10078/149‑1 (MQ=255)
                ggCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCggg                                                                                                                                  <  8:130747/149‑1 (MQ=255)
                        ctggctgCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTc                                                                                                                          <  5:96617/149‑1 (MQ=255)
                               ccGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGc                                                                                                                                            >  4:346264/1‑124 (MQ=255)
                               ccGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGc                                                                                                                                            <  3:346264/124‑1 (MQ=255)
                               ccGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTc                                                                                                                     <  6:162135/147‑1 (MQ=255)
                                    cTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCcagcgca                                                                                                              <  3:243147/149‑1 (MQ=255)
                                     tgtgGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCcagcgcag                                                                                                             <  2:7198/149‑1 (MQ=255)
                                       tggtgCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCcagcgcagc                                                                                                            <  7:5042/148‑1 (MQ=255)
                                            cATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGgtgt                                                                                                       <  1:386395/148‑1 (MQ=255)
                                             aTCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCAt                                                                                                                                               <  4:384247/107‑1 (MQ=255)
                                             aTCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCAt                                                                                                                                               >  3:384247/1‑107 (MQ=255)
                                             aTCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCAATCGCCTGCATGGCGTGCGCCgg                                                                                                                                   >  8:296362/1‑119 (MQ=255)
                                                  gcggcgGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGcaa                                                                                                >  6:278225/1‑149 (MQ=255)
                                                       ggTAGCGTCTGTTCAGGCATCGATCGAAGACCCGGTGCGTACCCACCAGAGTAACTTAATCGGTACCCTCAACCTTTGCCAAGCCATGCGCGTGCATGGCGTGCGCCCGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACCGCAACAATg                                                                                           >  4:312357/1‑149 (MQ=255)
                                                       ggTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATg                                                                                           <  7:213441/149‑1 (MQ=255)
                                                          aGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACg                                                                                                    >  1:8034/1‑137 (MQ=255)
                                                          aGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACg                                                                                                    <  2:8034/137‑1 (MQ=255)
                                                             gTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGATGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATg                                                                                           >  5:425169/1‑143 (MQ=255)
                                                             gTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGATGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATg                                                                                           <  6:425169/143‑1 (MQ=255)
                                                             gTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAg                                                                                     >  3:347206/1‑149 (MQ=255)
                                                                ggTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAggg                                                                                   <  5:278225/148‑1 (MQ=255)
                                                                 gTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAgggg                                                                                  >  6:154924/1‑148 (MQ=255)
                                                                        cATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGAtt                                                                          <  3:388779/149‑1 (MQ=255)
                                                                         aTCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGAttt                                                                         <  3:151089/149‑1 (MQ=255)
                                                                            ggTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTcc                                                                       >  4:88088/1‑148 (MQ=255)
                                                                             gTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCg                                                                      >  3:244163/1‑148 (MQ=255)
                                                                             gTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCg                                                                      >  1:99874/1‑148 (MQ=255)
                                                                               cGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCg                                                                      >  3:194107/1‑146 (MQ=255)
                                                                                 aaGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAAt                                                                               <  6:34938/135‑1 (MQ=255)
                                                                                 aaGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAAt                                                                               >  5:34938/1‑135 (MQ=255)
                                                                                      ccGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGaa                                                                    <  6:233964/141‑1 (MQ=255)
                                                                                      ccGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGaa                                                                    >  5:233964/1‑141 (MQ=255)
                                                                                        ggTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCa                                                           >  2:396563/1‑148 (MQ=255)
                                                                                               aCCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACgg                                                                                                   <  6:328925/101‑1 (MQ=255)
                                                                                               aCCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACgg                                                                                                   >  5:328925/1‑101 (MQ=255)
                                                                                                 ccaccaGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGt                                                 >  4:184482/1‑149 (MQ=255)
                                                                                                 ccaccaGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGt                                                 >  1:187576/1‑149 (MQ=255)
                                                                                                  caccaGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGtt                                                <  3:88088/149‑1 (MQ=255)
                                                                                                   accaGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGc                                                                                                                                            >  8:137373/1‑56 (MQ=255)
                                                                                                   accaGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGc                                                                                                                                            <  7:137373/56‑1 (MQ=255)
                                                                                                   accaGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGtt                                                >  6:214252/1‑148 (MQ=255)
                                                                                                   accaGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCACGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGtt                                                >  4:10531/1‑148 (MQ=255)
                                                                                                             cTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATg                                     >  1:118752/1‑149 (MQ=255)
                                                                                                                  tCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTg                                <  2:187576/149‑1 (MQ=255)
                                                                                                                  tCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTg                                <  3:326281/149‑1 (MQ=255)
                                                                                                                           tCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTAt                                      <  5:14710/134‑1 (MQ=255)
                                                                                                                           tCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTAt                                      >  6:14710/1‑134 (MQ=255)
                                                                                                                            cAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTgg                      >  5:29047/1‑149 (MQ=255)
                                                                                                                             aaCGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTgg                      <  5:147921/148‑1 (MQ=255)
                                                                                                                               cGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCa                   <  5:154924/149‑1 (MQ=255)
                                                                                                                               cGTTTGCGAAGCCATGCGCGTGCACGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCa                   <  3:10531/149‑1 (MQ=255)
                                                                                                                                gTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAg                  >  5:420837/1‑149 (MQ=255)
                                                                                                                                 tttGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGt                 <  4:194107/149‑1 (MQ=255)
                                                                                                                                     gcacgCCATGCACGCGCATGGCTTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCtt                                        >  3:179170/5‑122 (MQ=255)
                                                                                                                                     gcacgCCATGCACGCGCATGGCTTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCtt                                        <  4:179170/118‑1 (MQ=255)
                                                                                                                                     gcacgCCATGCACGCGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAg              >  7:78669/5‑148 (MQ=255)
                                                                                                                                     gcacgCCATGCACGCGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGc             <  8:78669/145‑1 (MQ=255)
                                                                                                                                       aaGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACa                            >  6:66587/1‑132 (MQ=255)
                                                                                                                                       aaGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACa                            <  5:66587/132‑1 (MQ=255)
                                                                                                                                        aGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTTCGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGt          <  4:244163/149‑1 (MQ=255)
                                                                                                                                            aTGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGcc                                                   <  2:90079/104‑1 (MQ=255)
                                                                                                                                            aTGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGcc                                                   >  1:90079/1‑104 (MQ=255)
                                                                                                                                            aTGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCt      >  7:316921/1‑149 (MQ=255)
                                                                                                                                            aTGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCt      <  4:347206/149‑1 (MQ=255)
                                                                                                                                              gcgcGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTg     >  1:153818/1‑148 (MQ=255)
                                                                                                                                               cgcgTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGa   <  2:153818/149‑1 (MQ=255)
                                                                                                                                                gcgTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGAc  >  6:37459/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CAGGGCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGAC  >  NC_002947/590935‑591227

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