Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,007,473 G→A 5.1% intergenic (+21/‑37) ygeW → / → ygeX putative carbamoyltransferase/2,3‑diaminopropionate ammonia lyase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,007,4730GA 5.1% 457.1 / 8.6 137intergenic (+21/‑37)ygeW/ygeXputative carbamoyltransferase/2,3‑diaminopropionate ammonia lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (64/66);  new base A (4/3);  total (68/69)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 137 total aligned reads displayed.

CTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCA  >  NC_000913/3007328‑3007619
                                                                                                                                                 |                                                                                                                                                  
ctacctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCg                                                                                                                                                 >  1:421174/1‑149 (MQ=255)
 tacctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCgg                                                                                                                                                >  3:706221/1‑149 (MQ=255)
   cctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGcc                                                                                                                                              <  2:760996/149‑1 (MQ=255)
   cctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGACCAGCTGATCGGcc                                                                                                                                              <  4:238851/149‑1 (MQ=255)
      acAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTc                                                                                                                                                  <  4:515083/142‑1 (MQ=255)
      acAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTc                                                                                                                                                  >  3:515083/1‑142 (MQ=255)
       cAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCctct                                                                                                                                          >  3:693834/1‑149 (MQ=255)
          agaagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAGTACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTAtt                                                                                                                                       <  3:310130/149‑1 (MQ=255)
           gaagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTAttt                                                                                                                                      <  2:527848/149‑1 (MQ=255)
             agCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTcc                                                                                                                                    <  3:175978/149‑1 (MQ=255)
               ccAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCag                                                                                                                                  >  2:222115/1‑149 (MQ=255)
                cAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCaga                                                                                                                                 >  3:202887/1‑149 (MQ=255)
                cAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCaga                                                                                                                                 <  1:709504/149‑1 (MQ=255)
                 aGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCagag                                                                                                                                >  3:596602/1‑149 (MQ=255)
                    tGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGcc                                                                                                                             <  1:282253/149‑1 (MQ=255)
                        aGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCaaaa                                                                                                                         <  1:559761/149‑1 (MQ=255)
                        aGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCaaaa                                                                                                                         >  3:802496/1‑149 (MQ=255)
                         gCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAag                                                                                                                        <  2:408648/149‑1 (MQ=255)
                          ccTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAagg                                                                                                                       <  2:376901/149‑1 (MQ=255)
                           cTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAagga                                                                                                                      <  4:1724/149‑1 (MQ=255)
                            ttATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAATCCCTCGGGCCGCTCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggat                                                                                                                     >  1:931126/1‑149 (MQ=255)
                            ttATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCaaaa                                                                                                                         <  4:239837/145‑1 (MQ=255)
                            ttATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggat                                                                                                                     >  2:56964/1‑149 (MQ=255)
                             tataTCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggata                                                                                                                    <  1:819542/149‑1 (MQ=255)
                             tataTCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATCTCCAGAGGCCAAAaggata                                                                                                                    <  4:353134/149‑1 (MQ=255)
                              atatCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggatag                                                                                                                   <  4:826225/149‑1 (MQ=255)
                              atatCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggatag                                                                                                                   <  4:693834/149‑1 (MQ=255)
                                atcatcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggatagga                                                                                                                 <  4:442977/149‑1 (MQ=255)
                                atcatcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggatagga                                                                                                                 <  2:544026/149‑1 (MQ=255)
                                 tcatcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggataggat                                                                                                                >  1:188889/1‑149 (MQ=255)
                                  catcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAggataggata                                                                                                               <  3:593358/149‑1 (MQ=255)
                                   atcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGAGCCGAACAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGatat                                                                                                              >  3:24673/1‑149 (MQ=255)
                                    tcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATg                                                                                                             >  4:40254/1‑149 (MQ=255)
                                         gcgATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGt                                                                                                        <  4:315047/149‑1 (MQ=255)
                                         gcgATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGt                                                                                                        >  2:661951/1‑149 (MQ=255)
                                           gatgatCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGttt                                                                                                      >  1:251908/1‑149 (MQ=255)
                                           gatgatCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGttt                                                                                                      <  4:197555/149‑1 (MQ=255)
                                               atCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTtctc                                                                                                  >  4:329771/1‑149 (MQ=255)
                                                tCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCaaaa                                                                                                                         <  3:145830/125‑1 (MQ=255)
                                                tCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCaaaa                                                                                                                         >  4:145830/1‑125 (MQ=255)
                                                tCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCa                                                                                                 >  2:888955/1‑149 (MQ=255)
                                                  cTGTCCCGTAAATGCGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGACCAGCTGATCGGcc                                                                                                                                              <  3:563020/102‑1 (MQ=255)
                                                  cTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCAtt                                                                                               >  1:269396/1‑149 (MQ=255)
                                                  cTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGACCAGCTGATCGGccc                                                                                                                                             <  1:170284/103‑1 (MQ=255)
                                                  cTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGACCAGCTGATCGGccc                                                                                                                                             >  2:170284/1‑103 (MQ=255)
                                                  cTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGACCAGCTGATCGGcc                                                                                                                                              >  4:563020/1‑102 (MQ=255)
                                                     tCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCaaaa                                                                                                                         >  1:200903/1‑120 (MQ=255)
                                                     tCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGaa                                                                                            <  2:431318/149‑1 (MQ=255)
                                                       ccGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggatagga                                                                                                                 <  2:412921/126‑1 (MQ=255)
                                                        cGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGAt                                                                                         <  4:802663/149‑1 (MQ=255)
                                                         gTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGAtt                                                                                        >  1:150917/1‑149 (MQ=255)
                                                            aattcaccAAACCAGGTGCACTGCTCGCGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGat                                                                                     <  4:596602/143‑1 (MQ=255)
                                                            aaTACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTTAAGATTGat                                                                                     >  1:837768/1‑149 (MQ=255)
                                                            aaTACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGat                                                                                     <  4:202887/149‑1 (MQ=255)
                                                            aaTACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGat                                                                                     <  4:609346/149‑1 (MQ=255)
                                                            aaTAAGACAACCCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCGAGCACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCAAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGat                                                                                     <  1:843235/149‑1 (MQ=255)
                                                                    aaaCCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGa                                                                             >  3:886223/1‑149 (MQ=255)
                                                                    aaaCCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGa                                                                             <  2:484489/149‑1 (MQ=255)
                                                                                 tGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTc                                                                <  4:452413/149‑1 (MQ=255)
                                                                                  gCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCa                                                                                                                                   >  3:369727/1‑81 (MQ=255)
                                                                                  gCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATaa                                                                          <  3:902952/138‑1 (MQ=255)
                                                                                  gCTCGAGCAACTGCGGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATaa                                                                          >  4:902952/1‑138 (MQ=255)
                                                                                   cTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCaa                                                              <  2:800840/149‑1 (MQ=255)
                                                                                    tCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAAc                                                             <  3:188844/149‑1 (MQ=255)
                                                                                           aCTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGaaa                                                      <  2:837768/149‑1 (MQ=255)
                                                                                              gctgAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACAt                                                   <  2:251908/149‑1 (MQ=255)
                                                                                               ctgAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATc                                                  >  4:812327/1‑149 (MQ=255)
                                                                                                tggaagaagCGCAAGAACGCGAGAATTAACCCCTCGGGCCGATCAGCTGGTCGGCTCTCTATTTCCAGAGGTCAAAAGGTTAGGATATGTCCGTTTTCTCATTGACGATTGAGATCGCCGATAGCAAATTTTTCAACGGCGAAACATCa                                                 >  4:574460/4‑149 (MQ=255)
                                                                                                tgAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGa                                                                                      <  3:283956/112‑1 (MQ=255)
                                                                                                tgAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCa                                                 <  2:150917/149‑1 (MQ=255)
                                                                                                  aaaGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCAcc                                               <  3:217542/149‑1 (MQ=255)
                                                                                                  aaaGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCAcc                                               >  4:52764/1‑149 (MQ=255)
                                                                                                  aaaGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCAcc                                               >  3:709440/1‑149 (MQ=255)
                                                                                                    agaagCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGc                                             >  2:484012/1‑149 (MQ=255)
                                                                                                    agaagCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATAGCTCCGATAACAAATTTTTCAACGGCGAAACATCACCGc                                             >  4:76805/1‑149 (MQ=255)
                                                                                                      aagCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTtct                                                                                                   >  4:136995/1‑93 (MQ=255)
                                                                                                      aagCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTtct                                                                                                   <  3:136995/93‑1 (MQ=255)
                                                                                                        gcgcAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAAc                                                             >  1:655604/1‑129 (MQ=255)
                                                                                                        gcgcAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAAc                                                             <  2:655604/129‑1 (MQ=255)
                                                                                                             aGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGc                                    <  4:782513/149‑1 (MQ=255)
                                                                                                                aCGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTc                                                                <  3:392472/118‑1 (MQ=255)
                                                                                                                aCGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTc                                                                >  4:392472/1‑118 (MQ=255)
                                                                                                                aCGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCaaa                                 <  2:188889/149‑1 (MQ=255)
                                                                                                                  gcgTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGc                               >  1:101255/1‑149 (MQ=255)
                                                                                                                     tGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAg                            >  2:703162/1‑149 (MQ=255)
                                                                                                                      gAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAgg                           <  3:52764/149‑1 (MQ=255)
                                                                                                                            aCCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAAc                     <  3:812327/149‑1 (MQ=255)
                                                                                                                              cccTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTg                   >  4:387687/1‑149 (MQ=255)
                                                                                                                                cgcGGGCCTATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGgc                 >  4:279756/3‑149 (MQ=255)
                                                                                                                                cTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGgc                 >  3:793106/1‑149 (MQ=255)
                                                                                                                                  cGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGgcgc               <  1:888955/149‑1 (MQ=255)
                                                                                                                                  agggCCGACCAGCTGATCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGgcgc               >  1:493855/2‑149 (MQ=255)
                                                                                                                                      ccGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCa           <  4:706221/149‑1 (MQ=255)
                                                                                                                                        gATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGt         >  1:434765/1‑149 (MQ=255)
                                                                                                                                           cAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGtttg      <  3:329771/149‑3 (MQ=255)
                                                                                                                                           cAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTcc      <  1:703162/149‑1 (MQ=255)
                                                                                                                                           cAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTcc      >  1:291212/1‑149 (MQ=255)
                                                                                                                                            aGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTcca     >  4:497850/1‑149 (MQ=255)
                                                                                                                                            aGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTATCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGGAAGGCCAGGCCAAACTGGCGCGCCAGTTcca     <  4:567211/149‑1 (MQ=255)
                                                                                                                                               tGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTccacca  <  1:56964/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCA  >  NC_000913/3007328‑3007619

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