Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 812,699 C→T L144L (CTG→TTG)  bioD → dethiobiotin synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913812,6990CT97.4% 248.5 / NA 76L144L (CTG→TTG) bioDdethiobiotin synthetase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (41/33);  minor base G (0/2);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACAGCAGCCTGCAGCTGGATTACGCAACAGTAAATCCTTACACCTTCGCAGAACCCACTTCGCCGCACATCATCAGCGCGCAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTT  >  NC_000913/812452‑812943
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
cAACAGCAGCCTGCAGCTGGATTACGCAACAGTAAATCCTTACACCTTCGCAGAACCCACTTCGCCGCACATCATCAGCGCGCAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACGCTGGCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTgg                                                                                                                                                                                                                                                   <  1:58220/251‑1 (MQ=255)
cAACAGCAGCCTGCAGCTGGATTACGCAACAGTAAATCCTTACACCTTCGCAGAACCCACTTCGCCGCACATCATCAGCGCGCAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTgg                                                                                                                                                                                                                                                   <  1:446525/251‑1 (MQ=255)
cAACAGCAGCCTGCAGCTGGATTACGCAACAGTAAATCCTTACACCTTCGCAGAACCCACTTCGCCGCACATCATCAGCGCGCAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTgg                                                                                                                                                                                                                                                   >  2:379233/1‑251 (MQ=255)
  aCAGCAGCCTGCAGCTGGATTACGCAACAGTAAATCCTTACACCTTCGCAGAACCCACTTCGCCGCACATCATCAGCGCGCAAGCGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGt                                                                                                                                                                                                                                                  >  1:38857/1‑250 (MQ=255)
         ccTGCAGCTGGATTACGCACCAGTAAATCCTTACACCTTCGCAGAACCCACTTCGCCGCACATCATCAGCGCGCAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGgtg                                                                                                                                                                                                                                          >  1:543859/1‑251 (MQ=255)
            gCAGCTGGATTACGCAACAGTAAATCCTTACACCTTCGCAGAACCCACTTCGCCGCACATCATCAGCGCGCAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGtgtg                                                                                                                                                                                                                                        >  1:406619/1‑250 (MQ=255)
              aGCTGGATTACGCAACAGTAAATCCTTACACCTTCGCAGAACCCACTTCGCCGCACATCATCAGCGCGCAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCCGATTGGGTAACACAGGAACAACTGCCGGGGATATTGGTAGTTGGGGTGaa                                                                                                                                                                                                                                      >  2:458154/1‑250 (MQ=255)
              aGCTGGATTACGCAACAGTAAATCCTTACACCTTCGCAGAACCCACTTCGCCGCACATCATCAGCGCGCAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATCTTGGTAGTTGGTGTGaa                                                                                                                                                                                                                                      >  2:92976/1‑250 (MQ=255)
                   gATTACGCAACAGTAAATCCTTACACCTTCGCAGAACCCACTTCGCCGCACATCATCAGCGCGCAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCgg                                                                                                                                                                                                                                >  1:521473/1‑251 (MQ=255)
                    aTTACGCAACAGTAAATCCTTACACCTTCGCAGAACCCACTTCGCCGCACATCATCAGCGCGCAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGc                                                                                                                                                                                                                               <  1:56562/251‑1 (MQ=255)
                        cGCAACAGTAAATCCTTACACCTTCGCAGAACCCACTTCGCCGCACATCATCAGCGCGCAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTa                                                                                                                                                                                                                           <  1:559765/251‑1 (MQ=255)
                          cAACAGTAAATCCTTACACCTTCGCAGAACCCACTTCGCCGCACATCATCAGCGCGCAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTAtt                                                                                                                                                                                                                         >  2:93814/1‑251 (MQ=255)
                                 aaaTCCTTACACCTTCGCAGAACCCACTTCGCCGCACATCATCAGCGCGCAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCAcg                                                                                                                                                                                                                  <  2:556046/251‑1 (MQ=255)
                                         acacCTTCGCAGAACCCACATCGCCGCACATCATCAGCGCGCAAGAGGGCAGACCGATAGCATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAAATGCCGGTGATATTGGTAGTTGGTGTGAAACTAGGCTGTATTAATCACGCGATGTTg                                                                                                                                                                                                          >  2:48760/1‑251 (MQ=255)
                                              ttCGCAGAACCCACTTCGCCGCACATCATCAGCGCGCAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGc                                                                                                                                                                                                     >  1:449083/1‑251 (MQ=255)
                                               tCGCAGATCCCACTTCGCCGCACATCATCAGCGCGCAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGATACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGc                                                                                                                                                                                                     <  2:121844/250‑1 (MQ=255)
                                                cGCAGAACCCACTTCGCCGCACATCATCAGCGCGCAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGcac                                                                                                                                                                                                   <  2:449083/251‑1 (MQ=255)
                                                cGCAGAACCCACTTCGCCGCACATCATCAGCGCGCAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGcac                                                                                                                                                                                                   >  1:503683/1‑251 (MQ=255)
                                                    gAACCCACTTCGCCGCACATCATCAGCGCGCAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAggg                                                                                                                                                                                               >  1:258779/1‑250 (MQ=255)
                                                              cgccgcACATCATCAGCGCGCAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAacac                                                                                                                                                                                     <  2:524398/251‑1 (MQ=255)
                                                                 cgcACATCATCAGCGCGCAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTAAGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACAAGcc                                                                                                                                                                                  >  2:187775/1‑251 (MQ=255)
                                                                  gcACATCATCAGCGCGCAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGCTATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCg                                                                                                                                                                                 >  1:294954/1‑251 (MQ=255)
                                                                      atcatcAGCGCGCAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGgact                                                                                                                                                                             >  1:596383/1‑251 (MQ=255)
                                                                        catcaGCGCGCAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGgactg                                                                                                                                                                            <  1:248076/250‑1 (MQ=255)
                                                                            aGCGTGCAAGAGGGCATACCGATAGAATCTGTGGTAATTAGCGCCGGATTACGCGCGCTTGAACAACAGACTGACTGGTTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGTCGCGTGTTTACTCTTTCACTTTTGCAGATTTTGTAACACAGTATCACACGCCGGTGATATTGTTAGTTGGTTTGAAACTCGGCTGTATTAATGACGCGATGTTCACTGCACTGGttattaaaatcgtaggactgacgct                                                                                                                                                                       <  2:396914/251‑25 (MQ=255)
                                                                               gcgcAAGAGGGCAGACCGATAGAATCATTGGGAATGAGAGCCGGATTACGCGCGCTTGACAAACAGGCTGACTGGGTGTTAGTGGAAGGTGCGGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACGCAGGAACACCTGCCGGTGATAGTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACGGACGCTGGc                                                                                                                                                                    <  2:303716/251‑1 (MQ=255)
                                                                                 gcAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCgg                                                                                                                                                                  <  2:40736/251‑1 (MQ=255)
                                                                                 gcAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGAGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCgg                                                                                                                                                                  >  2:37568/1‑251 (MQ=255)
                                                                                    agagGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGGGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCgggat                                                                                                                                                               >  1:299133/1‑249 (MQ=255)
                                                                                      agGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTgg                                                                                                                                                             <  2:556258/251‑1 (MQ=255)
                                                                                         gCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGGGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGggtgggtg                                                                                                                                                          >  1:363413/1‑251 (MQ=255)
                                                                                          cAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTgg                                                                                                                                                         <  2:482874/251‑1 (MQ=255)
                                                                                           aGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGTCGGCTGGTTTACGCCGCTTTCTGACACTTTCTCTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGTCGGGTTGGGTGGc                                                                                                                                                        <  2:252025/251‑1 (MQ=255)
                                                                                           aGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGc                                                                                                                                                        >  1:268854/1‑251 (MQ=255)
                                                                                             aCCGATAGAATCATTGGTAATGAGCGCCGTATTACGAGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGCCACTTTCACTTTTGCGGAGTGGGTAACACCGGAACACCGGCAGGGGATTTTGGTAGTTGGTGTGAAACTCGGCTGTCTTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGGTGGGTGGCgg                                                                                                                                                      >  2:61962/1‑250 (MQ=255)
                                                                                              ccGATAGAATCATTGGTAATTAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTGCGCCGTTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGACCAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTGGACTGCACATGTAATACAACACGCCGGACTGACTCTGGTGGGTTGGGTGGCGaa                                                                                                                                                     <  2:354111/251‑1 (MQ=255)
                                                                                                   aGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATg                                                                                                                                                <  1:320741/251‑1 (MQ=255)
                                                                                                     aaTCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGtt                                                                                                                                              <  1:350260/251‑1 (MQ=255)
                                                                                                                         ggATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCa                                                                                                                          <  1:248750/251‑1 (MQ=255)
                                                                                                                         ggATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGGTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGGTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTc                                                                                                                           >  1:303219/1‑250 (MQ=255)
                                                                                                                                 cgcgcTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGGAATACAACACGCCGGACTGACTCTGGCGGGGTGGGTGGCGAACGATGGTACGCCCCCGGGAAAACGGCACGCTggat                                                                                                                  >  2:564701/1‑248 (MQ=255)
                                                                                                                                     cTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAatatat                                                                                                              >  1:144343/1‑251 (MQ=255)
                                                                                                                                       tGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAatatat                                                                                                              <  1:269291/249‑1 (MQ=255)
                                                                                                                                                   tGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCAccc                                                                                                >  2:94569/1‑251 (MQ=255)
                                                                                                                                                   tGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGATTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCAccc                                                                                                >  1:559655/1‑251 (MQ=255)
                                                                                                                                                     aCTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGc                                                                                              >  2:568068/1‑251 (MQ=255)
                                                                                                                                                        gggTGTTAGTGGAAGTTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGTTGATATTGGTAGTTGGTGTGAAACTCGGCTTTATTAATCACGCGATGTTGCCTGCACAGGTAATACAACACGCCGGACTGACGCTGGCGGGTTGGGTGGCGAACGATGTTACGCGTCCGGGATAACGTCACGCTGAATATATGACCACGCTCACCCGCATg                                                                                           <  1:71981/251‑1 (MQ=255)
                                                                                                                                                            gtTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTc                                                                                       <  2:543859/251‑1 (MQ=255)
                                                                                                                                                               aGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCcg                                                                                    <  2:234145/251‑1 (MQ=255)
                                                                                                                                                                            ggcggcTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTgg                                                                        <  1:41735/250‑1 (MQ=255)
                                                                                                                                                                             gcggcTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGGCCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGga                                                                      <  2:406619/251‑1 (MQ=255)
                                                                                                                                                                               ggcTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTggg                                                                       >  1:189755/1‑248 (MQ=255)
                                                                                                                                                                                  tGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGCACTGACTCTGGCGGGTTGGGTGGCGAAGGATTTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCCCACCCGCATGATGCCCGCCCCGCCGCTggaag                                                                     >  2:237710/1‑244 (MQ=255)
                                                                                                                                                                                  tGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATAGTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGgagaga                                                                  <  1:220076/250‑1 (MQ=255)
                                                                                                                                                                                   ggTTTACGCCGCATTCTGACACTTTCACTTGTGCAGATTGGGTAACACAGGAACAACTTCCGGTGATATTGGTAGTTGGTGGGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGGCACGCTGAATAGATGACCACGCTCACCCGCATGATTCCCGCGACGCTGCTGGGAGAGAt                                                                 <  2:237417/250‑1 (MQ=255)
                                                                                                                                                                                              cTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGc                                                     >  1:450523/1‑251 (MQ=255)
                                                                                                                                                                                               tttCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAGCACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGc                                                     >  2:2426/1‑250 (MQ=255)
                                                                                                                                                                                                 tCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAg                                                   <  2:194396/250‑1 (MQ=255)
                                                                                                                                                                                                  cTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGa                                                  >  1:5309/1‑250 (MQ=255)
                                                                                                                                                                                                    gACACTTTCAATTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGaaa                                                >  1:451702/1‑250 (MQ=255)
                                                                                                                                                                                                      cactttcactttTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATc                                             <  2:299573/251‑1 (MQ=255)
                                                                                                                                                                                                      cactttcactgtTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATc                                             <  1:439962/251‑1 (MQ=255)
                                                                                                                                                                                                         tttcactttTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGCTTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATcc                                            >  2:203853/1‑249 (MQ=255)
                                                                                                                                                                                                           tcactttTGCAGATTGGGTAACACAGGAACAACTGACGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACCGGTAATACAACACGCCGGACTGACTCTGGCGGGGTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGATCACCCGCATGATTCCCGCGCCGCTCCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGaa                                        >  2:458513/1‑250 (MQ=255)
                                                                                                                                                                                                             actttTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGaaa                                       >  1:144648/1‑250 (MQ=255)
                                                                                                                                                                                                                tttGCAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATTATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAAt                                     <  1:187775/249‑1 (MQ=255)
                                                                                                                                                                                                                    cAGATTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCCGAAAATCCAGAAAATGCGGc                                >  1:269346/1‑250 (MQ=255)
                                                                                                                                                                                                                       aTTGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCCAcc                            >  2:181162/1‑251 (MQ=255)
                                                                                                                                                                                                                        ttGGGTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCg                           >  2:612054/1‑251 (MQ=255)
                                                                                                                                                                                                                          gggTAACACAGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCgg                          <  2:503683/250‑1 (MQ=255)
                                                                                                                                                                                                                                cacaGGAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTAc                   >  1:254966/1‑251 (MQ=255)
                                                                                                                                                                                                                                     gAACAACTGCCGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATaaa              <  1:92976/251‑1 (MQ=255)
                                                                                                                                                                                                                                               ccGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTc                                                                                                                           >  2:46887/1‑133 (MQ=255)
                                                                                                                                                                                                                                               ccGGTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTc                                                                                                                           <  1:46887/133‑1 (MQ=255)
                                                                                                                                                                                                                                                 ggTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATGTGAACACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCttgt   >  1:653219/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  gTGATATTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACGGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCttgtt  <  1:317607/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
CAACAGCAGCCTGCAGCTGGATTACGCAACAGTAAATCCTTACACCTTCGCAGAACCCACTTCGCCGCACATCATCAGCGCGCAAGAGGGCAGACCGATAGAATCATTGGTAATGAGCGCCGGATTACGCGCGCTTGAACAACAGGCTGACTGGGTGTTAGTGGAAGGTGCTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTT  >  NC_000913/812452‑812943

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: