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

TAATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATC  >  NC_000913/646598‑646891
                                                                                                                                                      |                                                                                                                                                   
tAATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGt                                                                                                                                                     <  2:438536/150‑1 (MQ=255)
 aaTATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTc                                                                                                                                                   <  2:923141/151‑1 (MQ=255)
     tttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTACAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGtgt                                                                                                                                                   >  1:675071/1‑145 (MQ=255)
      ttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAt                                                                                                                                              <  2:300121/150‑1 (MQ=255)
       tCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATg                                                                                                                                            <  2:949135/150‑1 (MQ=255)
             tttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGcccccc                                                                                                                                      >  1:1032385/1‑149 (MQ=255)
              ttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCGGAATCCCCCCtt                                                                                                                                     >  1:926570/1‑151 (MQ=255)
                aTCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCtttt                                                                                                                                   <  1:951239/150‑1 (MQ=255)
                aTCAAGTGCTTAAATAATTAAATCGGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACCCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCtttt                                                                                                                                   <  1:980666/150‑1 (MQ=255)
                 tCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCttttt                                                                                                                                  >  2:577369/1‑150 (MQ=255)
                  cAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGc                                                                                                                                <  1:141750/151‑1 (MQ=255)
                     gTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCaat                                                                                                                             <  2:471727/151‑1 (MQ=255)
                     gTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCACTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCTAATATCCCCCTTTTTgta                                                                                                                               >  1:262546/1‑148 (MQ=255)
                          tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATg                                                                                                                         <  2:687795/150‑1 (MQ=255)
                          tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCttttttgt                                                                                                                               >  2:96574/1‑141 (MQ=255)
                          tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCCTTTTTCAATAATGt                                                                                                                        >  2:583398/1‑151 (MQ=255)
                          tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGGCCCCCTTTTTCaaaaa                                                                                                                           >  2:1085351/1‑148 (MQ=255)
                             ataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTg                                                                                                                     <  1:318012/151‑1 (MQ=255)
                              taatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGTATGCCCCCTCTTTTTAATAATATTTgt                                                                                                                    >  1:457416/1‑150 (MQ=255)
                                atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCTATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTgcg                                                                                                                   <  2:851992/150‑1 (MQ=255)
                                         gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCt                                                                                                          <  1:19030/150‑1 (MQ=255)
                                           gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAAGCCCCCCTTTTTCCATAATGtttttgccccccgc                                                                                                           >  2:986012/1‑136 (MQ=255)
                                           gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCGCATGCCCCCCTTTTTCAATAATGTTTTTGCCTcccg                                                                                                            >  1:672668/1‑146 (MQ=255)
                                              aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGcc                                                                                                    <  2:405998/151‑1 (MQ=255)
                                               aaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGGCCCCCTTTTTTGATAAAGTTTGCGcccccccctggggac                                                                                                    >  1:685558/1‑137 (MQ=255)
                                               aaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATTTTTTTGCCCCCCGCTGTAGcc                                                                                                    >  2:927912/1‑150 (MQ=255)
                                                   ccAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAgcg                                                                                                <  2:570364/150‑1 (MQ=255)
                                                     aGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGTGGTTTCAATTGCCCCCTTTTTCCAAAAAGTTTTCGCCTCGCCCCGTAgt                                                                                                     >  1:315576/1‑142 (MQ=255)
                                                     aGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGGCCCCCTTTTTCTATAAATGttttgcctcccgccggt                                                                                                        >  1:807112/1‑127 (MQ=255)
                                                      ggTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAgg                                                                                           <  1:649629/151‑1 (MQ=255)
                                                        tAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAAGCCCCCCTTTTTCCATAAATTTTTTGCCTCCCGCCGGAGCCAGCGCCGGcc                                                                                         >  1:383037/1‑151 (MQ=255)
                                                        tAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCGCATGCCCCCCTTTTGTCATAATTTTTTCGCCTCccgcgg                                                                                                         >  1:201952/1‑134 (MQ=255)
                                                         aaGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCTGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCgcgc                                                                                                           <  1:959555/133‑1 (MQ=255)
                                                         aaGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGcc                                                                                         <  1:925400/150‑1 (MQ=255)
                                                         aaGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAAGGCCCCCTTTTTCAATAAAGTTTTCGCCCCCCGCCGTAGCCCGCGCCGGccc                                                                                        >  2:346717/1‑151 (MQ=255)
                                                          aGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAAGCCCCCCTTTTTTAATAAATTTTTCGCCCCGCGCCGTAGCCAGCGCAGGcccc                                                                                       >  2:883547/1‑151 (MQ=255)
                                                            gATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTTGGATGCCCCCCTTTTGTAATAAAGTTTTCGCCTCGCGCTGTAGACCGCGCAGGcccccc                                                                                     >  2:370351/1‑151 (MQ=255)
                                                               cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTc                                                                                   <  1:304303/150‑1 (MQ=255)
                                                               cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTc                                                                                   <  2:181921/150‑1 (MQ=255)
                                                                  tAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCGAATGCCCCCCTTTTTCAAAAATTTTTTTGCCCCCCCCTGTAGGCAGCGCACGccccccccccc                                                                                >  1:876639/1‑145 (MQ=255)
                                                                       cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgc                                                                           <  1:692097/150‑1 (MQ=255)
                                                                        aGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGCAATCCCCCCTTTTGTTATAATTTTTTCGCCTCCCGCTGGAGGCCGCGCACGGccccccccccccccccc                                                                          >  2:522663/1‑140 (MQ=255)
                                                                        aGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCTTTTTTTAATAATTTTTTCGCCTCGCGCTGTAGCCCGCGCCGGCCCCCCTCCCCAcccc                                                                           >  2:803451/1‑147 (MQ=255)
                                                                        aGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGATCGACTGGCCCCCCTTTGTTATAAATTTTTCCCCTCGAGCTCCAGCCACCGCAGGcgcccccccccccctc                                                                           >  1:830094/1‑139 (MQ=255)
                                                                           aCTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTGCAATAAATTTTTCGCCTCGCCCCGTAGCCAGCGCAGGCCCCCCCCCCCAcccccc                                                                         >  2:146885/1‑144 (MQ=255)
                                                                                   cTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCCTTTTTCAATAAAATTTTCCCCCCCCCCCGTAGCCACGGCACGCCCCCCCCCCCACCCCGCGCGGCCACGt                                                               >  1:731894/1‑150 (MQ=255)
                                                                                    tGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTCAATAATATTTGCGCCTCCCGCCGTAGGCCGCGCCGGCCCCCCCCCCCCCCGCCCGGTGCAACGtt                                                              >  1:577776/1‑150 (MQ=255)
                                                                                     ggACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGtt                                                              <  1:300436/149‑1 (MQ=255)
                                                                                        cTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAAGGCCCCCTTTTTTAATAAAATTTTCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCCCGCCACCGCGCGCTgggaaggtggg                                                           >  2:1026942/1‑139 (MQ=255)
                                                                                         tGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGGCCCCCTTTTTTAATAATGTTTTGGCCCCGCGCTGTAACCAGGGGAGGCCCCCCCCCCCACCCCGCGGTGCAAAGCTgg                                                            >  2:538120/1‑147 (MQ=255)
                                                                                          gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATc                                                        <  1:651700/150‑1 (MQ=255)
                                                                                          gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTTTATAAATTTTTCGCCTCCCGCTTGAGCCCGCGCCGGCCCCCCCCCCCCCcccgcggcgcgacgagcg                                                            >  1:817630/1‑134 (MQ=255)
                                                                                                 cgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGa                                                <  2:619916/151‑1 (MQ=255)
                                                                                                  gttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCCCATGCCCCCCTTTTTCAAAAAATTTTTTGCCCCCCCCCGTAGCCCGCGCAGGGCCCCCCCCCCACCCCCCGGTGCAACGCTGGGATCGcccgt                                                  >  1:175002/1‑145 (MQ=255)
                                                                                                     cgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCg                                             <  2:281635/150‑1 (MQ=255)
                                                                                                          cacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCa                                       <  1:144811/151‑1 (MQ=255)
                                                                                                               ccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGt                                     <  1:39595/148‑1 (MQ=255)
                                                                                                                cTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCGAATGGCCCCCTTTTTCCATAATTTTTTTGCCCCCCCCCGTTGCCAGCGCCGGCCCCCCCCCCCCCCCCCCGCTGCCACGCTGGTATCGCCCCTGGTCGCCACCCGTAgg                                  >  1:664686/1‑150 (MQ=255)
                                                                                                                  gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAg                                <  1:392332/150‑1 (MQ=255)
                                                                                                                   tCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGCATGCCCCCCTTTTTCAATTATTTTTTCGCCCCGCGGGGTAGCCCGCGCAGGCCCCCCCCCCCCCCCCCCGAAGGCACGGTGGTATTGCCGGTGATCGCCATCAGTAg                                   >  1:240063/1‑146 (MQ=255)
                                                                                                                   tCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGCAAGCCCCCTTTTTTCAATAAAGTTTTTGCCCCCCGCGGTAGCCAGCGCCGGCCCCCCCCCCCCCCCCCCGGAGCCAAGGTGGGATCGCCGCTGATCGCCCTCCGTAGGAg                                >  1:480774/1‑149 (MQ=255)
                                                                                                                     gAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAg                            <  2:373869/151‑1 (MQ=255)
                                                                                                                       aCTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGt                          <  1:167996/151‑1 (MQ=255)
                                                                                                                          gCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTAt                        <  1:115288/150‑1 (MQ=255)
                                                                                                                           ccGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCg                      <  2:802007/151‑1 (MQ=255)
                                                                                                                           ccGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCg                      <  2:21192/151‑1 (MQ=255)
                                                                                                                                 aTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAAt                <  2:416759/151‑1 (MQ=255)
                                                                                                                                   aaTCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGc              <  1:671553/151‑1 (MQ=255)
                                                                                                                                        agaTCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGtt         <  1:20539/151‑1 (MQ=255)
                                                                                                                                          aTCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCa       <  1:930407/151‑1 (MQ=255)
                                                                                                                                               cggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATc  >  1:127592/1‑151 (MQ=255)
                                                                                                                                                gggggTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATc  >  2:362470/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
TAATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATC  >  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: