Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913646,7470TG43.5% 101.2 / 29.1 70T255P (ACC→CCC) citG2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase
Reads supporting (aligned to +/- strand):  ref base T (4/35);  new base G (30/0);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.47e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATC  >  NC_000913/646598‑646891
                                                                                                                                                      |                                                                                                                                                   
tAATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGt                                                                                                                                                     <  2:438536/150‑1 (MQ=255)
 aaTATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTc                                                                                                                                                   <  2:923141/151‑1 (MQ=255)
     tttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTACAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGtgt                                                                                                                                                   >  1:675071/1‑145 (MQ=255)
      ttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAt                                                                                                                                              <  2:300121/150‑1 (MQ=255)
       tCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑g                                                                                                                                            <  2:949135/150‑1 (MQ=255)
             tttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTC‑GGAT‑Gcccccc                                                                                                                                      >  1:1032385/1‑149 (MQ=255)
              ttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGGCGGAAT‑CCCCCCtt                                                                                                                                     >  1:926570/1‑151 (MQ=255)
                aTCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCtttt                                                                                                                                   <  1:951239/150‑1 (MQ=255)
                aTCAAGTGCTTAAATAATTAAATCGGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACCCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCtttt                                                                                                                                   <  1:980666/150‑1 (MQ=255)
                 tCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTC‑GAAT‑GCCCCCttttt                                                                                                                                  >  2:577369/1‑150 (MQ=255)
                  cAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGc                                                                                                                                <  1:141750/151‑1 (MQ=255)
                     gTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCaat                                                                                                                             <  2:471727/151‑1 (MQ=255)
                     gTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCACTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTC‑TAATATCCCCCTTTTTgta                                                                                                                               >  1:262546/1‑148 (MQ=255)
                          tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATg                                                                                                                         <  2:687795/150‑1 (MQ=255)
                          tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTT‑GAAT‑GCCCCCttttttgt                                                                                                                               >  2:96574/1‑141 (MQ=255)
                          tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTT‑GAAT‑GCCCCCCTTTTTCAATAATGt                                                                                                                        >  2:583398/1‑151 (MQ=255)
                          tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTC‑GGAT‑GGCCCCCTTTTTCaaaaa                                                                                                                           >  2:1085351/1‑148 (MQ=255)
                             ataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTg                                                                                                                     <  1:318012/151‑1 (MQ=255)
                              taatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTC‑GTAT‑GCCCCCTCTTTTTAATAATATTTgt                                                                                                                    >  1:457416/1‑150 (MQ=255)
                                atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCTATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTgcg                                                                                                                   <  2:851992/150‑1 (MQ=255)
                                         gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCt                                                                                                          <  1:19030/150‑1 (MQ=255)
                                           gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTC‑GAAA‑GCCCCCCTTTTTCCATAATGtttttgccccccgc                                                                                                           >  2:986012/1‑136 (MQ=255)
                                           gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGGC‑GCAT‑GCCCCCCTTTTTCAATAATGTTTTTGCCTcccg                                                                                                            >  1:672668/1‑146 (MQ=255)
                                              aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGcc                                                                                                    <  2:405998/151‑1 (MQ=255)
                                               aaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTC‑GAAT‑GGCCCCCTTTTTTGATAAAGTTTGCGcccccccctggggac                                                                                                    >  1:685558/1‑137 (MQ=255)
                                               aaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTC‑GAAT‑GCCCCCCTTTTTCAATAATTTTTTTGCCCCCCGCTGTAGcc                                                                                                    >  2:927912/1‑150 (MQ=255)
                                                   ccAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAgcg                                                                                                <  2:570364/150‑1 (MQ=255)
                                                     aGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGTGGTTTC‑AATT‑GCCCCCTTTTTCCAA‑AAAGTTTTCGCCTCGCCCCGTAgt                                                                                                     >  1:315576/1‑142 (MQ=255)
                                                     aGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTC‑GGATGGCCCCCTTTTTCTATAAATGttttgcctcccgccggt                                                                                                        >  1:807112/1‑127 (MQ=255)
                                                      ggTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑gg                                                                                           <  1:649629/151‑1 (MQ=255)
                                                        tAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTC‑GAAA‑GCCCCCCTTTTTCCATAAATTTTTTGCCTCCCGCCGGAGCCAGCGCC‑GGcc                                                                                         >  1:383037/1‑151 (MQ=255)
                                                        tAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGGC‑GCAT‑GCCCCCCTTTTGTCATAATTTTTTCGCCTCccgcgg                                                                                                         >  1:201952/1‑134 (MQ=255)
                                                         aaGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCTGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCgcgc                                                                                                           <  1:959555/133‑1 (MQ=255)
                                                         aaGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑GGcc                                                                                         <  1:925400/150‑1 (MQ=255)
                                                         aaGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTC‑GAAA‑GGCCCCCTTTTTCAATAAAGTTTTCGCCCCCCGCCGTAGCCCGCGCC‑GGccc                                                                                        >  2:346717/1‑151 (MQ=255)
                                                          aGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTC‑GAAA‑GCCCCCCTTTTTTAATAAATTTTTCGCCCCGCGCCGTAGCCAGCGCA‑GGcccc                                                                                       >  2:883547/1‑151 (MQ=255)
                                                            gATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTT‑GGAT‑GCCCCCCTTTTGTAATAAAGTTTTCGCCTCGCGCTGTAGACCGCGCA‑GGcccccc                                                                                     >  2:370351/1‑151 (MQ=255)
                                                               cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑GGCCCCCCTc                                                                                   <  1:304303/150‑1 (MQ=255)
                                                               cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑GGCCCCCCTc                                                                                   <  2:181921/150‑1 (MQ=255)
                                                                  tAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGGC‑GAAT‑GCCCCCCTTTTTCAAAAATTTTTTTGCCCCCCCCTGTAGGCAGCGCA‑CGccccccccccc                                                                                >  1:876639/1‑145 (MQ=255)
                                                                       cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACcgc                                                                           <  1:692097/150‑1 (MQ=255)
                                                                        aGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTC‑GCAA‑TCCCCCCTTTTGTTATAATTTTTTCGCCTCCCGCTGGAGGCCGCGCACGGccccccccccccccccc                                                                          >  2:522663/1‑140 (MQ=255)
                                                                        aGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTC‑GAAT‑GCCCCCTTTTTTTAATAATTTTTTCGCCTCGCGCTGTAGCCCGCGCC‑GGCCCCCCTCCCCAcccc                                                                           >  2:803451/1‑147 (MQ=255)
                                                                        aGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGATC‑GACT‑GGCCCCCCTTTGTTATAAATTTTTCCCCTCGAGCTCCAGCCACCGCA‑GGcgcccccccccccctc                                                                           >  1:830094/1‑139 (MQ=255)
                                                                           aCTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTC‑GAAT‑GCCCCCCTTTTGCAATAAATTTTTCGCCTCGCCCCGTAGCCAGCGCA‑GGCCCCCCCCCCCAcccccc                                                                         >  2:146885/1‑144 (MQ=255)
                                                                                   cTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTT‑GAAT‑GCCCCCCTTTTTCAATAAAATTTTCCCCCCCCCCCGTAGCCACGGCA‑CGCCCCCCCCCCCACCCCGCGCGGCCACGt                                                               >  1:731894/1‑150 (MQ=255)
                                                                                    tGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTC‑GGAT‑GCCCCCCTTTTTCAATAATATTTGCGCCTCCCGCCGTAGGCCGCGCC‑GGCCCCCCCCCCCCCCGCCCGGTGCAACGtt                                                              >  1:577776/1‑150 (MQ=255)
                                                                                     ggACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGtt                                                              <  1:300436/149‑1 (MQ=255)
                                                                                        cTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTC‑GAAA‑GGCCCCCTTTTTTAATAAAATTTTCGCCTCGCGCTGTAGCCAGCGCA‑GGCCCCCCCCGCCACCGCGCGCTgggaaggtggg                                                           >  2:1026942/1‑139 (MQ=255)
                                                                                         tGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTC‑GAAT‑GGCCCCCTTTTTTAATAATGTTTTGGCCCCGCGCTGTAACCAGGGGA‑GGCCCCCCCCCCCACCCCGCGGTGCAAAGCTgg                                                            >  2:538120/1‑147 (MQ=255)
                                                                                          gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATc                                                        <  1:651700/150‑1 (MQ=255)
                                                                                          gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTC‑GGAT‑GCCCCCCTTTTTTTATAAATTTTTCGCCTCCCGCTTGAGCCCGCGCC‑GGCCCCCCCCCCCCCcccgcggcgcgacgagcg                                                            >  1:817630/1‑134 (MQ=255)
                                                                                                 cgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGa                                                <  2:619916/151‑1 (MQ=255)
                                                                                                  gttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTC‑CCAT‑GCCCCCCTTTTTCAAAAAATTTTTTGCCCCCCCCCGTAGCCCGCGCA‑GGGCCCCCCCCCCACCCCCCGGTGCAACGCTGGGATCGcccgt                                                  >  1:175002/1‑145 (MQ=255)
                                                                                                     cgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCg                                             <  2:281635/150‑1 (MQ=255)
                                                                                                          cacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCa                                       <  1:144811/151‑1 (MQ=255)
                                                                                                               ccTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGt                                     <  1:39595/148‑1 (MQ=255)
                                                                                                                cTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGGC‑GAAT‑GGCCCCCTTTTTCCATAATTTTTTTGCCCCCCCCCGTTGCCAGCGCC‑GGCCCCCCCCCCCCCCCCCCGCTGCCACGCTGGTATCGCCCCTGGTCGCCACCCGTAgg                                  >  1:664686/1‑150 (MQ=255)
                                                                                                                  gTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAg                                <  1:392332/150‑1 (MQ=255)
                                                                                                                   tCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTC‑GCAT‑GCCCCCCTTTTTCAATTATTTTTTCGCCCCGCGGGGTAGCCCGCGCA‑GGCCCCCCCCCCCCCCCCCCGAAGGCACGGTGGTATTGCCGGTGATCGCCATCAGTAg                                   >  1:240063/1‑146 (MQ=255)
                                                                                                                   tCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGGTC‑GCAA‑GCCCCCTTTTTTCAATAAAGTTTTTGCCCCCCGCGGTAGCCAGCGCC‑GGCCCCCCCCCCCCCCCCCCGGAGCCAAGGTGGGATCGCCGCTGATCGCCCTCCGTAGGAg                                >  1:480774/1‑149 (MQ=255)
                                                                                                                     gAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAg                            <  2:373869/151‑1 (MQ=255)
                                                                                                                       aCTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGt                          <  1:167996/151‑1 (MQ=255)
                                                                                                                          gCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTAt                        <  1:115288/150‑1 (MQ=255)
                                                                                                                           ccGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCg                      <  2:802007/151‑1 (MQ=255)
                                                                                                                           ccGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCg                      <  2:21192/151‑1 (MQ=255)
                                                                                                                                 aTAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAAt                <  2:416759/151‑1 (MQ=255)
                                                                                                                                   aaTCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGc              <  1:671553/151‑1 (MQ=255)
                                                                                                                                        agaTCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGtt         <  1:20539/151‑1 (MQ=255)
                                                                                                                                          aTCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCa       <  1:930407/151‑1 (MQ=255)
                                                                                                                                               cgg‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATc  >  1:127592/1‑151 (MQ=255)
                                                                                                                                                gg‑gggTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATc  >  2:362470/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
TAATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGG‑GGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATC  >  NC_000913/646598‑646891

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