Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 349,717 T→C 8.6% intergenic (+145/‑295) prpB → / → prpC 2‑methylisocitrate lyase/2‑methylcitrate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913349,7170TC 8.6% 201.4 / 8.3 97intergenic (+145/‑295)prpB/prpC2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base T (40/45);  new base C (5/4);  total (45/52)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

ATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCG  >  NC_000913/349569‑349865
                                                                                                                                                    |                                                                                                                                                    
aTAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGt                                                                                                                                                      >  2:1104101/1‑149 (MQ=255)
aTAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGt                                                                                                                                                      >  2:1139841/1‑149 (MQ=255)
 tAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTg                                                                                                                                                     >  1:59860/1‑149 (MQ=255)
  gaaaaaCGCCCGTTGATTGTATTCGACAGCCGATGCCTGCTGCGTCGCTGCCGCGACTTACCAGGCCTACGAGGTGCACGGAACTATCTGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGa                                                                                                                                                    <  2:62227/148‑1 (MQ=18)
   aaagaCGCCCGTTGATTGTATTCGGCGGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGAc                                                                                                                                                   <  2:905533/149‑1 (MQ=255)
    aaaaCGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCGCTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACt                                                                                                                                                  <  2:1077285/149‑1 (MQ=255)
       aCGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGGCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTAt                                                                                                                                               <  1:202010/149‑1 (MQ=255)
             gTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGGGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGcc                                                                                                                                         >  1:1099362/1‑149 (MQ=255)
             gTTGATTGTATTCAACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGcc                                                                                                                                         >  1:268869/1‑149 (MQ=255)
                 aTTGTCTTAGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGAGGCGCCGCGGGCGCGACTTAGCAGACCTACg                                                                                                                                     >  2:567601/1‑149 (MQ=11)
                  ttGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTACCAGGccaccgg                                                                                                                                    >  2:34598/1‑144 (MQ=255)
                   tGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGGTGCGCCGCGGCAGTGGATGATCAGGCCTACtgg                                                                                                                                   >  1:512219/1‑149 (MQ=17)
                   tGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACggg                                                                                                                                   <  1:921944/149‑1 (MQ=255)
                   tGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACggg                                                                                                                                   >  1:132061/1‑149 (MQ=255)
                      aTTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTg                                                                                                                                <  1:486907/149‑1 (MQ=255)
                           acctccgatgcccgatgcCGCGGTGACCCTACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTg                                                                                                                           <  2:1317470/145‑1 (MQ=14)
                           aCAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTg                                                                                                                           >  1:569939/1‑149 (MQ=255)
                             aGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGaa                                                                                                                         >  1:216073/1‑149 (MQ=255)
                             aGCCGATGCCCGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGaa                                                                                                                         <  2:777071/149‑1 (MQ=255)
                                  aTGCCTGATGCGTCGCTGACGCGACTTATCCGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTa                                                                                                                    <  2:1079498/149‑1 (MQ=255)
                                   tGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGACCGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAg                                                                                                                   <  1:719432/149‑1 (MQ=255)
                                       tGATGCGTCGCTGACGCGACTTATCAGGCCTAGGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCg                                                                                                               <  1:609633/149‑1 (MQ=255)
                                         aTGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCACATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGa                                                                                                             <  2:1107678/149‑1 (MQ=255)
                                            cGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATaa                                                                                                          >  2:452483/1‑149 (MQ=255)
                                            cGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACCTATCAGGCCTACGCGGTGCACTGAACTGTAGGTAGGATaa                                                                                                          >  2:663272/1‑149 (MQ=255)
                                            cGTCGCCGACGCGACTTATCGGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATaa                                                                                                          <  2:422216/149‑1 (MQ=255)
                                              tCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTCCGGGGGGCACTGAACTGTAGGTCGGATAAGa                                                                                                        >  1:786340/1‑149 (MQ=255)
                                              tCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTgaactgaact                                                                                                                       <  2:1314419/134‑1 (MQ=37)
                                              tCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTgaactgaact                                                                                                                       >  1:1314419/1‑134 (MQ=255)
                                               cGCTGCCCCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGAc                                                                                                       <  1:1298247/149‑1 (MQ=255)
                                                   gACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACgg                                                                                                                                    >  2:1182318/1‑116 (MQ=37)
                                                   gACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACgg                                                                                                                                    <  1:1182318/116‑1 (MQ=25)
                                                      gcgACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATGCGACGctatgc                                                                                                >  1:778225/1‑144 (MQ=255)
                                                              tCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACGTATCAGGCCTACGGGGGGGACTGAACTGGAGGGCGGATCGGACGCATAGCGGCGCGTc                                                                                        >  2:89319/1‑149 (MQ=14)
                                                                 ggCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGa                                                                                     <  2:246915/149‑1 (MQ=40)
                                                                     tACGAGGTGCACTGAACTGTAGGTCAGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTGGGATAAGACGCATAGGGTCGCATCCGACAAc                                                                                 >  1:677144/1‑149 (MQ=255)
                                                                      aCGAGGGGCACCGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTATGTCGGATAAGACGCATAGCGTCGCATCCGACAAcc                                                                                <  1:71360/149‑1 (MQ=255)
                                                                       cGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACCGAACTGTAGGGCGGATAAGACGCATAGCGCCGCATCCCACAAcat                                                                               >  2:923248/1‑147 (MQ=255)
                                                                       cGAGGTCCACTGAACCGGAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCg                                                                               <  2:1325442/149‑1 (MQ=255)
                                                                         aGGTGCACTGGTCTGTTGGTCGGAGAAGCCTGCTGTCGTCGTATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAt                                                                             <  2:386309/149‑1 (MQ=11)
                                                                         aGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAt                                                                             <  1:1104101/149‑1 (MQ=39)
                                                                         aGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAt                                                                             <  1:1139841/149‑1 (MQ=40)
                                                                         aGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCCCTGCACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCAAt                                                                             >  2:538942/1‑149 (MQ=255)
                                                                           gggCACGGAACTGTAGGTCGGAGAAGACCGAGGGCGTCGCATCCGCCAACCGATGCCTGATGCGCCGCTGACGGGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGc                                                                           <  2:216785/147‑1 (MQ=11)
                                                                           gTGCACTGAACTGTAGTTCGGATAAGACGCATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGc                                                                           <  2:677144/149‑1 (MQ=17)
                                                                            ggCACTGAACGGTAGGTGGGATAAGCGGGGGGGCGTCGCACACGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAAATGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGcc                                                                          <  2:358868/148‑1 (MQ=11)
                                                                              cACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCCCTGAACTGTAGGCCGGATAAGACGCATAGCGTCGAATCCGCCAACCgtcgccag                                                                        >  2:223331/1‑142 (MQ=11)
                                                                               ccTGAACTGTACGGCGAATAAGACGGATGGCCTCGCCTCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGa                                                                       <  2:1125875/148‑1 (MQ=11)
                                                                               aCTGAACTGGAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGa                                                                       <  2:1343619/149‑1 (MQ=37)
                                                                                 tGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACaa                                                                                  >  1:362663/1‑136 (MQ=40)
                                                                                 tGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATg                                                                     >  1:685439/1‑149 (MQ=38)
                                                                                 tGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACaa                                                                                  <  2:362663/136‑1 (MQ=18)
                                                                                 tGAACTGTAGGTCGGATAAGACGGATGGCGTCACATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTACACTGAACTGTAGGGCGGGGAAGACGGATGGCGACGCAACCGACGAGCGACGCCTGATg                                                                     >  2:969886/1‑149 (MQ=11)
                                                                                      tGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGCCAACCGGTGCCTGGAGcgccg                                                                >  2:1039308/1‑149 (MQ=255)
                                                                                      tGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccg                                                                >  2:1039154/1‑149 (MQ=38)
                                                                                      tGTAGGTCGGATAAGACGCATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccg                                                                <  2:1195775/149‑1 (MQ=17)
                                                                                      gggcggacGGATAAGACGGAGTGCGTCGCATCCGACAGCCGATGCGTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccg                                                                <  2:71497/142‑1 (MQ=11)
                                                                                          ggtctaaTAAGCCGGATGGCAGCGCATCCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGa                                                            <  2:264557/143‑1 (MQ=11)
                                                                                           gTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCAGCGGAc                                                           >  2:1144347/1‑149 (MQ=255)
                                                                                            tCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGGGTCGCATCGGCCAACAGATGCCTGCTGACCCACTGagg                                                          >  1:180258/1‑147 (MQ=11)
                                                                                            tCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAcg                                                          >  1:836496/1‑149 (MQ=38)
                                                                                             cGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTGGGATAAGACGCATAGCGTCGCATGCGACAACCGATGCCTGATTCGAAGCTGAcgc                                                         >  2:1293704/1‑149 (MQ=11)
                                                                                              ggATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGCCAACAGCTGCATGATACGCAGCTAAagcg                                                        >  2:326402/1‑149 (MQ=11)
                                                                                               gATAAGCCGGATGGCGTCGCATCCGACAACCGTTGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGTGCCGCTGACGCGa                                                       <  1:922246/149‑1 (MQ=255)
                                                                                               gATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGTTGGGCCGCTCCTTTTAGGGAGCAGTCCTACGTGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGa                                                       <  1:593311/149‑1 (MQ=11)
                                                                                               gATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGa                                                                                                        >  1:1301632/1‑100 (MQ=255)
                                                                                               gATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGa                                                                                                        <  2:1301632/100‑1 (MQ=18)
                                                                                               gATAAGACGGATGGCGGCGCACCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGa                                                       <  1:577377/149‑1 (MQ=255)
                                                                                               caTAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGGCTGCTTTGCCCCAAATTTTCGCATTCGACAACCGATGCGCGATGCGCCGGTGACGCGa                                                       <  1:991519/148‑1 (MQ=11)
                                                                                                aTAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGTTGCCTGCTGCCCCGCGGACGCGAc                                                      >  1:19648/1‑149 (MQ=255)
                                                                                                aTAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCCCGGACGCCAc                                                      >  2:133304/1‑149 (MQ=255)
                                                                                                 tatgACGGCGGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACt                                                     <  2:333111/146‑1 (MQ=16)
                                                                                                  aaGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGCTGCGCCGCTGACGCGCCtt                                                    >  1:1353252/1‑149 (MQ=255)
                                                                                                     aCGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGCCTTATc                                                 >  1:339795/1‑149 (MQ=38)
                                                                                                     aCGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATc                                                 <  1:623328/149‑1 (MQ=37)
                                                                                                     aCGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATc                                                 >  2:1083034/1‑149 (MQ=38)
                                                                                                      cGGATTGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCa                                                <  2:160030/149‑1 (MQ=21)
                                                                                                      cGGATGGCGTCGCATNCGANANCNCNTGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCa                                                <  1:938946/149‑1 (MQ=255)
                                                                                                      cGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTCCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCGGCTGACGCGAGTTATCa                                                >  1:1179259/1‑149 (MQ=255)
                                                                                                       ggATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAg                                               >  1:340526/1‑149 (MQ=35)
                                                                                                       ggATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAg                                               >  1:537456/1‑149 (MQ=37)
                                                                                                       agATGGCGTCGCTGTCGCCAGGCGATGCCGGGTGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAg                                               <  2:165628/148‑1 (MQ=11)
                                                                                                                cGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACggg                                      <  2:132061/149‑1 (MQ=32)
                                                                                                                    tCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGc                                  >  1:734176/1‑149 (MQ=33)
                                                                                                                     ccGACAGCGGATGCCGGGTGCGCCGCGGACGTGACTTAGCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAGGCCGGATGGCGTCGCATCCGTCAACCGATGCCTGATGCGCCGCTGACGCGACTTTTCAGGCCTACGGGGTGCa                                 <  2:1074804/149‑1 (MQ=11)
                                                                                                                      cGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAgg                                              >  2:843213/1‑135 (MQ=33)
                                                                                                                       gACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCGACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACt                               <  1:452483/149‑1 (MQ=21)
                                                                                                                          aaCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGaa                            <  2:833899/149‑1 (MQ=35)
                                                                                                                               aTGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACACGAGTTACCCGGCCTCCGGGGGGCACTGAACTGTa                       >  2:34378/1‑149 (MQ=11)
                                                                                                                                tGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAg                      <  1:541227/149‑1 (MQ=35)
                                                                                                                                    ggaggcgcCGCTGACGGGCCTTATCAGGCCTACGGGGTGCAGTGGACTGTAGGTGCGAAAAGCCGGATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCg                  <  2:524693/145‑1 (MQ=11)
                                                                                                                                    cgATGCGCCGCTGACGTGACTTATCAGGCCCACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCg                  <  2:59860/148‑1 (MQ=25)
                                                                                                                                     gATGCGCCGGGGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCgg                 <  2:443961/149‑1 (MQ=25)
                                                                                                                                      aTGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGa                <  1:1167318/149‑1 (MQ=35)
                                                                                                                                        gcgcCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACGGAACTGTAGGTCGGATa              >  1:1295299/1‑149 (MQ=21)
                                                                                                                                        gcgcCGCTGACGGGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAATTGTCGGTCGGATa              <  2:795524/149‑1 (MQ=12)
                                                                                                                                          gccgcTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAg            >  1:821380/1‑149 (MQ=35)
                                                                                                                                           ccgcTGACGTGACTTAGCCGGCCGACGGGGTGCACTAAACTGTAGGTCGGATAAGACGCATAGCGTCGCATTCGACAACCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGa           <  2:1249613/149‑1 (MQ=11)
                                                                                                                                                  cGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAGGCGACGTATCAGGCCGACGGGGTGCGCTGAACTGTAGGTCGAATACGACGCATAg    >  2:1142889/1‑149 (MQ=11)
                                                                                                                                                    tGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCg  >  2:624120/1‑149 (MQ=32)
                                                                                                                                                    |                                                                                                                                                    
ATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCG  >  NC_000913/349569‑349865

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