Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,356,446 A→T 7.1% intergenic (‑34/+25) dtpC ← / ← cadA dipeptide and tripeptide permease/lysine decarboxylase, acid‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,356,4460AT 7.1% 226.1 / 4.8 81intergenic (‑34/+25)dtpC/cadAdipeptide and tripeptide permease/lysine decarboxylase, acid‑inducible
Reads supporting (aligned to +/- strand):  ref base A (39/35);  new base T (2/4);  total (41/40)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAG  >  NC_000913/4356299‑4356593
                                                                                                                                                   |                                                                                                                                                   
cTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCaa                                                                                                                                                    <  1:1037952/149‑1 (MQ=255)
  gATGGGTGAGATAGAGAATGAGTAAGGCACGGATGCCGTAAAAACTGAAGGGCTCCCAGATTTGGATCGCCACGATAGAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGc                                                                                                                                                  <  1:1018691/149‑1 (MQ=255)
       gtgtGATTGAGAATGAGTGAGGCACGCATGCCGTAAAACCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGTGCAAAAAAGGGAAGTGGCTTGCCACtt                                                                                                                                             <  2:593044/149‑1 (MQ=255)
       gTGAGTTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCCGTGTGGGCGTGTTCATGTATTCTCCTTATGAGCAATAAAGGGAAGTGGCTAGcgcctt                                                                                                                                             >  2:973034/1‑144 (MQ=255)
       gTGAGATAGAGAACGAGTACGGCACGCATTCCGTAAAAACTGAAGGACTCCCATATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGGGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACtt                                                                                                                                             <  2:592464/149‑1 (MQ=255)
         gagaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTcc                                                                                                                                           <  1:324560/149‑1 (MQ=255)
         gagaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTCTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTTGCAAGCCACTTcc                                                                                                                                           <  1:1044522/149‑1 (MQ=255)
          agaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTTTGAGGGTGTTTTCATGTGTTCTCCTTATTAGCAAAAGCGGGAAGTGGCAGGcccagcccc                                                                                                                                          >  1:574214/1‑142 (MQ=255)
              ggaggaggAGGAAGGCGCGCGTGCCGGAAAAACTGAAGTACTCGCAGATTGGGAGCGCCACGATATAGTATATCGCGCGCGGCTGTGCGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGt                                                                                                                                      <  2:38316/142‑1 (MQ=255)
              agagAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCCCGCGGCCTTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGt                                                                                                                                      >  1:1091944/1‑149 (MQ=255)
               gagaATGAGTCAGGCACGCATGCCACAAAAACTGAAGTACTCTCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTa                                                                                                                                     <  2:1437014/149‑1 (MQ=255)
               gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGAGCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGGGTTTTCATGTGTTCTCCCTATGGGCAAAGAAGGGAAGTGGCAAGCCCCTTCCCTTGtc                                                                                                                                     >  2:787842/1‑148 (MQ=255)
                 gaATGAGTAAGGCACGCATGCCGTAAAAATTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACg                                                                                                                                   <  1:1353795/149‑1 (MQ=255)
                    tGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGGTTTCCTGTGTTCTCCTTATGAGCAAAAAAGGGACGTGGCAAGCAACCTCCCTTGCACGAGc                                                                                                                                >  2:72907/1‑149 (MQ=255)
                       gTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTaa                                                                                                                             >  2:1082698/1‑149 (MQ=255)
                        tAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAACTGGCAAGCCACTTCCCTTGTACGAGCTAat                                                                                                                            >  2:155525/1‑149 (MQ=255)
                              aCGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATTATCAAAATAGGTAAGTGGCAAGCCACTTCCCTTGTGCGCGCTActttcctt                                                                                                                      >  1:797775/1‑144 (MQ=255)
                                     ccGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCCGTGAGGGTGAAAACAAGTGTTCTGCGGATGAGCAAAAAAGGGAGGTGGCATTCCACTTCCCTTGGACGAGCTACTTATGTTCTGCtgt                                                                                                               >  1:794371/1‑147 (MQ=255)
                                      cGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGCGTTTTCATGAGTTCTCCTTATGGGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTGTTGCAttc                                                                                                              >  2:1343459/1‑149 (MQ=255)
                                        taacatcTGAAGTACTCGCAGATTTGGAGCGCCTCCATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCCTATTACCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttctt                                                                                                            <  1:465989/143‑1 (MQ=255)
                                        tAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttctt                                                                                                            <  1:1021224/149‑1 (MQ=255)
                                        tAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttctt                                                                                                            >  2:625793/1‑149 (MQ=255)
                                          aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTtcttct                                                                                                          <  1:1258326/149‑1 (MQ=255)
                                            aaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTGTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGGGCTAATTATTTTTTGCTTTCtttttt                                                                                                        >  2:1080479/1‑149 (MQ=255)
                                            aaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCttt                                                                                                        <  1:771029/149‑1 (MQ=255)
                                               tGAAGTACACCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGGATTGTCCTTATGAGCAAAACAGGGTAGTGTCAAGCCACTTCCCTTGTAAGAGCTAATTATTTTTTGGTTTCTTCTTTCaa                                                                                                     >  1:1064436/1‑149 (MQ=255)
                                                   gTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAcc                                                                                                 <  1:977466/149‑1 (MQ=255)
                                                   gTACTCCCAGATTTGGACCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAcc                                                                                                 <  1:123062/149‑1 (MQ=255)
                                                    tACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGCGTGTTTTCATGTGTGCCCCCCATGAGCAAAAAAGGGAAGAGGCAAGCCACTTCCCTCCTACGCGCTAATTATATTTTGCTGTCTTCTTTccccaccg                                                                                                >  1:1415833/1‑142 (MQ=255)
                                                     aCTTCCAGAGTTGGATCGCCACGATAGAGTATATCGCGCGCGGCTGTGAGGGTGGTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCGAATTATTTTTTGCTTTCTTCTTTCAATACCtt                                                                                               <  1:442876/149‑1 (MQ=255)
                                                       tCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGAGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCATTTCCCTGGTATGGGCTACTTATTTTTTGCTTTCTTAGTTCAAtccattta                                                                                             >  1:672606/1‑142 (MQ=255)
                                                        cccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAAc                                                                                            >  2:1117913/1‑149 (MQ=255)
                                                        cccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGCGTTATCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCGTGTACGAGCTAATTAGTGTTTGCGTGCTTCTTTCAATACCTTAAc                                                                                            >  2:849899/1‑149 (MQ=255)
                                                        cccAGATTTGGATCGCCACGATATAGTACATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAAc                                                                                            <  2:808530/149‑1 (MQ=255)
                                                         ccAGATTTCGAGCGCCACGATATGGTATATCGCGGGCGGCTGTGAGGGTGTTTTCATGTGTTCGCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACg                                                                                           <  1:220858/149‑1 (MQ=255)
                                                             aTTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTCGTACGAGCTAATTATTTTTTGTTTTCTTTTTTAAATTCCTTAACGGtat                                                                                       >  2:996446/1‑149 (MQ=255)
                                                                 ggATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACNAGCTAATTATTTTTTNCAGTCNTCTTTCANTACCTTAACGGTATAGCg                                                                                   >  1:1197765/1‑149 (MQ=255)
                                                                    tCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGcc                                                                                >  1:718046/1‑149 (MQ=255)
                                                                    tCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGcc                                                                                >  1:503653/1‑149 (MQ=255)
                                                                          cGATATAGTATATCGCGGGCGGCTGTGAGTGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGCGCTAATTATTTTTTGCTTTGTTGTGCCTGGACCTTAACGGTATAGCGGCCATCAGc                                                                          <  2:677186/149‑1 (MQ=255)
                                                                           gATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGTTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGcg                                                                         >  2:791136/1‑148 (MQ=255)
                                                                            atatAGCATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCATTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCt                                                                        <  1:848424/149‑1 (MQ=255)
                                                                               taGTATATCGCGCTCGGCTGGGAGGTTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTCCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGAc                                                                     <  2:765409/149‑1 (MQ=255)
                                                                                aGTTTATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACg                                                                    <  2:1233647/149‑1 (MQ=255)
                                                                                aGTATATCGGGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTACCGGTATAGCGGCCATCAGCTTGACg                                                                    <  2:1022545/149‑1 (MQ=255)
                                                                                 gTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACgg                                                                   >  1:1412733/1‑149 (MQ=255)
                                                                                  tataTCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGt                                                                  >  2:1400842/1‑149 (MQ=255)
                                                                                   atatCGCGGGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAACAGGGAAGTGGATATCCATTTTCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTa                                                                 <  2:882378/149‑1 (MQ=255)
                                                                                             cgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGGATGCACCGTGa                                                       >  2:590066/1‑149 (MQ=255)
                                                                                              ggCTGTGATCGTGTTCTCATGTGTTCTCCTTATCAGCAAAAAAGGGCAGTGGCAAGCCCCTTCCCTTGTACGAGCTAATTATTTTTTGCCTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGaa                                                      <  1:547492/149‑1 (MQ=255)
                                                                                                cTGTGAGGTGGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATGTTTTGCTTTCTTCTTTCAATACCTCAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAata                                                    <  2:415042/149‑1 (MQ=255)
                                                                                                   tgAGGGTGTTTTCATGTGTTCCCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCg                                                 >  2:1171374/1‑149 (MQ=255)
                                                                                                     aGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGt                                               <  2:680340/149‑1 (MQ=255)
                                                                                                      gggTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGtt                                              >  1:144602/1‑149 (MQ=255)
                                                                                                       ggTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGttt                                             <  1:615633/149‑1 (MQ=255)
                                                                                                       ggTGTTTTCATGTGATCTACATATGAGGAAAGAATGTGAGTGGGAGGCCACGTCCCTTGTACGAGCTAATTATTGTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGttt                                             >  2:1439241/1‑149 (MQ=255)
                                                                                                        gtgtTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTc                                            >  2:770742/1‑149 (MQ=255)
                                                                                                         tgtTTTCATGTGGTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCa                                           <  1:73458/149‑1 (MQ=255)
                                                                                                          gtTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCGCCGTGAATATCGGTTTCaa                                          <  1:1056238/149‑1 (MQ=255)
                                                                                                           ttttCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCaaa                                         >  1:141308/1‑149 (MQ=255)
                                                                                                                aTGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTTCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCCGCCTGACGGTATGCACCCTGAATATCGGCCTCCAAGCACg                                    >  2:199761/1‑149 (MQ=255)
                                                                                                                aTGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCg                                    >  1:1236540/1‑149 (MQ=255)
                                                                                                                aTGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTCTTGCTGTCTTCTTTCAATACCTTAACGGTACAGCGGCCATAAGCCTGACGGTATGCACCGTGAATATCTGTTTCAAAGCCCg                                    <  1:996446/149‑1 (MQ=255)
                                                                                                                 tgtgTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGCGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAACGGCCATCAGCCTGACGAAATACACCGTGAATATCAGTTTCAAAGCCCgg                                   >  1:448043/1‑149 (MQ=255)
                                                                                                                    gtTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGGACGAGCTAATTATTGTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATa                                <  1:53226/149‑1 (MQ=255)
                                                                                                                     ttCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTACCGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAACACCGGTTTCAAAGCCCGGAAAg                               >  1:376130/1‑149 (MQ=255)
                                                                                                                         ccTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAg                           <  2:1091944/149‑1 (MQ=255)
                                                                                                                           ttATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAgcg                         >  1:1411758/1‑149 (MQ=255)
                                                                                                                            tATGCGCAAAAAAGGGATGTGGCGAGCCCCTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAgcgc                        <  2:264331/149‑1 (MQ=255)
                                                                                                                               gAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTCTTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGa                     <  1:315315/149‑1 (MQ=255)
                                                                                                                                aGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGAt                    >  1:77053/1‑149 (MQ=255)
                                                                                                                                    aaggggggAAGCGGAAAGCCACTTCCCTTGGACGGGCTCATGATTTTTGGCTTTCTTCTTGCAATCCCTTAACGGTATAGCGGCCATCAGCCGGCCGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTca                <  2:634229/144‑1 (MQ=255)
                                                                                                                                      aaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTcaca              >  2:1079883/1‑149 (MQ=255)
                                                                                                                                       aaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTcacac             <  1:162838/149‑1 (MQ=255)
                                                                                                                                        aGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTcacaca            <  2:718046/149‑1 (MQ=255)
                                                                                                                                        aGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTcacaca            <  2:503653/149‑1 (MQ=255)
                                                                                                                                          ggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTTTTTTTTTTTTTCAATACCTTAACGGTATAGCGGCCATCATCCTGACGGTATGCACCGTGAATGTCGTTTTCAAACCGCGGATAGGGTGCGCCGCTTTCGCCCAGc          >  2:70566/1‑149 (MQ=255)
                                                                                                                                           gAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTCAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGGGAGCGCCGATTTAACacaaca         >  1:981259/1‑149 (MQ=255)
                                                                                                                                              gTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCACCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCt      <  1:1117913/149‑1 (MQ=255)
                                                                                                                                                ggCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGc    >  1:560632/1‑149 (MQ=255)
                                                                                                                                                ggCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAATGAGCGCCGATTTCACACACCATCTGc    >  1:331521/1‑149 (MQ=255)
                                                                                                                                                  cAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAg  <  1:1082698/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAG  >  NC_000913/4356299‑4356593

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