Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,402,924 T→C 5.0% H111R (CAC→CGC)  abgB ← p‑aminobenzoyl‑glutamate hydrolase, B subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,402,9240TC 5.0% 248.7 / 2.2 79H111R (CAC→CGC) abgBp‑aminobenzoyl‑glutamate hydrolase, B subunit
Reads supporting (aligned to +/- strand):  ref base T (37/38);  new base C (2/2);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCCCCTCGCGAACCATGAACGTTTTACCCGAGCCGCCTTCTTCGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGG  >  NC_000913/1402777‑1403072
                                                                                                                                                   |                                                                                                                                                    
aCCCCCTCGCGAACCATGAACGTTTTACCCGAGCCGCCTTCTTCGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATtgtg                                                                                                                                                     >  2:10477/1‑149 (MQ=255)
   cccTCGCGAACCATGAACGTTTTACCCGAGCCGCCTTCTTCGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTcc                                                                                                                                                  <  2:347766/149‑1 (MQ=255)
   cccTCGCGAACCATGAACGTTTTACCCGAGCCGCCTTCTTCGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTcc                                                                                                                                                  >  4:32192/1‑149 (MQ=255)
                  aaCGTTTTACCCGAGCCGCCTTCTTCGACAGGACAACCATAAAAGCGCACCGTGCCGCCTTGACCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCACAGCAAATAGTGTCAGAAACCGTGACCAta                                                                                                                                   >  6:58842/1‑148 (MQ=255)
                     gTTTTACCCGAGCCGCCTTATTCGCCAGGACAACCATTAAATCGCACCTTGCCGCGTTGCTCATATTGTTACAGCCATTTCTTGACGTCTATTGCAGCGGCAAAGGAGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTc                                                                                                                                <  4:352691/149‑1 (MQ=255)
                          aCCCGAGCCGCCTTCTTCGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCgg                                                                                                                           >  4:407324/1‑149 (MQ=255)
                            ccGAGCCGCCTTCTTCGCCAGGACTACCATAAAAGCTCAACGTGCCGCCTTGCAAATATTGTTCCAGCAATTTCTTGACGGCTATTGCAGCGGCAAAGTCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGTGc                                                                                                                         <  5:58842/149‑1 (MQ=255)
                                   gccTTCTTCGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACg                                                                                                                   >  3:244818/1‑148 (MQ=255)
                                    ccTTCTTCGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTc                                                                                                                                                   <  7:19053/115‑1 (MQ=255)
                                    ccTTCTTCGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTc                                                                                                                                                   >  8:19053/1‑115 (MQ=255)
                                        cttcGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGt                                                                                                              <  3:407324/148‑1 (MQ=255)
                                         ttcGCCAGGACAACCATAAAAGCGCACCGTACCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAg                                                                                                            >  8:342617/1‑149 (MQ=255)
                                          tcGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAgg                                                                                                           >  8:112499/1‑149 (MQ=255)
                                              cAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTgc                                                                                                       >  6:407658/1‑149 (MQ=255)
                                                           aaaGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCtgctt                                                                                          <  8:50502/149‑1 (MQ=255)
                                                           aaaGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGGTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCtaca                                                                                           <  3:221323/148‑4 (MQ=255)
                                                            aaGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCtgcttg                                                                                         >  7:36513/1‑149 (MQ=255)
                                                              gcgcACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCtgcttgct                                                                                       <  6:387166/149‑1 (MQ=255)
                                                                  aCCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAAc                                                                                                >  3:251929/1‑136 (MQ=255)
                                                                  aCCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAAc                                                                                                <  4:251929/136‑1 (MQ=255)
                                                                         cgccTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAAccg                                                                            >  2:350983/1‑148 (MQ=255)
                                                                           ccTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGc                                                                          <  8:36513/149‑1 (MQ=255)
                                                                            cTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGc                                                                          >  1:122306/1‑148 (MQ=255)
                                                                               gCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAg                                                                       <  3:32192/148‑1 (MQ=255)
                                                                                cccATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAggg                                                                     >  5:41423/1‑149 (MQ=255)
                                                                                  cATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTAGCTGACTTAAACCTGCCAGGGCg                                                                   >  1:43513/1‑149 (MQ=255)
                                                                                     atTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCa                                                                >  6:204674/1‑149 (MQ=255)
                                                                                      ttGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCa                                                                >  8:410803/1‑148 (MQ=255)
                                                                                        gTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCATCata                                                              <  1:10477/148‑1 (MQ=255)
                                                                                         ttCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCatat                                                             >  8:331159/1‑148 (MQ=255)
                                                                                         ttCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATAtt                                                            <  6:41423/149‑1 (MQ=255)
                                                                                         ttCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGAAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGAAACCTGCTTGCTGACTTATACCTGCCAGGGCGTCatat                                                             >  5:7695/1‑148 (MQ=255)
                                                                                            cAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATtctc                                                         <  7:331159/149‑1 (MQ=255)
                                                                                            cAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATtctc                                                         <  4:221589/149‑1 (MQ=255)
                                                                                             aGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCgg                                                                                                                           >  1:287065/1‑82 (MQ=255)
                                                                                             aGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCgg                                                                                                                           <  2:287065/82‑1 (MQ=255)
                                                                                              gCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTccc                                                       <  5:407658/149‑1 (MQ=255)
                                                                                              gCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTccc                                                       <  4:158822/149‑1 (MQ=255)
                                                                                               ccATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTcc                                                        >  7:162846/1‑147 (MQ=255)
                                                                                                cATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCcag                                                     <  2:447347/149‑1 (MQ=255)
                                                                                                cATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGCGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGAt                                                                                                                >  1:394875/1‑90 (MQ=255)
                                                                                                cATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGCGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGAt                                                                                                                >  1:401205/1‑90 (MQ=255)
                                                                                                cATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGCGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGAt                                                                                                                <  2:401205/90‑1 (MQ=255)
                                                                                                cATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGCGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGAt                                                                                                                <  2:394875/90‑1 (MQ=255)
                                                                                                    tCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAgg                                                 >  1:334090/1‑149 (MQ=255)
                                                                                                     cTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAggg                                                >  5:161037/1‑149 (MQ=255)
                                                                                                       tGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGc                                               >  8:145512/1‑148 (MQ=255)
                                                                                                        gACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCTTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGc                                               <  8:222394/147‑1 (MQ=255)
                                                                                                        gACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGa                                             >  4:298554/1‑149 (MQ=255)
                                                                                                         aCGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATGTTCACCGGGCGTCAAGGATGTAGGTTGCGCGCAACCTGCTTGCTCACTTCAACCTGCCAGGGCGTCATATTCTCCcagcagtcc                                               >  5:322983/1‑143 (MQ=255)
                                                                                                           ggCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTaaa                                                                               <  2:124429/112‑1 (MQ=255)
                                                                                                           ggCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTaaa                                                                               >  1:124429/1‑112 (MQ=255)
                                                                                                           ggCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATtctc                                                         >  7:243384/1‑134 (MQ=255)
                                                                                                           ggCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATtctc                                                         <  8:243384/134‑1 (MQ=255)
                                                                                                            gCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAAc                                         <  7:410803/149‑1 (MQ=255)
                                                                                                                  gCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGtt                                    <  3:298554/148‑1 (MQ=255)
                                                                                                                   cAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTg                                  >  2:336775/1‑149 (MQ=255)
                                                                                                                   cAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTg                                  >  7:42803/1‑149 (MQ=255)
                                                                                                                    aGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCCCGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAACCCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGc                                 >  4:101145/1‑149 (MQ=255)
                                                                                                                      cGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCt                               <  7:112499/149‑1 (MQ=255)
                                                                                                                        gCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCat                                                               >  2:215150/1‑115 (MQ=255)
                                                                                                                        gCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCat                                                               <  1:215150/115‑1 (MQ=255)
                                                                                                                          aaaGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGAc                           <  5:204674/149‑1 (MQ=255)
                                                                                                                              gcggcgGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGttt                                   <  4:253010/137‑1 (MQ=255)
                                                                                                                              gcggcgGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGttt                                   >  3:253010/1‑137 (MQ=255)
                                                                                                                              gcggcgGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCaaa                       <  4:244818/149‑1 (MQ=255)
                                                                                                                                  cggTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGaa                   <  1:336775/149‑1 (MQ=255)
                                                                                                                                  cggTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGaa                   <  6:161037/149‑1 (MQ=255)
                                                                                                                                        ccAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATa                                           >  8:141048/1‑119 (MQ=255)
                                                                                                                                         cAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATa             >  7:98369/1‑148 (MQ=255)
                                                                                                                                          aGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGa                            <  8:93231/132‑1 (MQ=255)
                                                                                                                                          aGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGa                            >  7:93231/1‑132 (MQ=255)
                                                                                                                                          aGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATaaa           <  1:350983/149‑1 (MQ=255)
                                                                                                                                           gCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTcc                                                        <  4:223953/103‑1 (MQ=255)
                                                                                                                                           gCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTcc                                                        >  3:223953/1‑103 (MQ=255)
                                                                                                                                           gCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAg          <  2:334090/149‑1 (MQ=255)
                                                                                                                                           gCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAg          <  2:122306/149‑1 (MQ=255)
                                                                                                                                            cAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAgg         <  8:162846/149‑1 (MQ=255)
                                                                                                                                                ttGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCtgcttgc                                                                                        >  2:145093/1‑66 (MQ=255)
                                                                                                                                                ttGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCtgcttgc                                                                                        <  1:145093/66‑1 (MQ=255)
                                                                                                                                                   tgtCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACtt                                                                                  <  3:31524/69‑1 (MQ=255)
                                                                                                                                                   tgtCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACtt                                                                                  >  4:31524/1‑69 (MQ=255)
                                                                                                                                                   tgtCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTgg  >  2:152701/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACCCCCTCGCGAACCATGAACGTTTTACCCGAGCCGCCTTCTTCGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGG  >  NC_000913/1402777‑1403072

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: