Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,459,658 T→C 5.4% N271S (AAC→AGC)  fadI ← beta‑ketoacyl‑CoA thiolase, anaerobic, subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,459,6580TC 5.4% 228.7 / 2.8 74N271S (AAC→AGC) fadIbeta‑ketoacyl‑CoA thiolase, anaerobic, subunit
Reads supporting (aligned to +/- strand):  ref base T (33/37);  new base C (2/2);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

GCTGGACCGAGCAACATGTCCTGCCAGACATCAATCGCAGTAAATGCGTAGCTGCGCAGATACCCCAGCGGCACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCTCTTCTTTGAG  >  NC_000913/2459511‑2459806
                                                                                                                                                   |                                                                                                                                                    
gCTGGACCGAGCAACATGTCCTGCCAGACATCAATCGCAGTAAATGCGTAGCTGCGCAGATACCCCAGCGCCACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATAACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGtt                                                                                                                                                     >  8:50372/1‑149 (MQ=255)
         aGCAACATGTCCTGCCAGACATCAATCGCAGTAAATGCGTAGCTGCGCAGATACCCCAGCGGCACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGt                                                                                                                                            >  3:228316/1‑149 (MQ=255)
             aCATGTCCTGCCAGACATCAATCGCAGTAAATGCGTAGCTGCGCAGATACCCCAGCGGCACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTAcc                                                                                                                                        >  5:89660/1‑149 (MQ=255)
                tGTCCTGCCAGACATCAATCGCAGTAAATGCGTAGCTGCGCAGATACCCCAGCGGCACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGtt                                                                                                                                     >  3:171401/1‑149 (MQ=255)
                 gTCCTGCCAGACATCAATCGCAGTAAATGCGTAGCTGCGCAGATACACCAGCGGCACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTc                                                                                                                                    >  1:26930/1‑149 (MQ=255)
                  tCCTGCCAGACATCAATCGCAGTAAATGCGTAGCTGCGCAGATACCCCAGCGGCACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTcc                                                                                                                                   <  6:31207/149‑1 (MQ=255)
                       ccAGACATCAATCGCAGTAAATGCGTAGCTGCGCAGATACCCCAGCGGCCCCAGCCCTAATTCTTTCGCCCGGGATTAAGTCATCAGGATCACCGCTGACGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACACTTCCGTGtt                                                                                                                              >  5:179535/1‑149 (MQ=255)
                       ccAGACATCAATCGCAGTAAATGCGTAGCTGCGCAGATACCCCAGCGGCACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCgtgt                                                                                                                               <  6:89660/148‑1 (MQ=255)
                       ccAGACATCAATCGCAGTAAATGCGTAGCTGCGCAGATACCCCAGCGGCACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGtt                                                                                                                              <  6:57776/149‑1 (MQ=255)
                        cAGACATCAATCGCAGTAAATGCGTAGCTGCGCAGATACCCCAGCGGCACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGttt                                                                                                                             <  3:107566/149‑1 (MQ=255)
                         aGACATCAATCGCAGTAAATGCGTAGCTGCGCAGATACCCCAGCGGCACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTg                                                                                                                            >  3:170325/1‑149 (MQ=255)
                             aTCAATCGCAGTAAATGCGTAGCTGCGCAGATACCCCAGCGGCACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGt                                                                                                                                            >  8:139040/1‑129 (MQ=255)
                             aTCAATCGCAGTAAATGCGTAGCTGCGCAGATACCCCAGCGGCACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGt                                                                                                                                            <  7:139040/129‑1 (MQ=255)
                                       gTAAATGCGTAGCTGCGCAGATACCCCAGCGGCACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCCTCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCgg                                                                                                              >  2:110618/1‑149 (MQ=255)
                                                   cTGCGCAGATACCCCAGCGGCACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGc                                                                                                  <  5:128417/149‑1 (MQ=255)
                                                   cTGCGCAGATACCCCAGCGGCACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGc                                                                                                  <  2:26930/149‑1 (MQ=255)
                                                      cgcAGATACCCCAGCGGCACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCgg                                                                                                              <  2:225163/134‑1 (MQ=255)
                                                      cgcAGATACCCCAGCGGCACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCgg                                                                                                              >  1:225163/1‑134 (MQ=255)
                                                       gcAGATACCCCAGCGGCACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTaa                                                                                              >  3:22901/1‑149 (MQ=255)
                                                        cAGATACCCCAGCGGCACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAAt                                                                                             <  4:171401/149‑1 (MQ=255)
                                                          gATACCCCAGCGGCACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCg                                                                                           <  2:161710/149‑1 (MQ=255)
                                                                ccAGCGGCACCAGCCCTAATTCTTTCGCCCGGGAGTCAGTCATCATGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGc                                                                                     <  6:40166/149‑1 (MQ=255)
                                                                  agaggCACCAGCCCTAATTCTTTCGACCGGGATTCAGTCATCAGGAGCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGGTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGa                                                                                   <  6:189427/146‑1 (MQ=255)
                                                                   gCGGCACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGTGAg                                                                                  <  8:178604/149‑1 (MQ=255)
                                                                   gCGGCACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAg                                                                                  <  8:166286/149‑1 (MQ=255)
                                                                       cACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAAt                                                                              >  7:210667/1‑149 (MQ=255)
                                                                          cAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTAc                                                                           >  3:44126/1‑149 (MQ=255)
                                                                           aGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCGAGCGAGGAATTAcc                                                                          <  4:228316/149‑1 (MQ=255)
                                                                             cccTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACcgc                                                                        <  3:234069/149‑1 (MQ=255)
                                                                               cTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACcgcg                                                                       >  7:187133/1‑148 (MQ=255)
                                                                               cTAATTCTTTCGCCCGGGATTCAGTCATCAGCATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTACGTGTTTGCGATCAAACTCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGTAATTACcgcg                                                                       >  6:223514/1‑148 (MQ=255)
                                                                                tAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAAt                                                                                             >  2:159848/1‑125 (MQ=255)
                                                                                tAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAAt                                                                                             <  1:159848/125‑1 (MQ=255)
                                                                                tAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGaa                                                                     >  3:60610/1‑149 (MQ=255)
                                                                                  aTTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAat                                                                    <  5:103619/148‑1 (MQ=255)
                                                                                      tttCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGACCCCGTTACCGTTCAGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATAttgt                                                               >  5:25154/1‑149 (MQ=255)
                                                                                        tCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATCACCGCGAATAttg                                                                <  7:225281/146‑1 (MQ=255)
                                                                                          gCCCGGGATTCAGTCATCACGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGACGCCGTTACCGTTCCGTGTTAGCGATCAAACGCCGGGCGCAGCTATGCGTAATCGGCAAGCGAGCAATTACCGCGAAAATAGTAGTc                                                           >  8:83014/1‑149 (MQ=255)
                                                                                            ccGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCtt                                                         >  1:184681/1‑149 (MQ=255)
                                                                                             cGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGc                                                                                                  >  8:148148/1‑107 (MQ=255)
                                                                                             cGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGc                                                                                                  <  7:148148/107‑1 (MQ=255)
                                                                                               ggATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGa                                                      <  4:44126/149‑1 (MQ=255)
                                                                                                gATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGa                                                      >  8:204636/1‑148 (MQ=255)
                                                                                                     aGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCg                                                <  3:140251/149‑1 (MQ=255)
                                                                                                      gTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCgg                                               >  5:45539/1‑149 (MQ=255)
                                                                                                       tcatcaGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGgt                                              >  7:81067/1‑149 (MQ=255)
                                                                                                         atcaGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGgttg                                            >  3:211363/1‑149 (MQ=255)
                                                                                                          tcaGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGgttgt                                           >  7:97930/1‑149 (MQ=255)
                                                                                                           caGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGCTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGgttgtt                                          >  5:216023/1‑149 (MQ=255)
                                                                                                           caGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGCTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGgttgtt                                          >  5:219680/1‑149 (MQ=255)
                                                                                                            aGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGgttgtt                                          <  4:170325/148‑1 (MQ=255)
                                                                                                            aGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGgttgtt                                          <  4:211363/148‑1 (MQ=255)
                                                                                                              gATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGc                                                                                                  <  8:217319/90‑1 (MQ=255)
                                                                                                              gATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGc                                                                                                  >  7:217319/1‑90 (MQ=255)
                                                                                                               aTCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTtata                                      <  8:97930/149‑1 (MQ=255)
                                                                                                               aTCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTtata                                      <  4:22901/149‑1 (MQ=255)
                                                                                                               aTCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTtat                                       >  6:55298/1‑148 (MQ=255)
                                                                                                                tCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTtat                                       >  4:168492/1‑147 (MQ=255)
                                                                                                                tCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATaa                                     <  5:55298/149‑1 (MQ=255)
                                                                                                                  aCCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGc                                                 <  1:220005/135‑1 (MQ=255)
                                                                                                                  aCCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGc                                                 >  2:220005/1‑135 (MQ=255)
                                                                                                                  aCCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAag                                    <  8:210667/148‑1 (MQ=255)
                                                                                                                   ccGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAagga                                  <  8:32588/149‑1 (MQ=255)
                                                                                                                    cGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGTGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAaggag                                 <  6:179535/149‑1 (MQ=255)
                                                                                                                     gCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAatgatg                                <  6:238610/149‑3 (MQ=255)
                                                                                                                     gCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAaggagg                                <  7:204636/149‑1 (MQ=255)
                                                                                                                       tGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAAGGAGGGa                              >  2:15812/1‑149 (MQ=255)
                                                                                                                        gCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGGTTGttt                                         <  5:487/137‑1 (MQ=255)
                                                                                                                        gCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGGTTGttt                                         >  6:487/1‑137 (MQ=255)
                                                                                                                               cATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAAGGAGGGATAAAGGCa                      <  2:184681/149‑1 (MQ=255)
                                                                                                                                 tCGGTCAGCGGCGTACTGCTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGt                    <  6:216023/149‑1 (MQ=255)
                                                                                                                                 tCGGTCAGCGGCGTACTGCTTGCCGCCGGTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGt                    <  6:219680/149‑1 (MQ=255)
                                                                                                                                       aGCGGCGTACTGTTTGCCGCCGTTACCGTTGCGTGTTTGCGATCAATCGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGAGAAGCTGTTGTTTATAAGGAGGGATAAAGGCAGTCATCAc              <  5:223514/149‑1 (MQ=255)
                                                                                                                                           gcgTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACctct          >  7:178359/1‑149 (MQ=255)
                                                                                                                                             gTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCtcttc        >  7:80094/1‑149 (MQ=255)
                                                                                                                                                   tttGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCTCTTCTTTGAg  >  1:119287/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCTGGACCGAGCAACATGTCCTGCCAGACATCAATCGCAGTAAATGCGTAGCTGCGCAGATACCCCAGCGGCACCAGCCCTAATTCTTTCGCCCGGGATTCAGTCATCAGGATCACCGCTGCCGCGCCATCGGTCAGCGGCGTACTGTTTGCCGCCGTTACCGTTCCGTGTTTGCGATCAAACGCCGGGCGCAGCTTTGCGTAATCGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCTCTTCTTTGAG  >  NC_000913/2459511‑2459806

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