Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 253,361 A→G 7.9% intergenic (‑146/+25) PP_0202 ← / ← PP_0203 CBS domain‑containing protein/dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947253,3610AG 7.9% 731.0 / 35.4 241intergenic (‑146/+25)PP_0202/PP_0203CBS domain‑containing protein/dipeptidase
Reads supporting (aligned to +/- strand):  ref base A (106/116);  new base G (5/14);  total (111/130)
Fisher's exact test for biased strand distribution p-value = 9.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.21e-01

Only 100 of 244 total aligned reads displayed.

CATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCGGCATG  >  NC_002947/253213‑253509
                                                                                                                                                    |                                                                                                                                                    
cATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAAc                                                                                                                                                     >  2:2635274/1‑150 (MQ=255)
     gggCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGa                                                                                                                                                >  2:913638/1‑150 (MQ=255)
     gggCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGa                                                                                                                                                >  2:4299966/1‑150 (MQ=255)
     gggCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGCTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGa                                                                                                                                                >  1:3720831/1‑150 (MQ=255)
      ggCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGaa                                                                                                                                               <  2:5003131/150‑1 (MQ=255)
        ccTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGTTTTAGCCGCGaaaa                                                                                                                                             <  2:695544/150‑3 (MQ=255)
         ctctGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTCTTCGCGGTTTTAGCCGCGAaaag                                                                                                                                            <  1:2301987/150‑1 (MQ=255)
           ctGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAACGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGc                                                                                                                                          <  1:1394294/150‑1 (MQ=255)
           ctGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGc                                                                                                                                          <  1:944726/150‑1 (MQ=255)
             gCGGTTTAATAAAGCACTGTCTCGGATTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAAGGCAGGAAACGGATTGGTGGCTTTTGTAAAACGTCTTGGCTGTTCGCTTGTACCGGCCTCTTCGCGGTTTTCGCCGCGAAAAGGCCg                                                                                                                                        <  1:3023225/150‑1 (MQ=255)
             gCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCg                                                                                                                                        <  2:3136327/150‑1 (MQ=255)
                ggTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTa                                                                                                                                     >  1:1837055/1‑150 (MQ=255)
                ggTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTa                                                                                                                                     <  2:1936169/150‑1 (MQ=255)
                    aaaaaaGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTgc                                                                                                                                 <  2:5178924/150‑1 (MQ=255)
                    aaaaaaGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTgc                                                                                                                                 <  1:4168677/150‑1 (MQ=255)
                    aaaaaaGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTgc                                                                                                                                 >  1:89359/1‑150 (MQ=255)
                     aaaaaGCACTGTCTCGGAGTAACTAAAGCATCCGGTACTGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTgcg                                                                                                                                >  2:3166263/1‑150 (MQ=255)
                      aaaaGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTgcgc                                                                                                                               >  1:2369205/1‑150 (MQ=255)
                         aGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATc                                                                                                                            >  1:5097540/1‑150 (MQ=255)
                         aGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTATAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATc                                                                                                                            <  2:543390/150‑1 (MQ=255)
                           cACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCaa                                                                                                                          >  2:226333/1‑150 (MQ=255)
                           cACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCaa                                                                                                                          >  2:3531669/1‑150 (MQ=255)
                            aCTGTCTCGTATTAAATAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCTATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCaaa                                                                                                                         <  1:1079549/150‑1 (MQ=255)
                            aCTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCaaa                                                                                                                         <  1:2392116/150‑1 (MQ=255)
                              tGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGc                                                                                                                       >  1:2087573/1‑150 (MQ=255)
                                 ctcGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTc                                                                                                                    <  2:4571372/150‑1 (MQ=255)
                                 ctcGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTc                                                                                                                    <  1:4215716/150‑1 (MQ=255)
                                 ctcGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTc                                                                                                                    <  2:2335643/150‑1 (MQ=255)
                                  tcGGAGTAACTAAAGCATGCGGTACGGCGCTCTGCACGGCAATGCAGGACGCGGCATGGAGGCGATTGAAAAACGTCATGGCTGTGCGCTTGTTGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGGCGGTACTGCGCATCAAAGCTTCg                                                                                                                   >  2:2531465/1‑150 (MQ=255)
                                  tcGGAGTAACTAAAGCATCCGTTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCg                                                                                                                   <  1:5362562/150‑1 (MQ=255)
                                  tcGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCg                                                                                                                   <  1:5395626/150‑1 (MQ=255)
                                       gTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATc                                                                                                              <  2:2087573/150‑1 (MQ=255)
                                        tAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCt                                                                                                             <  1:4245505/150‑1 (MQ=255)
                                           cTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGcc                                                                                                          >  1:4567159/1‑150 (MQ=255)
                                            tAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGccc                                                                                                         >  2:1449838/1‑150 (MQ=255)
                                              aaGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGa                                                                                                       >  1:95607/1‑150 (MQ=255)
                                                  aTCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAAc                                                                                                                                                     <  1:3448055/100‑1 (MQ=255)
                                                  aTCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCcg                                                                                                   <  1:2231657/150‑1 (MQ=255)
                                                     cGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCgcgc                                                                                                <  2:1561774/150‑1 (MQ=255)
                                                     cGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCgcgc                                                                                                >  2:314978/1‑150 (MQ=255)
                                                      ggTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGcc                                                                                               >  2:1250356/1‑150 (MQ=255)
                                                      ggTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGcc                                                                                               <  2:3256219/150‑1 (MQ=255)
                                                       gTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATAGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGcca                                                                                              >  1:5150291/1‑150 (MQ=255)
                                                            gcgcTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCAtt                                                                                         <  1:5591814/150‑1 (MQ=255)
                                                                 ctGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCt                                                                                    <  2:5244125/150‑1 (MQ=255)
                                                                 ctGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCt                                                                                    <  2:5620711/150‑1 (MQ=255)
                                                                  tGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTc                                                                                   >  2:4210568/1‑150 (MQ=255)
                                                                     aCGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTcacc                                                                                >  1:4979339/1‑150 (MQ=255)
                                                                       ggCAATGCAGGAAACGGCATGGAGGCGAATGAAAAACGACATGGCTGTTCGCTAGTAGCGGCCTATTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCATCGGCATCTGCCCGATCCGCGCCACCATTTCGCTcaccac                                                                              >  2:456445/1‑150 (MQ=255)
                                                                           aTGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGc                                                                          <  2:4040930/150‑1 (MQ=255)
                                                                              cAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAg                                                                       >  2:5488819/1‑150 (MQ=255)
                                                                                    aCGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGc                                                                 <  2:1539426/150‑1 (MQ=255)
                                                                                    aCGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCGTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGc                                                                 <  2:3721540/150‑1 (MQ=255)
                                                                                          tGGAGGCGATTGAAAAACGTCATGGGTTTTCGATTGTAGCGGCCTTTTCGCGGATAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAAc                                                           <  1:484946/150‑1 (MQ=255)
                                                                                          tGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAAc                                                           <  2:4768272/150‑1 (MQ=255)
                                                                                             aGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGc                                                        >  1:3749271/1‑150 (MQ=255)
                                                                                             aGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGc                                                        >  2:2022959/1‑150 (MQ=255)
                                                                                              ggCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCt                                                       <  1:322957/150‑1 (MQ=255)
                                                                                               gCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTc                                                      >  1:146410/1‑150 (MQ=255)
                                                                                                  aTTGTAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCAc                                                   <  1:5488819/150‑1 (MQ=255)
                                                                                                  aTTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCAc                                                   >  2:3699643/1‑150 (MQ=255)
                                                                                                    tGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACgg                                                 >  1:5163532/1‑150 (MQ=255)
                                                                                                     gAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGt                                                <  1:4299966/150‑1 (MQ=255)
                                                                                                     gAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGt                                                >  2:3988050/1‑150 (MQ=255)
                                                                                                      aaaaaCGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTc                                               <  2:5097540/150‑1 (MQ=255)
                                                                                                       aaaaCGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCt                                              >  2:1869393/1‑150 (MQ=255)
                                                                                                       aaaaCGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCt                                              <  1:1195452/150‑1 (MQ=255)
                                                                                                        aaaCGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCtt                                             <  2:3242995/150‑1 (MQ=255)
                                                                                                        aaaCGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCtt                                             <  2:1033328/150‑1 (MQ=255)
                                                                                                        aaaCGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCtt                                             <  1:3376109/150‑1 (MQ=255)
                                                                                                              cATGGCTGTTCGGTTGTAGGGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtc                                       <  1:3916538/150‑1 (MQ=255)
                                                                                                               aTGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtcg                                      <  2:1652196/150‑1 (MQ=255)
                                                                                                               aTGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtcg                                      <  1:3440707/150‑1 (MQ=255)
                                                                                                               aTGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtcg                                      <  1:754510/150‑1 (MQ=255)
                                                                                                                 ggCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtcgtt                                    >  2:4576667/1‑150 (MQ=255)
                                                                                                                     gTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGc                                                        >  1:1001069/1‑126 (MQ=255)
                                                                                                                       tCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtcgttgtcgtt                              <  1:3637260/150‑1 (MQ=255)
                                                                                                                       tCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtcgttgtcgtt                              <  1:3325314/150‑1 (MQ=255)
                                                                                                                         gCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCNGCATGAACTGCTCCACGGTCTTGAACTcgttgtcgttgt                            <  2:49606/150‑1 (MQ=255)
                                                                                                                         gCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTcgttgtcgttgt                            <  2:3881045/150‑1 (MQ=255)
                                                                                                                         gCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTcgttgtcgttgt                            <  2:223387/150‑1 (MQ=255)
                                                                                                                         gCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTcgttgtcgttgt                            >  2:5338252/1‑150 (MQ=255)
                                                                                                                         gCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTcgttgtcgttgt                            >  1:1543336/1‑150 (MQ=255)
                                                                                                                           ttGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGc                          <  1:2883113/150‑1 (MQ=255)
                                                                                                                                gCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGgt                     <  1:5261203/150‑1 (MQ=255)
                                                                                                                                gCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGgt                     <  2:5289452/150‑1 (MQ=255)
                                                                                                                                 cGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGgtg                    >  1:4466831/1‑150 (MQ=255)
                                                                                                                                 cGGCCTTTTCGCGGCTAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGgtgt                   <  2:488386/150‑1 (MQ=255)
                                                                                                                                   gCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTa                  >  1:3975479/1‑150 (MQ=255)
                                                                                                                                   gCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGGTGTCGTTGTGCCCGgtgtc                  >  2:1225635/1‑149 (MQ=255)
                                                                                                                                    ccTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTAc                 >  2:4429721/1‑150 (MQ=255)
                                                                                                                                    ccTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGGTGTCGTTGTGCCCGGTGTAc                 >  2:4429471/1‑150 (MQ=255)
                                                                                                                                     ctcttCGCGGTTTTAGCCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACt                >  1:1311098/4‑150 (MQ=255)
                                                                                                                                      ttttCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACtt               >  1:5315494/1‑150 (MQ=255)
                                                                                                                                      ttttCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACtt               >  2:85006/1‑150 (MQ=255)
                                                                                                                                      tcttCGCGGTTTTAGCCGCGAAAAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACtt               >  2:338695/3‑150 (MQ=255)
                                                                                                                                           gcgGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCt          >  2:1269275/1‑150 (MQ=255)
                                                                                                                                           gcgGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCt          >  1:1990572/1‑150 (MQ=255)
                                                                                                                                                tAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCGGc     >  1:669474/1‑150 (MQ=255)
                                                                                                                                                  aaaCCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCGGCAt   >  2:1741620/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCGGCATG  >  NC_002947/253213‑253509

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